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ABSTRACT 


The  LOGACS  (LOvv-Ci  Accelerometer  Calibration  System)  experiment, 
which  contained  a  miniature  electrostatic  accelerometer  (MESA),  was  placed 
in  a  low-altitude  polar  orbit  on  22  May  1967.  The  experiment  provided 
approximately  100  hr  of  acceleration  data  from  which  the  accelerometer 
scale  factor,  accelerometer  null  bias,  and  atmospheric  drag  on  the  satellite 
were  calculated. 

Many  reports  were  written  during  1967  and  1968  on  the  LOGACS 
experiment  and  on  the  analyses  that  were  performed  on  the  data.  These 
reports  have  been  combined  into  this  document,  which  will  be  the  definitive 
report  on  the  LOGACS  experiment  by  The  Aerospat  *  Corporation. 

The  mission  description,  the  characteristics  of  the  accelerometer, 
the  instrumentation  used  in  the  experiment,  and  a  summary  of  the  results  of 
the  inflight  accelerometer  calibration  are  presented  in  Volume  I.  The  accel¬ 
erometer's  sensing  environment  is  discussed,  as  are  the  data  reduction 
methods  used  to  extract  the  scale  factor  and  bias  of  the  accelerometer  and 
the  methods  of  determining  the  acceleration  due  to  aerodynamic  drag. 

Volume  II  contains  data  plots  of  the  reduced  accelerometer  time 
histories  and  of  the  deduced  atmospheric  densities  over  the  flight.  Peculi¬ 
arities  of  the  density  data  are  examined  in  some  detail  and  extensive  com¬ 
parisons  of  the  data  to  various  model  atmospheres  are  made.  These  anal¬ 
yses  are  the  more  interesting  because  of  the  intense  solar  and  geomagnetic 
activity  during  the  latter  part  of  the  LOGACS  flight,  which  enormously 
affected  the  atmosphere. 

Models  of  atmospheric  density  have  been  developed  over  the  altitudes 
for  which  LOGACS  data  were  available.  Both  midlatitudc  models  and  polar 
and  auroral  models  are  formulated  and  presented. 

An  analysis  was  also  performed  on  the  LOGACS  data  to  determine 
the  wind  magnitude  and  characteristics  normal  to  the  orbit  plane  of  the 
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experiment.  The  results  confirm  that  the  earth's  atmosphere  (up  to  the 
altitude  of  100  nmi)  rotates  with  the  earth  and  find  that  high-velocity  winds 
are  present  in  high-latitude  regions  as  a  result  of  the  great  magnetic  storm. 
An  analysis  of  these  wind  data  shows  they  can  be  described  as  convective 
motion  due  to  excessive  heating  of  the  polar  thermosphere. 
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•  #  •  • 
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6,  Away  from  Perigee  . 

•  •  •  • 
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7,  Toward  Perigee  .  .  . 

•  •  •  • 
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7,  Away  from  Perigee  . 

•  •  •  • 
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9,  Toward  Perigee  .  .  . 

•  •  «  • 
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9,  Away  from  Perigee  . 

•  •  •  • 
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10,  Toward  Perigee  .  . 

•  •  «  • 
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10,  Away  from  Perigee 

•  •  •  * 

B-ll 

11,  Toward  Perigee  .  . 

•  •  •  • 
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•  •  •  • 
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12,  Toward  Perigee  .  . 

•  •  •  • 
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12,  Away  from  Perigee 

•  •  •  • 
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13,  Toward  Perigee  .  . 

•  •  •  • 
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•  •  •  • 
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•  •  •  • 

B-18 
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18,  Away  from  Perigee 

•  •  •  • 
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20,  Away  from  Perigee 

•  •  •  • 

B-21 

21,  Toward  Perigee  .  . 

♦  •  *  • 

B-22 

21,  Away  from  Perigee 

•  •  •  • 

B-23 

22,  Toward  Perigee  .  . 

•  •  •  • 

B  -24 
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SECTION  I 


UPPER  ATMOSPHERE  DENSITY  DETERMINATION  FROM  LOGACS 

A.  THE  ORBIT 
1.  ORBITAL  PARAMETERS 

The  LOGACS  experiment  was  placed  into  a  nearly  polar  orbit  on  an 
Atlas-Agena  vehicle  from  the  Western  Test  Range  at  approximately  1830  GMT 
on  22  May  1967.  The  perigee  altitude  was  80  nmi  initially,  apogee  was 
193  nmi,  and  the  orbital  period  was  very  nearly  90  min,  resulting  in  about 
16  orbital  revs,  per  day.  During  the  eighteenth  rev.  one  day  later,  two 
small  solid  propellant  rockets  on  the  Agena  were  fired  near  perigee  and 
the  orbit  was  adjusted  such  that  its  apogee  altitude  was  raised  to  about 
218  nmi  to  increase  the  orbital  lifetime.  The  orbital  parameters  for 
LOGACS  at  injection  and  after  the  orbit  adjustment  arc  summarized  in 
Table  1. 


Table  1.  Orbital  Parameters  for  LOGACS 


Injection 

After  Orbit  Adjust 
at  Rev.  18.  4 

Date 

22  May  1967 

23  May  1967 

Time 

1839  GMT 

2124  GMT 

Perigee  Altitude 

hp 

80,  0  nmi 

79.  4  nmi 

Apogee  Altitude 

ha 

193,  0  nmi 

217.  8  nmi 

Inclination 

i 

91.5° 

91.  5C 

Eccentricity 

e 

0.015 

0.018 

Period 

P 

89.4  min 

89.  8  min 

Longitude  of  Asc  Node 

4>p 

43.3°  N 

40.  7°  N 

Argument  Perigee 

Wp 

136.7° 

139.3° 
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During  the  course  of  the  flight,  due  to  atmospheric  drag  and 
gravitational  anomalies,  apogee  and  perigee  decayed  to  about  160  and  76  nmi, 
respectively,  and  the  latitude  of  perigee  moved  northward  from  about  41°  N 
to  59°  N  as  shown  in  Fig.  1. 

2.  TRACKING  AND  ORBIT  DETERMINATION 

The  LOGACS  vehicle,  which  had  an  S-band  transponder  on  board, 
was  tracked  by  the  pulse  radar  network  of  the  U.S.  Air  Force  Satellite  Con¬ 
trol  Facility.  Approximately  eight  tracking  passes  per  day  were  obtained 
and  the  orbit  was  determined  from  the  data  by  two  independent  programs 
(Refs.  1  and  2). 

The  quantity  (and  quality)  of  the  tracking  data  used  permitted  solu¬ 
tion  of  an  ephemeris  which  was,  in  most  instances,  within  450  ft  rms  of  the 
radar  range  observations.  This  degree  of  accuracy,  of  course,  is  more 
than  adequate  for  effective  merging  of  the  ephemeris  with  the  accelerations 
measured  by  LOGACS. 

Any  ephemeris  determined  from  tracking  data,  however,  is  always 
more  accurate  in  the  vicinity  of  the  tracking  stations  and  less  accurate  at 
points  far  removed  from  tracking  stations.  In  this  case,  it  is  felt  that  an 
accuracy  of  2500  ft  was  achieved  corresponding  to  regions  in  the  southern 
hemisphere  where  there  was  no  tracking  coverage,  and,  in  a  few  instances 
in  the  northern  hemisphere  when  the  spacecraft  was  adequately  tracked  but 
the  data  were  of  poor  quality.  For  the  atmosphere  density  data  presented 
herein,  the  accuracy  limitations  imposed  by  the  ephemeris  are  not  significant. 

B.  AERODYNAMIC  COEFFICIENTS 

To  deduce  density  values  from  the  LOGACS  acceleration  measure¬ 
ments,  a  detailed  study  of  the  aerodynamic  coefficients  of  the  Agena  was 
performed  (Ref.  3).  The  coefficients  were  obtained,  in  the  free  molecule 
flow  regime,  from  the  analysis  of  Sentman  (Ref.  4),  using  molecular  speed 
ratios  determined  from  the  Jacchia  1965  model  atmosphere  (Ref.  5)  modified 
by  Bruce  (Ref.  6). 
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Perigee,  and  Latitude  of  Perigee  vs  Time 


For  purposes  of  computing  the  aerodynamic  coefficients,  the  Agena 
configuration  geometry  shown  in  Fig.  2  was  used.  The  frontal  section  was 
assumed  to  be  a  flat  plate,  whereas  in  actuality  it  is  composed  of  a  number 
of  instrument  components  having  shapes  like  spheres,  boxes,  an  open  section, 
and  some  truss  members.  The  assumption  that  this  surface  is  flat  instead  of 
being  composed  of  complex  concavities  is  not  considered  to  significantly 
affect  the  estimate  of  the  overall  drag  coefficient  (Ref.  7).  The  aft  portion 
of  the  Agena,  in  the  region  of  the  structural  elements  near  the  rocket  nozzle, 
was  assumed  to  be  a  closed  cylinder  for  aerodynamic  purposes.  The  surface 
of  the  Agena  is  considered  to  be  rough  on  the  molecular  scale;  therefore, 
diffuse  molecular  reflection  was  assumed. 

In  Fig.  3  the  drag  coefficient  at  zero  angle  of  attack  is  shown  for  a 

representative  value  of  exospheric  temperature  to  illustrate  its  variation 

with  altitude.  Free  molecule  flow  is  prevalent  above  75  nmi.  Similar 

curves  exist  for  other  values  of  exospheric  temperature  and  angle  of  attack. 

During  the  LOGACS  flight  the  vehicle  angle  of  attack  ranged  between  about 

±2°.  This  resulted  in  no  appreciable  effect  on  drag  coefficient  variation, 

although  the  variation  in  exospheric  temperature  and  a  titude  which  occurred 

during  the  flight  did  affect  the  drag  coefficient  by  a  measurable  amount. 

However,  to  avoid  any  possible  confusion  in  the  interpretation  of  the  density 

data  contained  in  this  report,  a  constant  value  of  the  drag  coefficient  (equal 

to  2.6)  has  been  purposely  adopted  throughout.  This  corresponds  to  a 

2 

ballistic  coefficient  C^A/W  equal  to  0.0265  ft  / lb.  For  the  density  data 
contained  in  this  report,  the  effect  of  not  accounting  for  a  variable  is 
considered  to  fall  within  the  uncertainty  associated  with  the  itself,  which, 
under  the  assumption  of  diffuse  reflection,  amounts  to  about  ±10  percent. 

C.  SOLAR  ACTIVITY 

During  the  latter  part  of  the  LOGACS  flight  the  atmosphere  was 
strongly  perturbed  by  extreme  solar  activity.  The  accelerometer  measure¬ 
ments  that  paralleled  the  several  solar  phenomena  have  revealed  enormous 
variations  in  atmosphere  density  never  before  observed. 
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Fig.  3.  Drag  Coefficient  at  Zero  Angle  of  Attack  vs  Altitude  for  Agena 


Two  commonly  used  indices  of  solar  activity  are  shown  in  Fig.  4 

for  the  period  of  the  LOGACS  experiment.  Geomagnetic  activity,  as  shown 

by  the  planetary  range  index  K  ,  was  fairly  quiet  for  the  first  40  revs,  of 

P 

the  satellite.  At  about  Rev.  42  a  great  geomagnetic  storm  commenced,  and 

near  the  end  of  25  May  the  K  index  reached  its  maximum  value  of  9  for  two 

P 

successive  3-hr  reporting  periods  --  a  similar  event  having  occurred  only 

once  before  in  the  last  30  years.  The  great  storm,  which  was  accompanied 

by  auroral  activity  as  far  south  as  New  Mexico,  resulted  from  two  class  2B 

flares  and  one  class  3B  flare  on  23  May,  during  which  time  solar  X  rays  were 

also  detected.  The  10.  7-cm  flux,  usually  associated  with  solar  EUV, 

2  -22 

increased  rapidly  during  this  period  from  170  (W/m  /Hz  X  10  )  to  220. 

These  levels  were  also  well  above  the  average  value  of  144  which  prevailed 
over  the  previous  90  days. 

Solar  protons  with  energies  of  0.5  to  2  MeV  were  detected  frequently 
during  this  period,  with  intensities  two  to  six  times  the  normal  background 
level.  On  25  May  at  1335  GMT,  an  extremely  intense  proton  shower  nearly 
20  times  higher  than  the  normal  level  was  recorded.  At  1600  GMT,  the 
protons  began  to  fall  off  and  within  0.  5  hr  had  reduced  by  a  factor  of  5  from 
the  peak  level. 

Various  types  of  ionospheric  disturbances  were  also  prevalent  during 
the  entire  period  of  the  LOGACS  flight. 

D.  ATMOSPHERE  DENSITY  DATA 

The  atmosphere  density  data  contained  herein  have  been  obtained 
from  hand-reduced  drag  acceleration  measurements  from  several  selected 
orbital  revs,  which  span  the  LOGACS  flight.  The  drag  measurements  were 
obtained  directly  from  the  Bell  MESA  accelerometer  data  during  the  fixed 
modes  of  operation.  In  the  rotating  modes,  the  drag  data  were  obtained  by 
linear  interpolation  with  time  of  the  acceleration  measured  at  each  instant 
the  accelerometer  was  pointing  along  the  drag  vector. 

Appendix  A  contains  plots  of  the  drag  accelerations,  while  Appendix B 
contains  plots  of  computer  reduced  densities.  Section  II  of  this  volume  deals 
with  some  analyses  of  these  computer-reduced  densities. 
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Fig.  4.  Solar  Indices  Observed  During  LOGACS  Experiment 
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The  drag  acceleration  values  are  considered  to  be  accurate  to  less 
than  1  percent  based  upon  the  stable  performance  of  the  Bell  MESA  accelerom¬ 
eter  which  was  demonstrated  throughout  the  experiment.  The  density  values 
deduced  from  the  acceleration  measurements  are  considered  accurate  to 
wuhin  ±10  percent,  due  primarily  to  the  uncertainty  associated  with  the  drag 
coefficient. 

Variation  of  the  density  data  with  altitude,  latitude,  and  time  is 
shown,  and  a  number  of  interesting  comparisons  are  presented. 

1.  ATMOSPHERE  DENSITY  BETWEEN  75  AND  120  NMI  BEFORE 
GREAT  MAGNETIC  STORM 

Atmosphere  density  values  obtained  from  satellite  Revs.  26,  27,  28, 
and  31  have  been  selected  to  illustrate  conditions  prior  to  the  great  magnetic 
storm.  These  revolutions  occurred  during  24  May  1967  (see  Fig.  4),  when 
the  Kp  level  was  less  than  4.  These  revs,  are  also  particularly  interesting 
because  they  pass  through  the  high  latitude  atmosphere  density  bulge  found 
and  later  modeled  by  Jacobs  (Refs.  8  and  9).  The  size  and  location  of  this 
density  bulge  is  shown  in  Fig.  5  together  with  the  gro  ind  traces  of  the 
LOGACS  satellite  for  Revs.  26,  27,  28,  31.  Note  that  the  points  of  mini¬ 
mum  altitude  of  the  orbit  are  near  the  vicinity  of  the  bulge  and  that  Rev.  28 
passes  directly  through  the  center,  and  most  dense  part,  of  the  distribution. 

The  LOGACS  densities  versus  altitude  obtained  during  these  revs, 
are  shown  in  Figs.  6,  7,  and  8.  As  indicated,  the  solid  circles  represent 
data  obtained  as  the  satellite  was  approaching  perigee  and  the  open  circles 
the  data  leaving  perigee.  For  the  most  part,  the  data  are  shown  at  1-min 
intervals,  although  some  gaps  exist.  During  the  "toward-perigee"  por¬ 
tion  of  the  trace,  the  satellite  was  moving  through  the  high  latitudes  of 
the  northern  hemisphere,  and  going  south.  For  the  "away-from-perigee" 
portion,  the  satellite  was  in  the  midlatitudes  and  approaching  the  equator. 

This  is  the  situation  for  all  the  data  presented  in  this  report.  Therefore,  for 
every  case,  the  toward-perigee  data  represent  high  latitude  measurements 
and  the  away-from-perigee  data  represent  mid  and  equatorial  latitudes. 
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Upon  close  comparison  of  the  data  from  Revs.  26,  27,  and  28, 
little  evidence  can  be  found  to  support  the  existence  of  the  high  latitude 
density  bulge  modeled  by  Jacobs.  If  it  did  exist  at  this  location,  the  densities 
obtained  toward  perigee  in  Rev.  28  would  be  significantly  higher  ( — 30  percent) 
than  those  obtained  in  Rev.  27,  which  in  turn  would  be  higher  than  the 
density  toward  perigee  in  Rev.  26.  This  is  not  the  case. 

It  cannot  be  concluded,  however,  that  such  a  bulge  does  not  exist  at 
some  other  location  on  the  earth  or  under  different  conditions  of  solar 
activity.  Nor  can  it  be  concluded,  as  did  Jacchia  and  Slowey  (Ref.  10),  that 
the  bulge -attributed  phenomena  are  entirely  of  nonatmospheric  origin,  as 
will  be  evident  shortly.  In  fact,  in  the  light  of  unpublished  density  data 
recently  obtained  from  satellite  tracking  data,  Jacobs  now  contends  that 
there  is  clearly  a  bulge  due  to  atmospheric  causes  but  that  it  sporadically 
changes  position  from  time  to  time  or  even  disappears  completely  (Ref.  11). 

Figure  9  shows  the  LOGACS  density  data  for  Rev.  31,  which  passes 
over  the  north  magnetic  dip  pole.  The  Rev.  31  data,  contrary  to  the  earlier 
data,  do  suggest  increased  density  at  the  high  latitudes.  This  is  best 
illustrated  by  comparison  of  Figs.  10  to  13,  which  show  normalized  density 
for  these  revs,  versus  latitude.  The  normalization  eliminates  the  dominant 
altitude  variation  and  was  accomplished  by  dividing  the  LOGACS  density 
values  by  the  reference  density  model  shown  in  Fig.  14  (Refs.  5  and  6). 
Note,  in  particular,  that  the  density  between  50°N  and  70°N  during  Rev.  31 
is  20  percent  greater  than  the  earlier  revs.,  which  clearly  indicates  the 
possibility  of  a  bulge  over  the  region  of  the  magnetic  dip  pole  even  during 
conditions  of  quiet  geomagnetic  activity. 

Further  analysis  of  the  data  in  the  region  of  the  dip  pole  may  provide 
better  definition  of  this  effect. 

It  is  also  evident  in  Figs.  10,  11,  and  13  that  there  is  a  systematic 
increase  in  density  in  the  region  of  the  sub-solar  point  at  about  10°  N  to 
20°  N  latitude.  This  is  very  likely  a  diurnal  bulge  effect,  in  part,  since  the 
orbital  plane  was  about  26°  west  of  the  noon  meridian. 
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Fig.  9.  Density  vs  Altitude  for  Rev.  31 
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Fig.  10.  Normalized  Density  vs  Latitude  for  Rev.  26 
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Normalized  Density  vs  Latitude  for  Rev.  27 
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Fig.  12.  Normalized  Density  vs  Latitude  for  Rev.  28 
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Fig.  13.  Normalized  Density  vs  Latitude  for  Rev.  31 
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2. 


ATMOSPHERE  DENSITY  BETWEEN  75  AND  120  NMI 
DURING  GREAT  MAGNETIC  STORM 


To  illustrate  conditions  in  the  atmosphere  during  the  great  magnetic 
storm,  dersity  values  obtained  from  satellite  Revs.  46,  47,  48,  and  57  were 
selected.  The  ground  traces  for  these  revs,  are  shown  in  Fig.  15. 

The  density  data  exhibit  enormous  variations  in  response  to  the 
extreme  solar  and  geomagnetic  activity  that  was  prevalent  on  25  and  26  May 
(see  Fig.  4)  as  shown  in  Figs.  16  to  21.  In  general,  the  toward-perigee  (high 
latitude)  data  are  more  oscillatory  than  the  away-from-perigee  (low  latitude) 
data,  although  regular  patterns  exist.  Density  variations  of  50  percent  at 
the  same  altitude,  but  different  latitudes,  which  occur  within  the  same  rev.  , 
are  not  uncommon,  as  shown  in  Fig.  16.  In  Rev.  47,  Fig.  17,  a  50  percent 
spatial  change  in  density  is  seen  to  occur  at  81  nmi  (150  km).  A  regular  pat¬ 
tern  of  behavior,  not  evident  in  Fig.  17,  is  clearly  seen  when  the  same  data 
are  separated  as  shown  in  Figs.  18  and  19. 

Rev.  48  data,  shown  in  Fig.  20,  indicate  the  same  tendency  for 
density  to  begin  decreasing  with  decreasing  height  at  about  85  nmi  as  was 
evident  in  Rev.  47. 

Rev.  57  data,  shown  in  Fig.  21,  have  been  included  to  show  a  fac¬ 
tor  of  2  increase  in  density  which  was  measured  near  79  nmi  at  two  points  over 
the  earth  separated  by  only  4°  along  the  orbital  path. 

The  density  values  for  these  revs. ,  normalized  and  plotted  versus 
latitude  in  Figs.  22  to  25,  clearly  reveal  the  large  fluctuating  spatial  and 
temporal  behavior  that  was  characteristic  during  the  storm. 

In  Fig,  22,  during  Rev,  46,  an  immense  two -pronged  high  latitude 
bulge  almost  70  percent  greater  than  the  reference  density  was  measured. 

The  diurnal  bulge,  although  small  in  comparison,  was  still  evident  between 
10°  N  and  20°  N.  By  Rev.  47,  1. 5  hr  later,  shown  in  Fig.  23,  the  high 
latitude  bulge,  still  two -pronged,  had  diminished  somewhat,  but  the  "diurnal" 
bulge  increased  (the  quotes  signifying  that  this  effect  may  have  been  influenced 
by  other  phenomena,  e.g. ,  the  equatorial  electrojet,  Refs.  12  and  13). 

Finally,  by  Rev.  48,  Fig.  24,  the  high  latitude  bulge  had  essentially 


1-21 
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Fig.  15.  LOGACS  Satellite  Ground  Traces  for 
Revs.  46,  47,  48,  and  57 
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Fig.  22.  Normalized  Density  vs  Latitude  for  Rev.  46 
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Fig.  23.  Normalized  Density  vs  Latitude  for  Rev.  47 


disappeared,  but  the  diurnal  bulge  had  doubled  in  amplitude.  Then,  by 
Rev.  57,  half  a  day  later,  shown  in  Fig.  25,  the  high  latitude  bulge  had 
reappeared,  and  the  diurnal  bulge  was  even  more  intense  than  before. 

These  results  suggest  the  existence  of  a  high  latitude  density  bulge 
that  is  dependent  upon  the  level  of  geomagnetic  activity,  and  contradict  the 
high  latitude  bulge  model  of  Jacobs  (Ref.  9),  which  is  not  dependent  upon 
magnetic  activity.  Also,  the  pronounced  diurnal  bulge  that  is  evident  con¬ 
flicts  with  most  atmosphere  models  (Refs.  5,  9,  14  and  15),  which  often  show 
decreases  in  density  in  the  region  of  the  sub-solar  point  at  low  altitudes. 

3.  VARIATION  IN  DENSITY  AT  90  NMI  DURING  INTENSE 
25  MAY  PROTON  SHOWER 

As  mentioned  previously,  an  intense  proton  shower  nearly  20  times 
greater  than  normal  levels  was  recorded  on  25  May  1967,  reaching  its  peak 
at  1335  GMT.  At  1600  GMT,  the  proton  intensity  began  to  decrease,  and 
within  0.5  hr  it  had  been  reduced  by  a  factor  of  5  from  its  peak  level. 

The  LOGACS  data  corresponding  to  that  time  period  were  examined 
to  determine  the  effects,  if  any,  that  might  be  correlated  with  the  proton 
shower.  A  marked  density  response  was  found,  at  about  90  nmi  and  80°  N 
latitude,  that  may  be  the  result  of  highly  localized  particle  precipitation 
related  to  the  proton  event  (Ref.  16).  This  is  illustrated  in  Fig.  26  where 
the  density  values  for  this  time  period  are  shown.  As  time  progressed,  it 
can  be  seen  that  the  density  at  about  90  nmi  increased  by  35  percent  shortly 
after  the  recorded  peak  of  the  proton  shower.  Then,  by  1630  GMT,  the 
density  returned  to  its  prestorm  levels  with  a  rapidity  totally  unexpected 
according  to  present  one-dimensional  theory  (Ref,  17). 

4.  RELATIVE  DENSITY  VS  TIME  FOR  SPECIFIC  ALTITUDES 

With  the  LOGACS  data,  it  has  been  possible  for  the  first  time  to 
observe  atmosphere  density  as  a  function  of  time  at  specific  altitudes. 
Examples  are  shown  here  for  100,  90,  and  80  nmi,  and  the  minimum  altitude 
experienced.  During  each  rev.,  each  of  the  selected  altitudes  (except 
the  minimum)  was  reached  twice  by  the  satellite.  Therefore,  the  density 
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Fig.  25.  Normalized  Density  vs  Latitude  for  Rev.  57 
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values  presented  have  been  separated  according  to  whether  they  were 
measured  going  toward  perigee  (actually  minimum  altitude)  or  measured 
going  away  from  perigee.  To  aid  in  the  interpretation,  the  density  data 
were  normalized  in  the  same  manner  as  was  done  previously. 

In  Fig.  27,  normalized  density  toward  perigee  is  shown  versus  time 

for  the  altitudes  indicated.  The  geomagnetic  planetary  amplitudes  a  is  also 

P 

shown  for  comparison.  The  maximum  relative  density,  shown  at  the  bottom 
of  the  figure,  corresponds  to  the  minimum  altitudes  reached  by  the  satellite, 
which  varied,  for  the  time  period  shown,  between  78. 1  and  75,6  nmi. 

Several  extremely  striking  features  are  immediately  evident, 
namely,  the  remarkable  correlation  with  magnetic  activity  and  the  varying 
response  with  altitude,  both  in  amplitude  and  in  phase.  Also,  there  are  some 
instances  when  the  densities  are  increasing  with  time  at  one  altitude  and 
decreasing  at  another.  This  is  clear,  for  example,  on  24  May  at  90  and 
100  nmi  and  suggests  that  the  energy  responsible  for  this  response  may  be 
deposited  between  these  two  altitudes. 

In  Fig.  28,  the  corresponding  normalized  density  away  from  perigee 
is  shown.  Here  it  is  evident  that  the  variations,  although  complex,  are  much 
less  so  than  those  seen  toward  perigee.  However,  it  is  interesting  to  note  that 
the  increase  in  density  due  to  the  great  magnetic  storm  is  greater  in  this 
case  than  it  was  in  the  toward-perigee  portion  of  the  orbit. 

5.  ATMOSPHERE  DENSITY  AMPLITUDE  AND  PHASE 
RESPONSE  DUE  TO  GREAT  MAGNETIC  STORM 

When  the  relative  density  data  shown  in  Figs.  27  and  28  are  examined 
in  detail,  it  can  be  seen  that  the  amplitude  and  phase  of  the  response  to  the 
great  magnetic  storm  vary  with  altitude.  The  amplitude  and  phase  are  also 
dependent  on  latitude  in  that  the  response  is  markedly  different  toward  perigee 
(high  latitude)  than  it  is  away  from  perigee  (low  latitude). 

The  results  of  a  preliminary  analysis  of  the  phase  response  are 
shown  in  Fig.  29.  The  phase  (or  time  lag)  is  shown  versus  altitude  for  the 
toward-perigee  and  away-from-perigee  cases,  and  was  determined  by  com¬ 
paring  the  onset  of  the  storm  with  the  onset  of  the  response.  There  is 
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fig.  27.  Normalized  Density  vs  Time  (Toward  Perigee) 


1-36 


Fig.  28.  Normalized  Density  vs  Time  (Away  from  Perigee) 
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considerable  variation  with  altitude  and  latitude,  ranging  between  almost  no 
time  lag  to  a  lag  of  greater  than  6  hr.  Therefore,  these  results  support  in 
a  general  sense  the  variable  time  lags  found  by  DeVries,  Friday,  and  Jones 
(Ref.  18),  but  conflict  somewhat  with  the  results  of  Jacchia,  Slowey,  and 
Verniani  (Ref.  19),  who  found  only  slight  variation  in  time  lag  with  latitude 
from  the  drag  data  on  four  high  inclination  satellites. 

The  amplitude  of  the  density  response  to  the  magnetic  storm,  shown 
in  Fig.  30,  has  considerable  variation  with  altitude  and  latitude  also.  Here, 
it  is  easy  to  see,  as  mentioned  previously,  that  the  density  increase  at  low 
latitudes  (away  from  perigee)  is  greater  than  at  high  latitudes  (toward 
perigee) --a  most  interesting  result. 

6.  COMPARISON  OF  LOGACS  MEASURED  DENSITY  WITH 
DENSITY  OBTAINED  FROM  ORBITAL  DECAY  OF 
LOGACS  SPACECRAFT,  SATELLITE  1967  50B 

The  tracking  data  from  the  LOGACS  spacecraft,  Satellite  1967  50B, 
have  been  reduced  by  Schusterman  (Ref.  2)  to  obtain  perigee  density  as  a 
function  of  time  based  upon  observed  orbital  decay  rates.  We  have  been  pro¬ 
vided,  therefore,  with  a  unique  opportunity  to  compare  the  density  values 
determined  in  this  manner  with  the  actual  density  measurements  obtained  by 

a  sensor  on  the  same  vehicle.  This  comparison  is  shown  in  Fig.  31.  In 

2 

Ref.  2  a  value  of  the  ballistic  coefficient  C^A/W  equal  to  0.0311  ft  /lb  was 
assumed.  We  believe  that  this  value  was  based  upon  an  erroneous  estimate 
of  the  spacecraft  weight  and  have  therefore  adjusted  the  density  values 
derived  by  Schusterman  to  correspond  to  the  ballistic  coefficient  used 
throughout  this  report. 

The  trend  of  the  orbital  decay  density,  as  can  be  seen  in  Fig.  31, 
parallels  the  LOGACS  data,  but  there  arc  discrepancies  considered  to  be 
extremely  significant.  Note,  in  particular,  that  the  sharp  peak  in  the  orbital 
decay  density  which  occurs  during  25  May  lags  the  corresponding  peak  in  the 
measured  density  by  0.  15  day  (3.6  hr).  It  would  appear,  in  fact,  that  much 
closer  agreement  would  be  accomplished  if  the  entire  orbital  decay  density 
curve  was  shifted  in  time  by  the  above  amount.  Since  precise  knowledge  of 
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Fig,  31,  Density  Comparison:  Orbital  Decay  vs  LOGACS 
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the  phasing  of  short  term  phenomena  in  the  atmosphere  is  required  to 
successfully  construct  short  term  density  models,  it  is  clear  that  great  care 
is  required  when  determining  density  values  from  orbital  decay. 

7.  COMPARISON  OF  LOGACS  MEASURED  DENSITY  WITH 
PREVALENT  MODEL  ATMOSPHERES 

The  LOGACS  data  obtained  throughout  the  flight  were  scanned  to 
obtain  representative  density  values  that  could  be  compared  with  currently 
used  model  atmospheres.  This  was  done  for  two  altitudes  and  two  levels  of 
magnetic  activity.  The  density  values  selected  as  representative  of  these 
conditions  are  shown  in  Table  2. 

Values  of  density  that  would  be  predicted  by  several  prevalent  model 
atmospheres  (Refs,  5,  6,  9,  14,  and  15)  for  these  altitudes  and  levels  of  mag¬ 
netic  activity  were  then  compared  with  the  selected  LOGACS  densities.  This 
comparison  is  shown  in  Table  3.  As  can  be  seen,  there  is  an  enormous 
spread  in  the  results,  indicating  how  unsuccessful  recent  modeling  attempts 
have  been.  In  general,  the  best  agreement  is  at  the  higher  altitude,  108  nmi, 
which  is  probably  due  to  the  fact  that  density  data  at  this  altitude  have  been 
more  abundant.  The  agreement  at  tile  lower  altitude,  81  nmi,  is  extremely 
poor- -the  most  serious  contradiction  occurring  with  the  U.S.  19''6  summer 
model.  Curiously,  much  better  agreement  would  have  been  obtained,  in  this 
case,  if  the  U.S.  1966  winter  model  had  been  used  instead.  The  LOGACS 
data  seem,  therefore,  to  be  in  considerable  conflict  with  the  seasonal  results 
found  by  Champion  (Ref.  20),  which  were  extended  to  form  the  bases  of  the 
U.S.  1966  seasonal  models  (Ref.  21). 
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Table  2.  Representative  LOGACS  Density  Measurements 


Quiet  Kp 

Active  Kp 

(Kp 

=  0-4) 

7-8) 

8 1  nmi 

108  nmi 

81  nmi 

108  nmi 

(150  km) 

(200  km) 

(150  km) 

(200  km) 

LOGACSp , 
gm/cin^ 

2.8  X  10‘12 

- 1 3 

3.  9  X  10  1 

- 12 

3.4  X  10 

a 

Not  shown  because  of  extreme  temporal  and  spatial  variability. 


Table  3.  Comparison  of  LOGACS  Data  with  Density  Models 


P/PM  •'< 

percent 

Quiet  Kp 

Tffl  =  1500°  K 

Active  Kp 

T  a,  =  1800°K 

81  nmi 
(150  km) 

108  nmi 
(200  km) 

81  nmi 
(150  km) 

108  nmi 
(200  km) 

LMSC  67a 

- 1 3  percent^ 

-12b 

-  2C 

d 

Jacchia  65  (Bruce 
mod. ) 

+13  percent 

-11 

+  31 

d 

Jacchia  65 

+31  percent 

-  3 

+60 

d 

Cira  65 

+  34  percent0 

+26e 

+60f 

d 

U.S.  66  summer 

+51  percent 

+  1 

+84 

d 

aVersion  used  here  is  an  unpublished  modification  of  Ref.  9  submitted  by 
Lockheed  Missiles  and  Space  Company  (LMSC), 

bS  =  200,  diurnal  maximum. 

CS  =  250,  diurnal  maximum. 

aNot  shown  because  of  extreme  temporal  and  spatial  variability. 
eModel  5,  Hour  14. 

^  Model  8,  Hour  14. 
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SECTION  II 
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ANALYSIS  OF  UPPER  ATMOSPHERE  DENSITY 

A.  COMPARISONS  OF  LOGACS  DENSITIES  WITH  MODELS 

One  of  the  most  striking  features  of  the  LOGACS  density  data  is  the 
extreme  variability.  This  has  been  commented  upon  with  regard  to  the  hand- 
reduced  data  reported  in  the  previous  section.  These  hand-reduced  cases 
are  not  atypical,  as  a  glance  at  Appendix  B  will  confirm.  In  many  places 
(for  examples  see  Rev.  48,  approaching  perigee,  and  Rev.  53,  approaching 
perigee),  density  decreases  with  decreasing  altitude.  Further,  in  many 
places  pronounced  density  differences  can  be  seen  at  the  same  altitude  but 
different  latitudes  by  comparing  the  two  plots  for  each  rev. 

A  quick  comparison  of  the  LOGACS  densities  with  model  atmospheres 
deserves  mention.  Extrapolation  of  the  plots  in  Appendix  B  to  120  km  will 
confirm  that  the  invariant  boundary  conditions  assumed  at  this  altitude  for 
most  model  atmospheres  are  incorrect.  As  a  result,  these  plots  differ 
markedly  from  plots  of  model  atmospheres  derived  from  similar  conditions, 
and  show  one  reason  why  model  atmospheres  have  great  difficulty  in  "fitting" 
the  LOGACS  data. 

Comparison  of  the  LOGACS  densities  with  models  having  the  same 
values  of  exospheric  temperature  reveals  pronounced  differences  in  slope 
and  in  other  characteristics.  The  LOGACS  data  for  Rev.  5,  leaving  perigee, 
are  displayed  as  Curve  1  in  Fig.  32.  This  rev.  was  characterized  by  a 
computed  exospheric  temperature  T^  of  1164°K.  Curve  2  in  Fig.  32  is  the 
Jacchia  static  diffusion  model  for  Tm  =  1150’K;  while  representing  the  density 
very  well  near  perigee,  it  is  inadequate  beyond  170  km.  However,  a  Jacchia 
static  diffusion  model  with  T  10()0*K  (Curve  3)  represents  the  data  much 
better.  These  Jacchia  models  are  obtained  from  Ref.  5. 
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The  problem  is  not  one  of  a  simple  difference,  however.  Curve  1 
of  Fig.  33  shows  the  LOGACS  data  for  Rev.  59,  leaving  perigee.  A  much 
higher  value  of  T^,  =  1408°  K  is  associated  with  this  revolution.  Curve  2  in 
this  figure  shows  the  data  from  Jacchia's  static  diffusion  model  for 
Tgj  =  1400°K.  It  is  readily  apparent  that  this  curve  in  no  way  represents 
the  data;  it  is  consistently  lower  than  the  data.  A  drastic  increase  in  the 
model  temperature,  to  Tm  =  2100°K  (Curve  3),  does  not  improve  the  agree¬ 
ment  to  a  significant  degree.  These  two  examples  serve  to  show  clearly  that 
the  Jacchia  model  dependence  on  temperature  is  not  correct. 

The  behavior  of  model  atmospheres  in  predicting  LOGACS  densities 
is  consistently  poor  for  all  models.  The  Jacchia  model,  used  above  for  com¬ 
parison,  was  chosen  because  it  is  considered  to  be  representative  of  the 
state  of  the  art.  In  fact,  comparisons  shows  that  the  Jacchia  model  is,  if 
anything,  slightly  better  than  the  others  considered. 

Ratios  of  LOGACS  density  to  model  atmosphere  density  can  be  com¬ 
puted  by  the  Computer  Printout  of  Densities  (CPD)  program  mentioned  earlier. 
Figures  34  and  35  present  portions  of  a  plot  of  such  ratios  against  time,  using 
the  Walker -Bruce  modification  of  Jacchia's  model  atmosphere  (Ref.  6).  Also 
plotted  against  time  is  vehicle  altitude,  for  visual  reference.  These  graphs 
show  a  clear  correlation  of  the  density  ratio  with  altitude.  They  further 
show  that  this  model  does  a  better  job,  relatively  speaking,  of  predicting 
densities  near  perigee  than  near  apogee.  In  fact,  the  ratio  becomes  very 
poor  near  apogee.  In  some  cases  near  apogee,  large  "spikes"  occur  (see 
1870  min  on  Fig,  34,  for  example);  these  are  as  yet  unexplained.  The  larger 
ones  occur  at  intervals  with  periods  that  are  near  integral  multiples  of  6  hr; 
smaller  peaks  of  this  sort  occur  on  nearly  every  orbit  (period  1.  5  hr).  These 
spikes,  which  occur  near  apogee  where  the  data  quality  was  poor,  may  be  a 
result  of  the  data  reduction  process,  and  may  not  represent  actual  behavior. 

In  general,  all  data  obtained  near  apogee  are  suspect;  however,  the  data 
obtained  below  250  km  are  of  excellent  quality.  Figure  35  shows  the  portion 
of  the  ratio  data  in  which  the  effects  of  the  storm  were  most  pronounced;  the 
failure  of  the  model  to  adequately  respond  is  very  evident  in  this  region. 
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®  REV  59,  LEAVING  PERIGEE ,  T®  =  1408 0  K 

(2)  JACCHIA  STATIC  DIFFUSION  MODEL, T®  =  I400°K 

®  JACCHIA  STATIC  DIFFUSION  MODEL, T®  =  2IOO#K 


ALTITUDE  h,  km 


Fig.  33.  Comparison  of  LOGACS  Data  (Rev.  59,  Leaving  Perigee) 
with  Jacchia  Static  Diffusion  Model 
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Figure  35.  Ratio  of  LOGACS-Derived  Density  to  Walker-Bruce  Model  Density  for 
a  Portion  of  25  May  1967;  Vehicle  Altitude  Included  for  Reference 


For  comparison,  Figs.  36  and  37  present  the  ratio  of  the  LOGACS 
density  to  the  1962  U.S.  Standard  Atmosphere  density  (Ref.  22).  This  atmo¬ 
sphere,  being  a  table  lookup  model,  incorporates  no  response  to  varying 
solar  and  geomagnetic  parameters.  Its  behavior  is  much  worse  than  the 
Jacchia  model  in  predicting  LOGACS  densities.  Figure  37,  which  contains 
the  time  interval  of  the  storm's  effects,  shows  the  complete  failure  of  this 
model  under  these  conditions. 

The  two  models  discussed  above  are  the  best  and  worst  cases, 
respectively,  of  all  models  considered  in  attempting  to  predict  LOGACS 
densities.  Clearly  the  comparisons  show  that  the  models  are  fundamentally 
in  error  in  the  lower  altitude  region  (140  to  200  km)  in  particular  and  that  the 
behavior  of  the  models  under  conditions  characterized  by  high  solar  and 
geomagnetic  indices  at  all  altitudes  is  inadequate. 

The  ratios  of  LOGACS  density  to  the  density  of  the  Jacchia  -  Walker  - 
Bruce  model  were  examined  then  at  fixed  altitudes.  Figures  38  and  39  show 
the  approaching-perigee  and  away  from  perigee  data,  respectively. 

It  is  obvious  that,  especially  during  the  storm,  the  model  represents 

the  data  rather  poorly.  The  minimum  altitude  case  shows  a  single  large  peak 

on  26  May  which  lags  the  peak  in  the  planetary  range  index  K  (see  Fig.  28) 

P 

by  some  10  hr. 

The  data  at  all  altitudes  both  toward  and  away  from  perigee  show  a 
number  of  large  peaks  during  25  and  26  May.  The  relative  size  and  separa¬ 
tion  of  these  peaks  vary  with  altitude  and  latitude  in  a  manner  as  yet  not 
understood. 

B.  INTERNAL  COMPARISONS  OF  LOGACS  DENSITIES 
1.  SORTING  BY  TIME 

The  LOGACS  data  is  broken  into  seven  groupings,  called  vectors,  as 
a  result  of  the  postflight  ephemeris  computation  (Ref.  1).  These  vectors  are 
of  differing  lengths,  but  serve  as  convenient,  albeit  arbitrary,  data  spans  for 
analysis  purposes.  Figures  40  through  46  present  linear  plots  of  density 
versus  altitude  for  each  of  the  seven  vectors;  these  plots  have  been  annotated 
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Fig.  36.  Ratio  of  LOGACS- Derived  Density  to  1962  U.  S.  Standard 
Atmosphere  Density  for  a  Portion  of  23  May  1967; 
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Fig.  37.  Ratio  of  LOGACS-Derived  Density  to  1962  U.  S.  Standard 
Atmosphere  Density  for  a  Portion  of  25  May  1967; 
Vehicle  Altitude  Included  for  Reference 


TOWARD  PERIGEE 


I  23  MAY  I  24  MAY  I  25  MAY  I  26  MAY 

Fig.  38.  The  Density  Ratio  LOGACS/ Jacchia- Walker-Bruce 
Model  for  High  Latitude  (Toward  Perigee)  Data 
at  Various  Altitudes 
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Fig.  39.  The  Density  Ratio  LOGACS/ Jacchia- Walker-Bruce 
Model  for  Low  Latitude  (Away  from  Perigee)  Data 
at  Various  Altitudes 
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VECTOR  I 
(REVS.  5-14) 

23  MAY  1967:  1/45/9 -16/1/43 
MAX  flp  3.22  MAX  F|0  195.2 
MIN  op  5.65  MIN  F,0  189.3 
Fl0  144.0 
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Fig.  40.  Density  vs  Altitude  for  Vector  1 
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Fig.  41.  Density  vs  Altitude  for  Vector  2 
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Fig.  42.  Density  vs  Altitude  for  Vector  3 
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VECTOR  4 
(REVS.  33-44) 
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Fig.  43.  Density  vs  Altitude  for  Vector  4 
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Fig.  44,  Density  vs  Altitude  for  Vector  5 
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Fig.  46.  Density  vs  Altitude  for  Vector  7 
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to  reflect  the  time  span,  and  the  solar  and  geomagnetic  conditions  which 
existed  during  the  time  span,  of  each  vector. 

In  Fig.  40  (Vector  1),  there  is  a  spread  in  density  values  of  33  per¬ 
cent  at  150  km.  This  vector  is  characterized  by  relatively  constant  solar 
and  geomagnetic  activity  and  is  well  before  the  disturbances.  A  part  of  this 
spread  is  attributable  to  latitudinal  effects,  but  a  good  portion  of  it  is  due  to 
local  temporal  perturbations  of  the  density.  The  effects  of  the  storm  are 
most  clearly  seen  in  Figs.  43  and  44  (Vectors  5  and  6).  Corresponding  to 
the  above,  the  spread  in  density  values  at  150  km  is  55  percent  in  Fig.  43. 

2.  SORTING  BY  ALTITUDE 

The  interpolative  features  of  the  CPD  program  we  re  utilized  to 
determine  density  and  other  parameters  at  reference  altitudes  throughout  the 
flight.  To  make  the  variation  discussed  in  the  previous  section  more  quan¬ 
titative,  some  simple  statistics  were  performed  on  the  densities  at  the  refer¬ 
ence  altitudes.  A  mean  and  rms  of  the  density  for  each  vector  were  com¬ 
puted  at  each  altitude.  These  results  are  graphically  displayed  in  Figs. 

47  through  54.  The  mean  was  plotted  at  the  time  midpoint  of  its  vector,  and 
error  bars  are  used  to  denote  a  ±lff  variation.  Vector  2  was  extremely  short 
and  therefore  is  not  a  valid  sample. 

Readily  apparent  from  these  plots  is  the  increase  in  the  spread  of 
density  values  caused  by  the  storm,  as  well  as  the  increase  in  the  average 
density  caused  by  the  storm.  Apparently  the  storm  caused  different 
phenomena  at  various  levels  in  the  atmosphere.  For  example,  at  148  km 
and  150  km  the  average  density  was  still  increasing  with  time  in  Vector  7 
(although  the  la-  spread  was  somewhat  reduced).  This  can  be  compared  with 
the  data  beyond  185  km,  where  an  apparent  return  to  unperturbed  conditions 
commenced. 
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Fig.  47.  Mean  and  RMS  of  Dtnsi  mes  at  148  km  (80  nmi) 

During  LOGACS  Flight 
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Fig.  48.  Mean  and  RMS  of  Density  Values  at  150  km  (81  nmi) 
During  LOGACS  Flight 
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Fig.  52.  Mean  and  RMS  of  Density  Values  at  203  km  (110  nmi) 
During  L.OGACS  Flight 
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Fig.  53.  Mean  and  RMS  of  Density  Values  at  222  km  (120  nmil 
During  LOGACS  Flight 


SECTION  III 


ATMOSPHERIC  PREDICTABILITY 

A.  DENSITY  PREDICTION 

The  LOGACS  experiment  provided  several  days  of  near-continuous 
data  about  the  forces  the  Agena  experienced.  From  these  data,  it  was  possi¬ 
ble  to  extract  density  along  the  orbit,  based  on  theoretical  aerodynamic  con¬ 
siderations  of  the  Agena.  The  LOGACS  density  data  have  made  it  possible  to 
determine  some  idea  of  the  predictability  of  atmospheric  density  at  various 
altitudes.  Along  these  lines,  the  following  analysis  was  performed.  In  order 
to  simulate  the  average  drag  in  a  manner  analogous  to  what  is  done  in  a 
seven-parameter  orbit  determination,  LOGACS  densities  were  averaged  at 
selected  altitudes  over  suitable  rev.  spans.  Twelve -rev.  averaging  was 
chosen  for  this  study.  The  average  values  determined  in  this  way  were  con¬ 
sidered  to  correspond  to  the  average  drag  determined  from  a  fit  to  tracking 
data.  The  "predicted"  density  values  were  obtained  from  the  LOGACS  data 
for  revs,  immediately  following  the  spans  used  for  averaging.  Each  of  the 
four  revs,  adjacent  to  the  averaged  spans  was  examined  at  the  selected 
altitudes.  The  differences  between  each  of  the  density  values  in  the  "predict" 
interval  and  those  obtained  by  averaging  were  considered  to  represent  a 
measure  of  drag  predictability.  The  results  of  this  analysis  are  contained 
in  this  section. 

B.  DISCUSSION 

The  LOGACS  flight  covered  a  period  of  four  days,  totaling  some  66 
orbits.  The  data  from  49  of  these  orbits  were  sufficiently  complete  for 
inclusion  in  this  analysis.  Four  altitudes,  110,  100,  90,  and  80  nmi,  were 
chosen  as  representative  altitudes  for  study.  The  density  data  for  these 
altitudes  were  collated  by  rev.  ;  those  revs,  having  two  data  points  or  less 
were  rejected. 
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r ho  density  data  for  a  given  altitude  contained  in  all  12- rev.  spans 
were  collated  and  examined.  The  mean  density  value  and  the  rms  of  the  data 
(expressed  as  a  percentage  of  the  mean  valued  were  computed  for  each 
12 -rev.  group  at  each  altitude,  as  shown  in  Table  4. 


Fable  4.  Mean  and  RMS  Density  for  All  12-Rev.  Groups 


Altitude,  nmi 

Mean,  gm/cm3 

Rms  (Percent  of  Mean) 

110 

3.  876  X  10'13 

11.91 

100 

6.  971  X  10' 13 

11. 52 

90 

1.432  X  10'12 

10.96 

80 

3.  094  X  10~12 

11.09 

Overall 

11.  37 

Data  from  each  of  four  revs,  following  each  12 -rev.  group  were 
then  compared  with  the  mean  for  the  12-rev.  groups.  The  difference  of  this 
value  from  the  mean  was  computed  and  expressed  as  a  percentage  of  the 
mean.  These  percentage  values  were  averaged  for  each  altitude  and  for 
all  the  data;  as  a  measure  of  the  dispersion  of  the  data,  the  rms  about  the 
average  was  computed  for  each  case.  These  results  are  shown  for  the  first, 
second,  third,  and  fourth  rev.  after  the  12 -rev.  groups  in  Tables  5,  6,  7, 
and  8,  respectively. 

C.  CONCLUSION 

As  can  be  easily  seen  from  the  following  tables,  the  average  pre¬ 
dicted  dispersions  varied  between  approximately  6  and  10  percent,  with  an 
rms  of  approximately  14  percent.  The  most  obvious  feature  of  the  table  is 
the  expected;  that  the  prediction  grows  worse  as  time  increases.  Another 
expected  feature  is  the  increase  in  dispersions  with  altitude. 
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Table  5.  First  Rev.  after  12-Rev.  Mean,  Relative  Difference 


Altitude,  nmi 

Average  of  Difference, 
Percent  of  Mean 

Rms  about  Average 
of  Difference, 
Percent  of  Mean 

110 

7.  125 

±19.  534 

100 

5.  421 

±12.  738 

90 

5.  426 

±1 1.  558 

80 

5.  740 

±  8.  975 

Overall 

5.  987 

±13.  655 

Table  6.  Second  Rev.  after  12 -Rev.  Mean, 

Relative  Difference 

Altitude,  nmi 

Average  of  Difference, 
Percent  of  Mean 

Rms  about  Average 
of  Diffe  rence. 
Percent  of  Mean 

110 

8.  679 

±20. 727 

100 

6.  328 

±13.  936 

90 

6.  945 

±11.  625 

80 

6.  921 

±  9.437 

Overall 

7.  207 

_ 

±14. 516 

Table  7.  Third  Rev.  after  12 -Rev.  Mean, 

Relative  Difference 

Altitude,  nmi 

Average  of  Difference, 
Percent  of  Mean 

Rms  about  Average 
of  Difference, 
Percent  of  Mean 

110 

9.  707 

±21.  428 

100 

7.  485 

±14.  370 

90 

8.  094 

±12.  108 

80 

8.  001 

±  9.  607 

Ove  rail 

8.  303 

±14.  977 
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Table  8.  fourth  Rev.  after  12-Rev.  Mean,  Relative  Difference 


Altitude,  nmi 

Average  of  Difference, 
Percent  of  Mean 

Rms  about  Average 
of  Difference, 

Percent  of  Mean 

110 

11.  390 

±21.903 

100 

8.994 

±15. 293 

00 

9.434 

±12. 364 

80 

9.081 

±  9.503 

Overall 

9.705 

±15.415 

The  data  in  Tables  5  through  8  present  the  difference  of  the  pre¬ 
dicted  rev.  from  the  base  rev.  as  a  measure  of  prediction  accuracy.  The 
rms  about  this  difference  is  sufficiently  large  to  allow  the  predicted  value 
to  be  either  larger  or  smaller  than  the  base  case.  The  analysis  here 
implicitly  assumes  that  the  data  are  random,  when  in  fact  they  are  not.  The 
LOGACS  flight  was  characterized  by  extreme  values  of  solar  and  geomagnetic 
parameters;  see,  for  example,  Ref.  23.  These  parameters  were  increasing 
over  most  of  the  flight.  Since  density  responds  to  changes  in  these  param¬ 
eters,  all  or  part  of  the  worsening  prediction  with  time  can  be  attributed  to 
this  change. 
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SECTION  IV 


MID  LATITUDE  DENSITY  MODEL  ATMOSPHERE  BETWEEN 

120  AND  17  5  km 


A.  INTRODUCTION 

Examination  of  the  LOGACS  density  data  revealed  three  separate 
geographic  regions  having  an  apparently  very  different  behavior.  For 
purposes  of  discussion  these  regions  are  referred  to  as  "low  latitude," 
"polar,"  and  "high  altitude"  regimes.  The  low  latitude  region  lies  between 
+40  and  -30*  geocentric  latitude,  at  altitudes  up  to  about  220  km;  the  polar 
region  lies  be!  ween  +40°  latitude  and  the  North  Pole  up  to  220  km;  the  high 
altitude  region  includes  the  whole  atmosphere  above  this  altitude.  These 
limits  were  delineated  by  the  availability  of  the  density  data,  not  necessarily 
by  physical  boundaries. 

Perigee  (or,  more  accurately,  the  minimum  altitude  point)  occurred 
within  the  low  latitude  regime,  and  varied  from  about  +30  to  +43  latitude 
(Ref.  1),  with  most  of  the  data  being  obtained  when  the  minimum  altitude  point 
was  near  +35*  latitude. 

B.  DATA  REDUCTION 

In  the  low  latitude  regime,  the  density  data  were  sorted  according 
to  geocentric  latitude.  Previous  models  suggested  that  a  "bulge"  would  be 
seen  either  near  the  latitude  of  the  subsolar  point,  or  near  the  equator;  there¬ 
fore,  the  density  data  were  examined  at  5°  increments.  When  no  measure¬ 
ments  were  available  at  an  exact  latitude  increment,  the  density  was  inter¬ 
polated  between  the  nearest  two  measurements  on  either  side,  provided 
neither  measurement  was  more  than  2.5°  from  the  desired  position.  Orbital 
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characteristics  were  such  that  the  altitude  was  nearly  constant  at  each 
latitude  throughout  the  flight. 

The  following  form  was  fitted  to  the  data  at  each  latitude 

p(z,t,4>)  =  A(4>)  +  B(4>)g(z,t  -  T,ap)  (1) 


where 

P  =  mass  density 
z  ::  altitude 
<j>  =  latitude 

t  =  time  of  measurement 

t  =  time  delay  between  changes  in  the  planetary  magnetic  index 
a^  and  changes  in  p. 

All  time  variations  in  p  were  absorbed  into  g,  and  the  form  of  both 
g  and  t  were  determined.  A  and  B  were  assumed  to  be  constant,  and  several 
forms  of  g  were  tried;  the  form  that  gave  the  largest  correlation  coefficient 
between  g  and  p  was  considered  "best."  After  selecting  the  best  form  of  g, 
it  was  found  that  to  first  order,  p  had  no  significant  latitudinal  dependence; 
the  apparent  latitudinal  changes  in  A  and  B  could  be  more  accurately  ascribed 
to  changes  in  altitude.  Equation  (1)  was  then  rewritten  as 

p(z,t)  =  A(z)  +  B(z)g(z,t  -  T,ap)  (2) 

To  proceed  further,  an  analytic  form  was  established  for  A  and  B  by 
making  the  following  physically  reasonable  assumptions: 

1.  If  there  were  no  disturbances  (g  =  0),  the  atmosphere  would 
obey  hydrostatic  constraints 

2.  If  g  -  0,  the  atmosphere  is  in  diffusive  equilibrium  above 
some  reference  level  (arbitrarily  chosen  to  be  120  km) 

3*.  If  g  =  0,  mean  free  paths  are  sufficiently  long  that  the  tem¬ 
perature  of  any  one  constituent  is  constant  above  the  reference 
level 
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4.  Atoms  are  warmer  than  molecules  since  they  are  created 
with  excess  kinetic  energy.  (Detailed  theoretical  calculations 
of  the  temperature  difference  have  been  made.  ) 

5.  To  a  first  approximation,  the  quiet-time  atmosphere  between 
120  and  220  km  can  be  considered  to  be  composed  of  only 
oxygen  atoms  and  molecules  with  a  mass  of  about  29  amu 

6.  During  disturbed  times  (g  i  0),  several  forms  of  atmospheric 
heating  may  be  operative;  joule  heating  and  dumping  of  high 
energy  (2  to  10  keV)  particles  from  the  exospheric  ring 
current  (Refs.  24,  25)  are  examples. 

These  six  conditions  impose  the  following  form  on  A  and  suggest  the 
form  of  B  that  will  be  used 


A(z)  =  a^  exp  [-Q^(z  -  120)]  +  a.^  exp  [~a^(z  -  120)]  (3) 

B(z)  =  b j  exp  [-Pj(z  -  120)]  +  b 2  exp  [-  P2(z  -  120)]  (4) 


It  was  expected  that  the  values  of  a^  and  b^  would  change  when  ap 
exceeded  some  critical  value,  indicating  that  the  ring  current  had  become 
unstable  and  was  dumping  particles  into  the  atmosphere  (Ref.  6). 


C.  RESULTS 


The  numerical  results  reported  in  this  section  were  derived  from 
low  latitude  data  obtained  during  the  entire  LOGACS  flight.  Discussions  of 
extensions  of  the  model  and  applications  to  future  flights  are  presented  in 
Section  D  below. 


1.  FORMS  OF  g  AND  t 


in 

is 


The  change  in  density  at  midlatitudes  correlates  well  with  changes 
raluated  at  a  time  t  earlier.  In  p 
greater  than  0.  9  when  g  is  of  the  form 


ap,  evaluated  at  a  time  t  earlier.  In  particular,  the  correlation  coefficient 


g(T, ap)  =  ln[ap(t  -  t)  +  1] 


(5) 
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and  t  is  represented  by 


42.759  -  1.406  (z  -  120) 
5.5 

0.0625  +  0.  145{z  -  120) 
7.9 


130.  0  <  z  <  146.  5 
146.  5  <  z  <  157.  5 
157.  5  <  z  <  174.  0 
174. 0  s  z 


(6) 


where  t  is  measured  in  hours,  and  z  in  kilometers;  the  decimal  places  have 
been  carried  only  to  ensure  continuity. 

The  altitude -dependent  time  delay  is  of  a  form  consistent  with 
localized  heating  near  140  km  (Ref.  17),  and,  when  averaged  over  the  entire 
flight,  is  accurate  to  about  10  percent,  which  does  not  justify  an  expression 
more  complicated  than  this  linear  form. 

2.  DENSITY  COEFFICIENTS 

Equation  (2),  using  the  functional  forms  given  by  Eqs,  (3)  through  (6), 

was  used  to  fit  the  atmospheric  density  measured  on  all  orbits  available 

(about  40  were  useful  in  the  low  altitude  regime).  The  residuals  showed  a 

sudden,  apparently  discontinuous,  increase  when  g(T, a  )  become  greater  than 

P 

about  4.  0;  therefore,  the  data  were  divided  into  two  groups.  One  group 
included  the  densities  measured  during  moderately  disturbed  times 
[a^(t  -  t)  <  59];  the  other  group  included  densities  measured  during  more 
severely  disturbed  times  [a^t  -  t)  >  59].  The  numeral  59  is  accurate  to 
about  20  percent.  The  free  parameters  implicit  in  Eq,  (2)  were  then  fitted 
independently  in  each  group.  The  results  are  shown  in  Table  9. 

A  measure  of  the  goodness  of  fit  can  be  obtained  by  forming  the 

quantity 


a 


P.l  P.J 
l  M 


1/2 


(7) 
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Table  9.  Values  of  Adopted  Parameters  for  Present  Model 


-aj(z-120)  -a^z-^O) 

p(z,t)  =  aj  e  +  a2  e 

-P. (z - 120) 

l[bl'  +b2 

-p2(zl20)l 
e  g(z*ap) 

3p(t  -  r  )  s  59,  gm/cm3 

3 

ap(t  -  t  )  >59,  gm/cm 

al 

103. 05  X 10'13 

115. 53  x  10'13 

a2 

12.45  X  10"13 

9.  18  X  10'13 

bl 

6.49  X  10"  13 

13. 56  X  10"13 

b2 

1.79  X  10‘13 

2.  66  x  10" 13 

0.06064  km'1 

°2 

0.02145 

h 

0.08410 

02 

0.02005 

where  p^  is 

th 

the  i  measured  density  and 

is  the  density  calculated  from 

the  model;  N  -  M  is  the  number  of  degrees  of  freedom.  The  ratios 
are  plotted  in  Fig.  55.  The  density  ratios  and  standard  deviations  derived 
from  this  model  for  the  period  of  the  LOGACS  flight  can  be  compared  with 
those  obtained  using  other  model  atmospheres  that  have  been  proposed.  The 
comparison,  both  in  an  overall  sense  and  as  a  function  of  altitude,  is  shown 
in  Table  10. 


r 
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150 


210 


230 


170  190 

ALTITUDE,  km 


Fig.  55.  Ratio  of  Measured  Densities  to  Those  Predicted  by  Present 
Model  for  Entire  Flight  at  Low  Latitudes 


4-6 


Table  10.  Comparison  of  Present  Model  with  Four  Others  During  LOGACS  Flight 
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Current  operational  model,  essentially  Jacchia  1960 

The  Bruce  modification  to  Walker's  analytical  Jacchia  model  (Ref.  6) 

See  Ref.  27.  Fractional  density  changes  and  time  delay  at  150  km  have  been  applied  at 
all  altitudes 


D. 


DISCUSSION 


1'ho  U.S.  Standard  1962  atmosphere  is  completely  static,  but  the 
"L-DENSIT  Y"  model  has  a  variation  due  to  change  in  From  results 

shown  in  Table  10,  for  these  two  models  alone,  it  would  appear  that  some 
type  of  Fjq  ^  variation  is  useful  above,  but  not  below,  about  170  km.  The 
present  model  has  no  short-term  radio  flux  variation  designed  into  it,  yet  it 
has  a  standard  deviation  of  the  residuals  considerably  less  than  that  of 
L-DENSITY;  thus,  the  apparent  short-term  Fjq  ^  dependence  in  previous 
models  is  not  sufficient  to  completely  reflect  actual  changes  in  density. 

The  J-W-B  model  eliminates  the  density  inversion  near  150  km  that 
troubles  the  more  recent  Jacchia  models,  but  retains  the  exospheric  tempera¬ 
ture  variations  with  Fjq  local  time,  magnetic  activity,  and  other  charac¬ 
teristics  of  the  Walker -Jacchia  formulation  (Ref.  6).  A  dramatic  improvement 
at  low  altitudes  is  the  result. 

The  Champion  model  (Ref.  27),  as  used  here,  introduces  a  constant 
time  delay  between  geomagnetic  variations  and  the  resulting  density  changes. 
The  magnetic  disturbance  not  only  increases  the  exospheric  temperature,  but 
also  increases  the  calculated  density  by  an  empirically  derived,  disturbance - 
dependent  fraction  at  all  altitudes.  The  improvement  of  a  results  primarily 
from  the  time  delay. 

The  present  model  differs  in  several  ways  from  the  other  four 
models  of  Table  10.  It  is  not  simply  a  modification  to  an  existing  model,  but, 
rather,  is  an  attempt  to  "start  from  scratch,"  incorporating  the  relevant 
physical  processes  that  we  believe  dominate  the  upper  atmosphere.  The  two 
most  significant  results  of  this  approach  are  (1)  incorporation  of  an  altitude- 
dependent  time  delay,  and  (2)  recognition  that  the  atomic  and  molecular 
constituents  of  the  atmosphere  do  not  necessarily  share  a  common  tempera¬ 
ture.  The  latter  makes  the  reasoning  behind  the  present  one-parameter 
models,  i.e.,  Tffl,  somewhat  suspect  at  low  altitudes  where  the  number 
density  of  molecules  becomes  significant  relative  to  that  of  the  atoms. 
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The  present  model  is  not  restricted  by  the  artificial  constraint  that 
the  boundary  conditions  at  120  km  remain  constant;  in  fact,  the  density  at 
120  km  can  vary  by  as  much  as  a  factor  of  2  during  a  severe  magnetic  storm, 
Krankowsky,  et  al.  (Ref.  28)  observe  similarly  large  variations  near  120  km 
on  a  day-to-day  or  day-to-night  basis,  even  in  the  absence  of  exceptionally 
great  magnetic  activity.  Density  variations  are  minimal  in  this  model  near 
150  km,  and  even  there  the  variations  can  be  as  much  as  80  percent. 

1.  LIMITATIONS  OF  MODEL 

The  form  of  the  present  model  is  believed  adequate  to  describe  and 
predict  density  changes  in  the  atmosphere  at  midlatitudes  between  140  and 
200  km.  It  must  be  emphasized,  however,  that  the  numerical  values  of  the 
coefficients  apply  to  conditions  during  the  LOGACS  flight.  In  particular, 
long-term  variations  in  Fjq  ^  and  diurnal  variations  were  not  considered. 

From  an  operational  viewpoint,  this  model,  alone,  is  not  adequate 
to  describe  the  drag  history  of  any  satellite  that  is  at  low  altitudes  above 
about  50°  latitude  for  a  significant  time  period.  From  a  physical  viewpoint, 
the  model  is  incomplete  because  temporal  relaxation  effects  and  the  horizontal 
transfer  of  energy  from  the  polar  and  auroral  regions  have  not  been 
included,  but  will  be  incorporated  later. 

The  quiet-time  temperatures  derived  from  the  model  are  near  810° 
and  320° K  for  the  atomic  and  molecular  constituents,  respectively  (a  mass 
of  29  was  used  for  the  molecules  to  correct  for  the  oxygen  actually  present). 
Even  for  extremely  undisturbed  times,  a  300° K  temperature  difference  is 
surprisingly  large,  and  brings  into  question  the  original  functional  form.  If 
assumption  3  (Section  B),  that  the  temperature  of  any  one  constituent  is 
constant  above  the  reference  level,  is  replaced  by  the  assumption  that  all 
constituents  have  the  same  temperature  at  a  given  altitude  but  that  the 
temperature  profile  is  altitude  dependent,  another  plausible  form  for  the 
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model  can  be  generated.  In  particular,  a  temperature  profile  of  the  form 
shown  in  Eq.  (8)  is  being  used  to  fit  the  data 

T  =  T m  1 1  -  a  exp  (  (z  •  zQ)/b]J  (8) 

where  T^  ,  a,  and  b  are  constants,  and  z  is  the  reference  altitude.  The 
results  of  this  fitting  procedure  will  be  discussed  in  Section  V. 

2.  EXTENSIONS  OF  MODEL 

The  correlation  between  changes  in  density  and  changes  in  solar 
radio  flux  is  well  known.  Presumably,  a  change  in  radio  flux  is  accompanied 
by  a  change  in  solar  EUV  flux;  therefore,  the  heat  input  to  the  atmosphere 
should  be  similar  to  that  due  to  particle  dumping.  Equation  (5)  can  be 
rewritten  as 


g(T,a  ,  F)  =  ln[a  (t  -  r  )  +  1]  +  0.  027(F  -  183)  (9) 

P  P 

where  F  is  the  five-day  average  of  the  Ottawa  10.  7  cm  solar  flux  ending  on 

the  day  for  which  p  is  to  be  calculated;  F  is  measured  in  solar  flux  units  of 
22  -  2  -  1 

10-  W-m  -Hz  .  The  numerical  factors  were  chosen  to  agree  with  the 
average  F  observed  during  the  LOGACS  flight,  and  to  match  the  density 
variations  of  CIRA  1965  (Ref.  28)  near  220  km. 

Diurnal  effects  must  eventually  be  included.  Several  investigators 
have  measured  these  changes  (Refs.  29  and  30),  but  the  large  discrepancies 
make  more  complete  sampling  necessary  before  they  can  be  included 
quantitatively.  Other  effects,  such  as  horizontal  energy  transport  and  the 
latitude  dependence  of  the  solar  "bulge,  "  will  be  included  when  a  more 
sophisticated  model  becomes  necessary. 

An  independent  model  applicable  to  the  polar  and  auroral  regions 
is  being  developed.  It  will  be  joined  to  this  midlatitude  model  in  a  future 
report. 
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3. 


EXAMPLES  OF  MODEL 


If  the  altitude -dependent  time  delay  is  ignored  (or,  alternatively,  if 

the  value  of  a^  has  been  constant  for  about  24  hr),  the  present  model  can  be 

compared  with  other  models  currently  in  use.  Figure  56  shows  the  total  mass 

density  plotted  against  altitude  for  a  “  0,  20,  59,  and  200,  calculated  from 

P 

the  present  model  for  F  =  175  flux  units;  and  the  density  plotted  from  the 

CIRA  1965  model  for  a^  5  and  200.  When  compared  with  the  CIRA  model, 

the  major  characteristics  of  this  new  model  arc  immediately  apparent: 

variable  density  at  120  km,  generally  higher  densities  between  140  and  220 

km,  and  much  larger  variability  with  a^. 

Figure  57  shows  how  p,  calculated  from  the  extended  model,  depends 

on  F,  when  F  varies  from  125  to  200  flux  units  for  a  =5. 

P 

Figure  58  shows  the  response  of  the  atmosphere  to  an  impulsive 

change  in  a  ,  as  reflected  by  the  present  model.  An  increase  of  a  from 
P  P 

zero  to  100  is  made  to  occur  from  time  t  =  0  to  t  =  0.  5.  The  density  of  the 

atmosphere  is  unaffected  until  t  =  5.  5  hr,  when  a  "bump"  appears  between 

147  and  157  km.  The  bump  then  splits  into  two.  The  high-altitude  bump 

moves  to  successively  higher  altitudes,  and  becomes  increasingly  broad 

until,  by  t  =  8  hr,  the  whole  atmosphere  above  172  km  is  disturbed.  Between 

146  and  171  km,  however,  the  density  has  already  returned  to  its  quiet-time 

value;  by  t  =  8.  5  hr,  the  density  above  171  km  has  also  returned  to  its 

prestorm  value.  The  low-altitude  bump  moves  very  slowly  toward  lower 

altitudes,  and  becomes  somewhat  narrower;  by  t  =  8.  5  hr,  this  bump  is  the 

only  indication  that  a  storm  has  occurred. 

The  traveling  density  bulges  are  the  result  of  a  slowly  moving 

disturbance,  accompanied  by  mass  flow,  propagating  away  from  the  150  km 

region  where  the  initial  heating  occurred.  The  existence  of  density  bulges 

considerably  complicates  the  interpretation  of  any  measurements  of 

atmospheric  density.  For  example,  in  the  extreme  case  shown  in  Fig.  58, 

a  measurement  of  density  scale  height  made  vertically  upward  at  t  =  7.  5  hr 

would  show  the  scale  height  to  be  reasonable  up  to  about  145  km,  then  vary 
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Comparison 
of  Present  1\ 


PRESENT  MODEL,  F  «  175 
CIRA  1965,  MODEL  6,  1000  hr 
CIRA  1965,  MODEL  10,  1000  hr 


tween  Steady-State  Density  Prediction 
el  and  That  of  CIRA  1965 


DENSITY,  gm/cm 


i 


wildly  from  positive  to  negative  through  minus  infinity  and  back  to  positive, 
all  in  the  space  of  a  few  kilometers;  it  would  then  settle  down  until  an  altitude 
of  about  168  km  was  reached,  then  repeat  the  same  gyrations,  although  some¬ 
what  less  violently.  The  model  predicts  that  negative  scale  heights  rarely 
occur  in  the  absence  of  particle  dumping;  therefore,  for  small  disturbances 

(a  <  40),  density  measurements  would  indicate  only  that  there  were  regions 
P 

where  the  deduced  scale  height  became  very  large,  then  suddenly  small 

before  returning  to  a  more  normal  value. 

An  example  of  the  density  response  to  a  very  moderate  magnetic 

storm  is  shown  in  Fig.  59.  It  is  assumed  that  a  was  zero  until  time  t  =  0 

P 

when  it  linearly  increased  to  100  in  a  period  of  3  hr,  remained  stable  for 
another  3  hr,  then  linearly  decayed  to  zero  in  the  succeeding  15  hr.  The 
density  is  shown  as  a  function  of  altitude  for  t  =  5  to  10  hr.  Again,  the  effect 
of  the  traveling  disturbances  is  evident.  Since  previous  model  atmospheres 
cannot  reproduce  this  very  complicated  response  of  the  atmosphere  to  even 
moderate  storms,  it  is  to  be  expected  that  this  new  model  will  allow  much 
more  accurate  and  detailed  prediction  of  the  drag  force  experienced  by  low- 
altitude  satellites. 


4 


1 
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Fig.  59.  Response  of  Present  Model  to  Realistic 
Moderate  Magnetic  Storm 


SECTION  V 


ALTERNATE  FORMULATION  OF  MID  LATITUDE  MODEL 

A.  INTRODUCTION 

Section  IV  proposed  a  novel  physical  explanation  for  the  midlatitude 
results  of  the  LOGACS  experiment,  namely,  that  during  magnetically  quiet 
times  each  individual  constituent  of  the  atmosphere  has  its  own  characteristic 
temperature.  It  was  found  that  atomic  species  were  several  hundred  degrees 
warmer  than  molecular  species.  Since  atomic  oxygen  is  created  with  a 
considerable  excess  of  kinetic  energy,  it  must  be  somewhat  hotter  than  either 
molecu’ar  oxygen  or  nitrogen,  but  the  temperature  difference  necessary  to 
fit  the  data  was  surprisingly  large. 

Accordingly,  the  exact  opposite  assumption,  that  all  constituents 
share  a  common  altitude -dependent  temperature,  is  made  in  this  section. 

The  other  characteristics  of  the  previous  model  are  retained;  most  impor¬ 
tantly,  the  altitude -dependent  time  delay  between  variations  in  the  magnetic 
field  and  the  corresponding  variations  in  atmospheric  density  has  been  left 
unchanged. 

B.  TEMPERATURE/DENSITY  PROFILE 

The  following  temperature  profile  was  assumed 

T  =  Tjl  -  a  exp  [-(z  -  (10) 

where  T  is  the  temperature  at  any  altitude  z,  T^,  is  the  exospheric  temperature, 
a  and  H^,  are  positive  constants,  and  z^  is  the  reference  altitude.  At  z^,  one 
finds 

T=  TQ  =  (1  -  a)Ttt  (11) 
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Assuming  diffusive  equilibrium,  it  can  be  shown  that  the  number 
density  of  any  one  constituent  ic  given  by 

1+Ht/Hm 

n  =  n0(TQ/T)  exp  [ - (z  -  zQ)/HM]  (12) 

where  ng  is  the  number  density  at  Zg  and  is  the  exospheric  scale  height. 
We  have  defined 

I^skT./mg  (13) 

where  k  is  Boltzmann's  constant,  m  is  the  mass  of  the  constituent,  and  g  is 
the  acceleration  due  to  gravity  (assumed  constant  between  120  and  250  km 
where  the  model  applies). 

This  temperature  profile  can  be  made  to  fit  that  used  by  CIRA  by 
adjusting  the  values  of  a,  Rj,,  and  Tq,  and  is  of  the  same  form  as  that  which 
Jacchia  assumes.  It  has  several  properties  which  are  very  useful  in  a  wide 
variety  of  atmospheric  problems.  The  temperatures  gradient,  given  by 

If  *  '  V/ht>  <14> 

is  everywhere  positive,  and  both  its  magnitude  and  curvature  are  adjustable 
by  varying  a  and  H^,.  The  temperature  profile  makes  the  transition  to  the 
isothermal  case  smoothly  as  a  — 0,  and  there  is  no  discontinuous  change  of 
form  in  the  number  density  profile.  The  number  density  of  any  one  constituent 
in  a  vertical  column  above  some  height  z  is  simply  given  by 


N=  nH 


(15) 
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where  n  is  evaluated  at  z  and 


H  -  H(z)  =  kT(z)/mg 
Similarly,  the  pressure  due  to  that  constituent  is 


(16) 


p  =  nkT  =  pQ 


(1  -  a)  exp  [-(z  -  zq)/Ht] 
1  -  a  exp  f-(z  -  zq)/Ht] 


ht  /hm 


(17) 


which  is  useful  in  employing  the  analytical  isobaric  models  of  Thomas  and 
Ching  (Refs.  7  and  31),  for  instance. 

The  total  mass  density  is  given  by 


P=£  M.n. 

~  li 
i 


(18) 


p  =  <To/T>  KoI^o.n^1'75  +po,o/] 

where  pQ  Q,  pQ  and  pQ  ^  are  the  mass  densities  of  the  three  major 

constituents  at  the  120  km  reference  level,  and 


q  =  (TQ/T)T  16  exp  [-(z  -  zq)/H16]  (20) 

where  is  the  exospheric  scale  height  of  atomic  oxygen. 

C.  METHOD  OF  ANALYSIS 

The  model  derived  in  the  earlier  section  represented  the  density  at 
any  altitude  by  a  sum  of  four  exponentials  with  scale  heights  independent  of 
the  level  of  magnetic  activity;  the  preexponential  factors  were  functions  of  the 
magnetic  index  a  and  altitude  only. 
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In  deriving  the  present  model,  the  time -independent  densities  from 

the  previous  model  were  used  as  inputs  to  a  least- squares  routine  which 

sought  the  best  values  of  Tq,  a,  H,p,  and  the  individual  mass  densities  at  the 

reference  level  pn  ~  ,  and  prt  for  eleven  values  of  a  between 

K0,O’  m0,C>2  0,N2  p 

zero  and  400.  The  fits  in  each  case  showed  standard  deviations  of  less  than 
2  percent  in  the  logarithmic  density  ratios,  which  were  well  within  the 
experimental  uncertainties.  Several  runs  were  then  made  allowing  z^  to 
vary  from  its  arbitrarily  assigned  value  of  120  km.  The  results  gave  a 
least  -  squares  fit  to  the  data  that  always  yielded  a  value  of  z^  within  a  few 
kilometers  of  120  km,  and  the  standard  deviation  was  not  significantly 
changed.  Consequently,  z^  was  held  fixed  at  120  km  for  the  computations 
reported  here. 


D. 


RESULTS 


Figures  60  through  62  show  the  quality  of  the  least  -  squares  fits. 
The  points  arc  the  density  values  calculated  from  the  previous  model;  the 
solid  curves  represent  the  best  fit  densities  of  this  new  formulation. 

Figures  63  and  64  show  how  each  of  the  six  fitted  parameters  men¬ 
tioned  above  vary  as  a  function  of  ln(a^  +1).  The  sharp  break  at  4.  1 

(a  =  59)  reflects  the  discontinuous  nature  of  the  prior  model,  which  we  have 
P 

attributed  to  the  onset  of  significant  particle  dumping. 

The  variations  in  the  temperature  parameters  can  be  approximated 
analytically,  as  follows 


x  =  In  (a  +1) 

P 

(21) 

a  =  0.6719  +  0.  027x 

(22) 

TqPK) 

|  536.  3  -  2.  31x 

a 

P 

<  59 

(23a) 

<  505.  5  -  2.  03x 

a 

p 

>  59 

(23b) 

H  ( km ) 

(305.  6  -  2.  87 x 

r 

a 

p 

<  59 

(24a) 

1 

<308.  0  -  2.  98x 

a 

p 

>  59 

(24b) 

Tb(°K) 

=  (1612  +  175. 4x 

r 

a 

p 

<  59 

(25a) 

<1658  +  204. 6x 

r 

a 

>  59 

(25b) 

P 
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Fig.  62,  Mass  Density  vs  Altitude  for  ap  =  200,  300,  and  400, 
(Points  taken  from  Section  IV;  solid  curves  are 
best-fit  solutions  reported  here.  ) 
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Fig.  63.  Temperature  Distribution  Parameters  Ht#  Tq,  and  a. 
[Solid  lines  represent  Eqs.  (22),  (23),  and  (24).  ] 
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01  2345678 

Ln(op+ 1) 

Fig.  64.  Mass  Densities  of  N2,  O2,  and  O  at  120  km.  (Solid 
lines  represent  Eqs.  (26),  (27),  and  (28).  ] 
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The  densities  of  the  three  individual  constituents  O,  O^,  and  N2 
at  the  120  km  reference  level  vary  in  much  the  same  way 


p0>o(l°-12gm/cm3) 


p0,N2(l0'12gm/cm3) 


|  2.  537  +  0.  286x 
I  2.  405  +  0.449x 

|  6.  173  +  0.  4 1 3x 
I  6.  447  +  0.  972x 


a  <59 
P 

a  >59 
P 


a  <  59 
P 

a  >59 
P 


(26a) 

(26b) 

(27a) 

(27b) 


’0,o  (10  12  grn/cm3) 


|  3.  103  +  0.  201x 
1 4.  160  +  0.  273x 


a  2  59 
P 

a  >59 
P 


(28a) 

(28b) 


Equations  (22)  through  (28)  were  derived  using  a  linear  least- squares  fit. 
Equations  (22),  (23),  and  (24)  are  shown  as  solid  lines  in  Fig.  63;  Eqs.  (26), 
(27),  and  (28)  are  similarly  shown  in  Fig.  64. 

1.  DISCUSSION 

It  is  much  more  convenient  to  use  the  model  form  given  in  Section  IV 
(a  sum  of  four  exponentials)  to  calculate  the  total  density  at  any  given  altitude 
than  it  is  to  use  the  present  formulation.  However,  the  sum  of  exponentials 
cannot  easily  be  used  to  obtain  any  information  on  the  detailed  behavior  of 
the  individual  components  of  the  atmosphere,  nor  is  it  adequate  to  calculate 
the  pressure,  which  is  necessary  to  determine  the  dynamical  response  of  the 
atmosphere  to  changes  in  the  heating  rate. 

The  present  form  of  the  midlatitude  model  is  constructed  so  that  it 
reproduces  the  same  densities  and  density  variations  as  the  previous  one,  and 
in  addition,  makes  explicit  the  relationships  among  all  the  temperatures  and 
composition  parameters. 
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Figure  65  shows  the  steady -state  temperature  profile  for  five 
•'alues  of  a^,  as  calculated  by  the  present  model.  A  quiet  and  a  disturbed 
profile  taken  from  CIRA  arc  shown  for  comparison.  At  the  120  km  reference 
level  the  present  model  indicates  that  the  temperature  is  150  to  200°  hotter 
than  CIRAs  fixed  355°  K,  and  decreases  with  increasing  magnetic  activity. 
Above  approximately  150  km,  the  temperatures  indicated  by  this  model 
increase  monotonically  with  a^,  but  at  a  much  slower  rate  than  CIRA  pre¬ 
dicts.  The  temperature  gradient  is  about  half  that  of  CIRAs:  for  a  =0 

P 

and  400  we  have  (dT/dz)  150  km  ~3.  2  and  6.  7° /km,  respectively,  while  for 
the  two  extreme  CIRA  models  the  gradients  are  6.  1  and  13.7°/km.  The 
exospheric  temperatures  predicted  by  the  present  model  are  considerably 
higher  than  CIRAs;  since  this  model  was  derived  from  data  below  220  km 
only,  not  much  credence  should  be  placed  on  extrapolated  values,  including 
T  ,  far  above  about  250  km. 

Figures  66  through  68  show  the  steady-state  altitude  profiles  of 
the  number  densities  of  atomic  oxygen,  molecular  nitrogen,  and  molecular 
oxygen,  respectively,  for  the  same  five  values  of  a^.  The  altitude  dependence 
of  all  three  constituents  in  this  model  is  much  closer  to  being  exponential 
than  is  the  CIRA's  dependence,  which  is  a  consequent  of  the  smaller  tem¬ 
perature  gradient  in  this  model.  The  number  density  of  molecular  nitrogen 
is  fairly  close  to  that  predicted  by  CIRA,  although  below  140  km  it  is  some¬ 
what  less.  Both  atomic  and  molecular  oxygen  have  densities  that  arc  con¬ 
siderably  larger  than  the  CIRA  densities.  None  of  the  constituents  in  the 
present  model  exhibit  the  inversion  below  about  180  km  that  is  a  characteristic 
of  fixing  the  total  mass  above  a  reference  level,  as  is  the  practice  in  most 
other  models. 

Tables  11  through  21  list  the  temperature  (°K),  mean  molecular 

_3 

mass  (amu),  number  densities  of  O,  N^,  and  (cm  ),  and  pressure  scale 
height  (km)  every  5  km  between  120  and  220  km  for  a^  =  0,  5,  10,  20,  40,  59, 
60,  100,  200,  300,  and  400.  The  discontinuous  change  in  the  model  between 
a^  =  59  and  60  is  an  artifact  representing  the  onset  of  particle  dumping; 
LOGACS  saw  a  large  change  in  density  near  a  =  60,  but  the  data  were  not 
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Fig.  66.  Number  Density  Profiles,  n,  of  Atomic  Oxygen  for 
ap  =  0,  10,  40,  60,  and  400.  (Dashed  curves  are 
extreme  CIRA  number  densities. ) 


cm' 


Fig.  67.  Number  Density  Profiles,  n,  of  Molecular  Nitrogen  for 
ap  =  0,  10,  40,  60,  and  400.  (Dashed  curves  are 
extreme  CIRA  number  densities.  ) 
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Table  11.  Atmospheric  Parameters  for  F  =  183 
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Table  12.  Atmospheric  Parameters  for  F  =  183 
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Table  13.  Atmospheric  Parameters  for  F  =  183 
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Table  15.  Atmospheric  Parameters  for  F  =  183 
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Table  16.  Atmospheric  Parameters  for  F  -  183 
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Table  17.  Atmospheric  Parameters  for  F  =  183 
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Table  18.  Atmospheric  Parameters  for  F  18%  a  100 
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Table  19.  Atmospheric  Parameters  for  F  =  183,  a  200 
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Table  21.  Atmospheric  Parameters  for  F  =  183,  a  =  400 
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sufficient  to  allow  a  more  realistic  (i.  e.  ,  smoother)  representation  than 
has  been  used  here. 

2.  ENERGY  CONSIDERATIONS 

A  considerable  quantity  of  energy  is  required  to  change  the 

atmosphere  from  its  quiet -time  distribution  to  one  which  corresponds  to 

2 

values  of  a^  greater  than  zero.  The  total  energy  content  per  cm  above 
120  km  is  given  approximately  by 


E(a 


3  r 

^  I  fm.gz  +  1.  5  kT(a^)]n.(ap) 


dz 


i=l  120 


(29) 


where  the  index  runs  over  the  three  major  atmospheric  constituents. 

Figure  69  shows  the  quantity  E(a  )  -  E(0)  plotted  against  a  ;  the  dashed 

P  P 

line  represents  a  best -fit  curve 


E(a  )  -  E(0)  =  1.425  +  0.  658  Ja~ 
P  v  P 


(30) 


and  the  solid  line  represents  a  best-fit  curve  to  the  next  higher  order 


E(a  )  -  E(0)  =  0.  015  +  1.  142  /a~-  0.024  a 
P  P  P 


(31) 


4  2 

where  E  is  measured  in  units  of  10  erg/cm  . 

Thomas  and  Ching  (Ref,  17)  find  that  the  total  heat  input  at  150  km 
in  Jacchia's  models  (Ref.  5)  is  proportional  to  although  there  is  some 

scatter.  We  find  that  in  this  model  the  total  energy  change  in  the  upper 
atmosphere  can  be  represented  as  being  proportional  to  sfaT>  [Eq.  (30)],  but 
the  representation  is  roughly  three  times  better  when  the  expansion  is 
carried  to  the  next  higher  order,  as  in  Eq.  (31). 
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Of 


Fig.  69.  Left  Ordinate:  Energy  Content  of  Atmosphere  Above  120  km. 

Right  Ordinate:  Minimum  Energy  Deposition  Rate  Above  120  km. 
[Points  calculated  from  model;  dashed  and  solid  curves 
represent  Eqs.  (30)  and  (31),  respectively.] 
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The  rate  of  energy  deposition  cannot  be  obtained  directly,  since 

loss  processes  have  not  been  included.  It  is  possible,  however,  to  find  a 

lower  limit  to  this  rate.  According  to  the  time  delays  built  into  the  model 

(Section  IV),  it  takes  7.9  hr  for  the  upper  atmosphere  to  become  completely 

accommodated  to  a  step  change  in  the  value  of  a^.  Assuming  that  the  energy 

deposition  were  constant  over  this  entire  time  and  that  a^  changed  from  zero 

to  a  and  then  remained  constant,  the  average  energy  deposition  rate  is 

Pi  max  *“* 


dE 

dt 


0.  396/ 


a 

p,  max 


0.  0083  a 

p,  max 


erg/cm  -sec 


(32) 


This  is  the  same  form,  up  to  a  constant  term,  as  in  Eq.  (31);  the  scale  on 
the  right  side  of  Fig.  69  refers  to  this  minimum  average  energy  deposition 
rate. 

The  change  in  upper  atmospheric  energy  with  changes  in  a^  applies 
specifically  to  this  particular  model,  which  is  valid  for  Flf.  _  =  183  flux  units 

22  2  i  v#  i 

(10"  erg/m  -Hz)  at  about  1030  local  time.  It  was  suggested  in  Section  IV 
that  the  response  of  the  atmosphere  to  a  variation  of  FjQ  y  from  183  could  be 
treated  as  being  similar  to  its  response  to  a  variation  in  ln(a^  +  1) 

A?10>7  *  37  X  A[ln(ap  +  1)]  (33) 

This  relationship  is  a  consequence  of  requiring  that  the  densities  predicted  by 
this  model  vary  with  Fjq  y  in  the  same  way  as  do  the  CIRA  densities  -tt  about 
220  km. 

If  Eq.  (3  3)  holds,  we  can  obtain  some  knowledge  of  the  variation 
in  solar  EUV  flux  with  changes  in  F^q  y,  since  by  assuming  the  equality  in 
Eq.  (32)  and  combining  it  with  Eq.  (33)  we  find 

=  [2.  18  -  0.  230Q  -  (4.723  -  Q)1/2]/[-  8.  51  +  18.  5(4.  723  -  Q)1/2] 

AF 

(34) 
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where  Q  is  the  energy  input  in  erg/cm  -sec.  The  value  AQ/AF  varies  from 

0.  014  for  Q  =  1.  0  to  a  maximum  of  0.  032  at  Q  =  3.  5,  then  decreases  with 

increasing  energy  input.  During  solar  quiet  times,  Hinteregger  and  Hail 

2 

(Ref.  32)  find  that  the  solar  flux  is  about  3.0  erg/cm  -sec  between  100  and 

O 

1027  A;  this  is  the  wavelength  region  that  is  thought  to  be  most  effective  in 

heating  the  atmosphere  above  120  km.  Examination  of  Eq.  (34)  shows  that  Q 

2 

can  be  no  larger  than  ~4.  7  erg  'cm  -sec.  If  this  value  is  adopted  for  times 

of  maximum  solar  radio  flux  (F  =  225)  and  Eq.  (34)  is  used  to  determine  Q 

at  F  =  100,  we  find  Qjqq  =  1.45  erg/cm  -sec.  Comparing  this  with  the  3.0 

erg/cm^-sec  incident,  the  heating  efficiency  is  -50  percent.  The  ratio  of 

active-to-quiet-time  Q  is  about  3.2,  Since  even  during  solar  quiet  times  the 

atmosphere  is  optically  thick  in  this  wavelength  region,  if  we  assume  that  eddy 

diffusion  is  constant  or  small,  then  the  heating  efficiency  should  be  nearly 

constant.  This  strongly  suggests  that  the  solar  EUV  flux  itself  varies  by  more 

than  a  factor  of  three  when  T  jq  ^  varies  from  100  to  225  flux  units. 

This  conclusion  is  supported  by  some  results  from  OSO-1,  reported 

by  Bourdeau,  et  al.  (Ref.  34).  They  found  that  over  a  solar  rotation,  when 

Flr.  _  changed  from  about  80  to  125,  the  solar  emission  intensity  between  170 
l  u.  <  0  2 

and  370  A  changed  from  0.37  to  0.55  erg/cm  -sec,  yielding 

(1/Q)AQ/AF  =  0.  009.  Over  the  wavelength  range  included  in  Eq.  (34),  we 

find  that  for  F._  _  =  100,  (1/Q)AQ/AF  =  0.014  is  in  remarkably  good 

1  vi  * 

agreement. 

Hinteregger  and  Hall  (Ref.  34)  have  reported  some  measurements 

of  the  variations  with  F.n  _  of  the  intensities  of  six  solar  lines  between  304 
o  1  u.  r 

and  1025  A.  The  measurements  were  made  between  the  middle  of  April  and 
the  middle  of  June  1967,  a  period  which  included  the  time  of  LOGACS  flight 
when  F1q  ^  reached  a  maximum  of  209.  The  maximum  value  of  the  previous 
solar  rotation  was  135;  on  the  following  rotation  it  was  132.  The  minimum 
value  of  Fjq  ^  during  the  period  reported  was  105  in  June  and  108  in  early 
May.  The  energy  average  of  the  ratio  of  the  maximum  EUV  intensity  near 
the  26  May  maximum  to  the  average  of  the  minima  before  and  after  was  about 
1.  5,  and  (1/Q)AQ/AF  =  0.  004.  This  is  less  than  was  observed  by 
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Bourdeau,  et  al. ,  in  1962,  and  considerably  less  than  would  have  been 
predicted  by  the  present  model.  The  differences  may  be  due  to  the  differing 
wavelength  intervals  used  in  the  three  calculations,  or  it  may  be  that  the 
anomalously  high  values  of  F^q  ^  in  the  middle  of  May  were  not  related  to 
changes  in  the  EUV  flux  in  the  same  way  as  the  more  normal  changes  in 
Fjq  ^  during  most  solar  rotations.  It  is  interesting  to  note  that  if  the  maxi¬ 
mum  Fjq  y  during  the  rotations  preceding  and  following  the  26  May  event 
were  used,  (1/Q)AQ/AF  =  0.  017  is  in  much  better  agreement  with  the  other 
two  determinations.  Since  the  maximum  solar  EUV  was  measured  on  26  May, 
however,  this  apparent  agreement  may  be  fortuitous. 
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SECTION  VI 


AURORAL  AND  POLAR  ATMOSPHERE  DENSITY  MODEL 
BETWEEN  120  AND  175  km 


A.  INTRODUCTION 

The  previous  two  sections  described  a  density  and  composition 
model  of  the  upper  atmosphere  that  was  valid  between  120  and  220  km,  and 
that  applied  only  to  midlatitude  regions.  This  section  presents  a  new  density 
model  that  is  designed  to  represent  the  polar  atmosphere.  Since  both  of 
these  models  have  been  developed  using  the  LOGACS  density  data  exclusively, 
the  inputs  have  been  limited  to  a  single  4-day  period.  The  actual  use  of  these 
models  under  a  variety  of  other  conditions  will  probably  lead  to  some  modifi¬ 
cation  of  the  values  of  the  parameters  quoted,  but  it  is  hoped  that  the 
analytical  form  will  remain  an  accurate  representation  of  the  density,  and 
density  variations,  in  the  upper  atmosphere. 

B.  DATA  REDUCTION 

The  only  LOGACS -derived  density  data  used  in  this  model  were  those 
obtained  on  the  "day"  side  of  the  orbit  between  ~87  and  35° N  geocentric 
latitude.  The  data  were  first  used  in  an  interpolation  routine  to  estimate  the 
atmospheric  density  every  5°  of  geomagnetic  latitude  between  75°  and  40°N; 
no  estimate  was  made  unless  there  was  at  least  one  measured  density  value 
on  either  side,  and  within  2.  5°,  of  the  desired  latitude.  There  were  311 
interpolated  densities  so  generated;  these  were  used  in  the  initial  model 
calculations. 

The  parameters  of  an  analytic  density  model  similar  to  that  used  in 
the  midlatitude  case  (a  sum  of  exponentials)  were  adjusted  to  fit  the  data  in  a 
least-squares  sense  at  each  latitude  independently.  The  temporal  variation 
in  the  densities  was  assumed  to  be  correlated  with  magnetic  activity,  and 
various  indices  (a^,  ^Ygie^son»  ^VThule^  an^  ^me  delays  were  used  to  find 
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the  combination  that  best  fit  the  interpolated  densities.  No  combination  gave 
a  very  convincing  fit. 

The  original  density  data  were  again  used  in  the  interpolation 
routine,  this  time  generating  densities  at  each  5*  increment  of  geocentric 
latitude.  There  were  321  such  points  between  85°  and  50°N.  The  fits  were 
much  improved;  in  particular,  the  use  of  (AYj^eison/20  +  AY^^g/S)  as  the 
magnetic  index  and  a  time  delay  of  8  hr  gave  the  best  overall  fit. 

The  coefficients  at  each  latitude  were  very  different,  however,  and 
it  was  not  possible  to  find  a  single  set  of  coefficients  that  adequately  fit  the 
data.  A  detailed  examination  of  the  residuals  indicated  two  reasons  for 
this.  During  disturbed  times  there  was  a  marked  latitudinal  variation  in  the 
density,  with  some  longitudinal  asymmetry;  and  during  both  quiet  and  dis¬ 
turbed  times  there  was  a  local  and  very  large  increase  in  density  over  the 
region  near  Baffin  Island.  This  will  be  called  a  "bulge"  for  convenience  in 
discussion. 

Both  of  these  variations  had  to  be  modeled  to  obtain  a  good  fit,  but 
this  meant  the  introduction  of  several  more  free  parameters;  to  minimize  the 
total  number  of  parameters,  the  basic  sum-of-exponentials  model  was  recast 
into  a  slightly  different  form.  Equation  (35)  is  the  sum-of-exponentials 
model 


P  =  a  j  exp  ( -o^Z)  +  a ^  exp  (  -o^7.)  +  (b ^  exp  (- Pj Z)  +  b^  exp  ( - Z ) ]  ln(M  +  1) 

(35) 

The  a.,  b.,  ff.,  and  0.  are  free  parameters,  p  is  the  density  at  time  t,  7.  is  the 
altitude,  and  M  is  the  magnetic  index  evaluated  at  some  time,  t,  earlier. 
Excluding  t,  there  are  eight  parameters.  This  form  was  replaced  by 


p  =  aj  exp  jV^Z  -  °2^2J  +  bl  exP  [^lZ  ‘  ^2^]  +  ^ 


(36) 
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The  number  of  parameters  has  been  reduced  to  six,  but  at  some  loss  of 
generality.  Above  some  altitude,  roughly  o^/ Oj,  this  would  predict  that  the 
density  increases  without  limit  as  Z  approaches  infinity,  so  the  model  must 
be  restricted  to  altitudes  below  what  in  practice  turns  out  to  be  p^/Pj* 

All  the  interpolated  data  were  fit  using  Eq.  (36),  with  the  magnetic 
index  M  being  the  weighted  average  of  the  Eielson  and  Thule  AY  ranges,  and 
using  an  altitude  -  independent  time  delay  of  8  hr.  As  in  the  midlatitude  model, 
two  sets  of  coefficients  were  required:  one  for  quiet  and  moderately  active 
times  and  one  for  extremely  disturbed  times. 

The  influence  of  the  bulge  had  to  be  removed  before  the  latitudinal 
variations  could  be  studied.  The  bulge  seemed  to  be  localized  and  symmetric 
during  quiet  times,  centered  at  about  67°N  and  290°E.  During  more  disturbed 
times  it  moved  northward  and  eastward,  becoming  an  elongated,  somewhat 
ridge -like  structure.  The  easy  way  to  model  such  movement  and  shape  change 
is  to  transform  to  a  coordinate  system  whose  origin  moves  with  the  center  of 
the  bulge  along  one  cartesian  axis,  make  a  symmetric,  magnetic -activity- 
dependent  bulge  in  that  system,  then  transform  back. 

When  this  was  done  a  deep  density  trough  appeared  in  the  residuals 
at  about  68°N  over  Asia  and  the  western  Pacific,  and  at  about  60°N  over 
Canada  and  Europe.  This  suggests  that  geomagnetic  coordinates  would  be 
the  appropriate  system  to  use,  but  the  trough  in  any  one  section  seemed  to 
follow  a  line  of  constant  geocentric  latitude,  with  a  fairly  abrupt  transition 
between  the  sections.  This  is  consistent  with  the  results  of  the  initial  fits, 
where  the  use  of  geocentric  coordinates  gave  a  far  better  overall  fit  than  did 
the  use  of  geomagnetic  coordinates.  This  trough  was  then  modeled  with  a 
simple  cosine  function. 
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C.  POLAR  MODEL 

1.  SYMBOLS  LIST 

The  symbols  applicable  to  the  polar  model  are: 

B  correction  to  the  model  due  to  the  bulge 

F  5 -day  average  of  F^q  ending  on  the  day  of 

interest 

M  magnetic  index,  evaluated  8  hr  prior  to  the  time 

of  interest 

TR  correction  to  the  model  due  to  the  trough 

Z  altitude,  km 

A\g  Eiclson  AY  range 

AY^  Thule  AY  range 

AG'  effective  width  of  the  bulge,  deg 

0  geocentric  colatitude 

Gq  geocentric  colatitude  of  center  of  bulge 

G'  colatitude  in  bulge  system 

3 

P  density,  gm/cm 

4>  geocentric  east  longitude 

4>0  geocentric  east  longitude  of  center  of  bulge 

<i>'  east  longitude  in  bulge  system 

2.  REGION  OF  VALIDITY 

Because  of  the  variability  of  the  polar  atmosphere,  this  model 
applies  only  when  conditions  are  similar  to  those  when  the  LOGACS  flight 
was  made.  Some  extensions  of  the  model  are  suggested,  but  these  must  be 
validated  in  practice  before  too  much  confidence  can  be  placed  in  them. 
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Altitude:  135  £  Z  <  175 


Latitude:  2  50°N 

(0  <  6  <40°) 


Lower  limit  imposed  by  available  data, 
but  can  probably  be  extrapolated  to 
120  km;  upper  limit  imposed  by  the 
form  of  the  model;  for  higher  altitudes 
use  the  midlatitude  model  with  a 
replaced  by  M.  ^ 

For  lower  latitudes,  to  50°S,  use  the 
midlatitude  model. 


Local  time:  ~1030 


F:  ~180 


The  model  presumably  can  be  extended 
to  about  3  hr  on  either  side  fairly 
safely;  no  information  is  available  at 
other  times. 

For  difference  values  of  F,  replace  M 
by  {M  +  exp  [0.  027(F - 183)]  -  1 1  ; 
this  is  in  analogy  to  the  midlatitude 
model. 


Season:  -summer  solstice 


There  are  two  models.  One  includes  the  bulge  and  is  more  accurate 
in  reproducing  the  LOGACS-de rived  densities.  Because  the  bulge  is  so 
variable  in  both  position  and  extent,  and  because  a  badly  modeled  bulge  is 
worse  than  none  at  all,  a  second  model,  which  averages  the  effects  of  the 
bulge,  is  also  presented. 

3.  BULGE  MODEL 

Parameters  for  the  bulge  model  are  evaluated  as  follows: 
a.  Magnetic  Index 


M  =  AYe/20  +  AYt/8 


(37) 


b.  Evaluation  of  B 


If  4>  is  not  between  280  and  340°,  B  =  0. 
If  4>  is  between  280  and  340° 


0  =  23°  -0,  18M 
o 


(38) 
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(39) 


4>o  =  290°  +  0.  18M 

_j  (*  Gsin  (4>  -  115°)  -  6o  sin  (<j>Q  -  115°)1 
=  tan  [ecos  (4>  -  115°)  -  eo  cos  (4>o  -  1 15° )J 

0*  =  [e  cos  (4>  -  115°)  -  0q  cos  (4>o  -  1 15° )]/cos  <t>' 


If  M  <  16.  33 


(40) 


AG'  =  exp  (0.  177M) 


(41a) 


M  >  16.  33 


AG'  =  18° 


If  |  G'|  >  AG' 


B  =  0 


If  |6’  |  <  AG'  and  M  <  59 

B  =  0.  276  cos  [(G7aG')(450)] 
If  Ig’I  <  AG'  and  M  >  59 


B  =  0.  222  cos  [(Q'/AG'  )(4 5  °)  ] 


c.  Evaluation  of  TR 

If  M  <  59,  TR  =  0. 

If  M  >  59,  and  50°  <  d>  <  230° 

TR  =  cos  [7.  9(6  -  22.  5°)] 
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(41b) 


(42a) 


(42b) 


(42c) 


(43a) 


If  M  >  59,  230°  <<$>  <  [(360° )  +  50°  ],  and  6  <  13° 


TR  =  -0.743 


If  M  >  59,  230°  <  4>  <  [(360°)  +  50°  ],  and  6  >  13° 


TR  =  cos  [8.  1(0  -  30°)] 


(43b) 


(43c) 


d.  Evaluation  of  p 
If  M  <  59 


P  =  3.  89  X  10'14 


Z  -  324. s\Z 
87.08  / 


+  I.  635  exp 


(z  -  215. 
V  55.  II 


ln(M  +  1) 


X(1  -  0.0047  M)  X  (1  +  B) 


(44a) 


If  M  >  59 


P  = 


5.605  x  10  *4  exp 


X  [1  +  ln(M  +  1)](1  +  B)[  1  +  0.0396  TR  (1  -  0.0515  M)]  (44b) 

4.  BULGELESS  MODEL 

The  correction  B  is  not  calculated;  TR  is  evaluated  in  the  same  way 
as  for  the  bulge  model,  and  the  coefficients  in  the  p  calculation  are  changed 
slightly. 
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If  M  <  59 


P  = 


3.96  X  10'14 


exp 


-  324.5 
87.08 


+  1.  635  exp 


In(M  +  1) 


X  (1  -  0.0047M) 


(45a) 


If  M  >  59 


P  =  5.  75 


- 14 

X  1 0  exp 


/Z  -  241. 7\2 

[\  59.3  ;  J 


X[1  +  ln(M  +  l)]x[l  +  0.  0442  TR(1  -  0.  0515M)]  (45b) 


D.  DISCUSSION 

Since  only  LOGACS  data  were  used  in  the  development  of  this  polar 
model,  it  is  difficult  to  estimate  its  accuracy  when  applied  to  other  flights. 

If  attention  is  restricted  to  the  LOGACS  vehicle  only,  however,  one  may 
compare  the  rms  density  residuals  obtained  by  using  this  model  with  those 
from  several  other  model  atmospheres.  Such  a  comparison  is  given  in 
Table  22. 

The  first  three  models,  L-DENSITY,  ARDC-59,  and  US -62,  were 
derived  from  midlatitude  data  near  the  peak  of  the  previous  solar  cycle. 

Even  so,  they  underestimate  the  average  polar  density  below  200  km  by  about 
30  percent.  (A  similar  situation  occurred  when  the  midlatitude  LOGACS 
data  were  compared  with  these  models;  see  Section  IV.)  Evidently,  either 
the  density  boundary  conditions  at  120  km  or  the  vertical  temperature 
gradient,  or  both,  are  in  error.  The  J-W-B  atmosphere  maintains  constant 
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Table  22.  Comparison  of  Density  Models 


Model 

Ra 

a 

b 

100  cr/R 

L-DENSITY0 

1.  50 

0. 

32 

21.  5 

ARDC-59 

1.62 

0. 

31 

19.2 

US-62 

1.56 

0. 

28 

18.0 

J- W-Bd 

1.21 

0. 

20 

16.2 

Champion0 

Prag  (no  bulge) 

1.11 

1. 

17 

15.  0 

1.00 

0. 

12 

11.7 

Prag  (bulge) 

1. 00 

0. 

10 

9.9 

R  =  U/NJVi /pM 
.  ->  1/2 
b<r  =  j[l/(N  -  1)]2.(1  -  P./PM)  f 


Q 

Essentially  Jacchia  I960. 

^Bruce's  modification  to  Walker's  analytical  Jacchia  model 
(Ref.  6). 

0 

See  Ref.  7.  Fractional  density  changes  and  the  time  delay  at 
150  km  have  been  applied  at  all  altitudes. 

^The  models  proposed  here. 

boundary  conditions,  but  has  a  considerably  smaller  temperature  gradient 
than  the  three  previous  models;  the  Champion  model  modifies  both  the 
boundary  conditions  and,  to  a  lesser  extent,  the  temperature  distribution. 

The  Prag  models  proposed  here  were  adjusted  to  fit  the  measured 
densities;  in  particular,  it  was  required  that  the  average  ratio  of  the 
measured  density  to  the  model  density  be  unity.  As  a  consequence,  the 
average  model  density  is  considerably  higher  over  the  polar  regions  than  the 
earlier  models  had  predicted. 

The  last  column  in  Table  22  gives  the  overall  rms  percentage  error 
of  the  various  models  when  compared  with  the  observations.  The  major  por¬ 
tion  of  the  error  in  all  seven  models  occurred  in  the  last  day  and  a  half  of 
the  flight,  during  the  magnetic  storm.  The  average  overall  percentage  error 
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of  the  first  three  models  was  19.6  percent;  these  models  have  no  provision 
for  including  magnetic  activity.  The  J-W-B  model  uses  the  planetary  mag¬ 
netic  index  (derived  from  midlatitude  observations)  to  modify  its  predicted 
density,  and  the  Champion  model  uses  the  College,  Alaska,  K  index,  which 
is  a  measure  of  auroral  activity.  There  is  a  17  and  a  25  percent  improve¬ 
ment,  respectively,  over  the  no-index  models.  The  Prag  models  use  the 
weighted  average  of  the  Eielson  and  Thule  AY  ranges.  The  Eielson  range  is 
essentially  equivalent  to  the  College  K  index,  and  is  also  a  measure  of  auroral 
activity.  The  Thule  range  is  a  measure  of  polar  cap  activity,  which  is  often 
quite  different  than  that  in  the  auroral  regions.  There  is  a  40  percent  (no 
bulge)  and  a  49  percent  (with  bulge)  improvement  over  the  no-index  models. 

Additionally,  the  Prag  models  introduce  a  considerable  amount  of 
structure  in  the  atmosphere  during  disturbed  times.  (During  reasonably 
quiet  times  the  density  of  the  atmosphere  below  200  km  had  little  or  no 
dependence  on  either  latitude  or  longitude,  except  for  the  very  localized 
bulge.  )  A  "trough"  appears  in  the  region  of  the  auroral  oval,  which 
reduces  the  density  thereby  as  much  as  30  percent  (M  =  150)  at  all  altitudes 
(135  _  Z  5  200  km).  A  possible  explanation  for  this  is  the  explosive  heating 
due  to  a  very  large  influx  of  high  energy  particles  in  this  region.  This  heating 
is  presumably  the  source  of  high-velocity  winds  in  the  polar  regions,  as 
reported  in  Section  VII,  so  some  care  must  be  taken  in  interpreting  the 
trough  as  a  true  density  minimum;  about  20  percent  of  the  apparent 
decrease  in  density  could  be  due  to  the  winds  "pushing"  the  vehicle,  thereby 
reducing  the  measured  deceleration.  Since  the  effect  on  drag  measurements 
is  the  same  whether  there  are  tailwinds  or  a  true  density  decrease,  it  was 
decided  to  include  the  trough  in  the  density  model  and  not  bring  in  the  added 
complication  of  dynamic  mass  motion. 

Figure  70  shows  the  average  density  of  the  polar  atmosphere  between 
120  and  190  km  according  to  the  models  presented  here.  The  dashed  lines 
represent  two  extreme  CIRA  models.  For  low  to  moderate  magnetic  activity 
the  new  models  are  lower  than  CIRA  below  140  km  and  roughly  comparable  at 
higher  altitudes;  during  highly  active  times  CIRA  everywhere  underestimates 
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the  density  above  125  km.  As  with  the  midlatitude  model  (Section  IV),  the 
low-altitude  scale  height  is  much  larger  than  previous  model  atmospheres 
ha  ve  predicted. 

Figure  71  is  a  three-dimensional  projection  of  the  atmosphere  near 
the  bulge.  The  "fractional  density"  is  the  ratio  of  the  model  density  to  the 
average  density;  the  figure  has  been  drawn  for  M  =  60  and  applies  to  all 
altitudes  between  120  and  175  km.  The  density  trough  between  350°  and  30°E 
longitude,  and  again  at  270°,  emphasizes  the  highly  localized  nature  of  the 
bulge  even  for  the  fairly  active  time  under  consideration. 

No  altitude -dependent  time  delay,  such  as  was  found  in  the  mid¬ 
latitude  regime,  is  included  here.  Such  a  delay  is  undoubtedly  present,  but 
the  data  were  not  extensilve  enough  to  adequately  determine  its  nature.  As 
more  information  becomes  available,  it  will  be  included  and  will  probably 
significantly  improve  the  performance  of  these  polar  models. 
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SECTION  VII 


LOGACS  WIND  ANALYSIS 


A.  INTRODUCTION 

The  processing  of  LOGACS  data  (Volume  I,  Section  II)  revealed 
large  anor.  alies  in  drag  acceleration  following  the  occurrence  of  the  large 
magnetic  storm.  Correlated  with  this  phenomenon  was  a  change  in  polarity 
of  the  yaw  control  force  and  consequently  a  corresponding  polarity  change  of 
the  yaw  control  moment.  From  this  response  characteristic  and  the  assump¬ 
tion  it  is  externally  caused,  it  follows  that  during  this  interval  the  yaw  aero¬ 
dynamic  moment  changed  polarity.  This  in  turn  leads  to  the  conclusion  that 
there  occurred  a  large  change  in  the  yaw  angle  of  attack  and  the  correspond¬ 
ing  hypothesis  of  large  horizontal  wind  components  normal  to  the  orbit  plane. 

Wind  calculations  could  not  have  been  performed  had  the  LOGACS 
experiment  not  been  instrumented  for  control  gas  firing  data.  Also,  the 
rotating  accelerometer  mode  provided  data  necessary  to  calibrate  the  control 
force  scale  factor  (see  Section  E),  The  attitude  data,  although  not  explicitly 
used  in  the  wind  calculation,  provided  required  information  for  evaluating  the 
control  system  performance. 

B.  FORCE  AND  MOMENT  ALGORITHMS 

Figure  72  presents  a  force /moment  diagram  of  the  Agena  configura¬ 
tion  and  defines  the  symbols  used  in  the  analysis.  The  force  and  moment 
equations  can  be  written  by  inspection  as 

Force 


mA 

y 


f  ) 

cy' 


(46) 


Moment 


M  =  Id>  =  F  X 
z  z  cy  gv 


F  X  -  F  (Y  -  Y  ) 
ay  cp  ax  eg  cp 


(47) 
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ATTITUDE  REFERENCE 
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L 


ATMOSPHERIC  VELOCITY  NORM 


where  the  individual  terms  on  the  right-hand  side  are  defined  in  Fig.  72. 
When  the  side  acceleration  is  measured,  Eq.  (46)  can  be  used  to  calculate 
the  yaw  angle  of  attack'  as  follows 

Fay  =  CNPQSP 

F  =C,QS-QS=F  /C  (product  of  dynamic  pressure  Q 
ax  A  ax  A  and  reference  area  S) 

Fay  =  <CNP/CA|Fa^  =  Fy  '  Fcy 
from  which  p  is  calculated  as 


P  = 


(48) 


where 


CA/C 


NP 


=  ratio  of  aerodynamic  axial  to  normal  force  slope 
coefficients 


F  =  total  side  force  determined  from  processing  LOGACS 
^  accelerometer  data  in  the  rotating  modes 

F  =  control  force  determined  from  processing  LOGACS 
cy  control  valve  firing  data 

F  =  axial  force  determined  from  the  LOGACS  accelerometer 
ax  ,  . 


When  the  accelerometer  is  in  the  fixed  fore  or  aft  mode,  Eq.  (47) 
can  be  used  to  calculate  the  yaw  angle  of  attack.  In  the  perigee  region  the 
control  system  is  operating  on  the  control  deadband  so  that  the  time  average 


When  it  is  assumed  there  is  no  wind  component  in  the  in-track  coordinate, 
the  yaw  wind  velocity  is  calculated  as  Vj  tan  p  «  V^P. 
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of  angular  acceleration  is  zero  (see  Volume  1,  Section  II).  Therefore,  Eq.  (47) 
relates  a  moment  balance  condition  as  follows 


FX  =  F  X  +  F  ( Y  -Y) 
cy  gv  ay  cp  ax  eg  cp 


When  it  is  assumed  that  Y  is  a  linear  function  of/3,  i.  e,, 

cp 


Y  =  Kp 
cp 


where  K  may  be  zero,  the  above  equation  reduces  to 


F  X  =  (C..-X  -  CaK)QSP  +  C  aQSY 

cy  gv  '  NP  cp  A  A  eg 


We  define  the  coefficient  of  p  as 


CMX  -  C  aK  =  CMflX  =  C 
N  cp  A  Np  epe  mp 


where  X  is  defined  as  the  effective  center  of  pressure  and  C  „  is  the 
epe  r  m/3 

moment  coefficient  about  the  vehicle  center  of  gravity.  Using  this  relation¬ 
ship  and  solving  the  (i  results  in 


_  L  X  -  Y 

Fax  c8 


) 


(49) 


Also,  assuming  vehicle  symmetry,  the  pitch  angle  of  attack  can  be  calculated 
as 


a  = 


(50) 


With  the  assumption  of  symmetry,  Cm«  =  CmQ;  and  C^jp  =  Cnq>  The  sub¬ 
script  symbols  a  and  p  are  used  interchangeably  throughout  the  report. 
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During  the  period  when  Eq,  (48)  can  be  used  to  solve  for  p,  a 
consistency  check  can  be  made  between  the  results  obtained  by  Eqs,  (48)  and 
(49).  If  they  do  not  agree  (as  was  the  case  initially,  see  Section  E),  then 
the  following  relationship  is  used  to  determine  the  source  of  inconsistency 


The  terms  on  the  right-hand  side  of  Eq.  (51)  can  be  adjusted  to  produce  an 
expected  (predicted)  value  of  X  g.  Since  the  acceleration  data  errors  are 
small,  the  terms  most  suspect  are  a  scale  factor  on  the  control  force  level 
and  the  value  of  Y 

eg 

C.  AERODYNAMIC  COEFFICIENTS 

The  data  of  Ref.  3  were  used  to  determine  the  aerodynamic 
coefficients  required  in  this  analysis.  Since  force  coefficients  were  pre¬ 
sented  in  terms  of  wind  axis  (lift  and  drag),  they  were  transformed  to  body 
axis  coefficients  using  the  following  transformation 


The  ratio  of  C^  to  C^  was  calculated  as 
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rhe  effective  center  of  pressure  was  determined  from 


X 

epe 


-  L  AC  /CM)  + 
ref  m  N 


eg 


who  re 


eg 


reference  length  (5  ft) 

moment  coefficient  referenced  to  Agena  Station  247 
distance  from  Station  247  to  the  vehicle  eg  (8.7  ft) 


Since  the  coefficients  depend  on  the  molecular  speed  ratio,  which  in 
turn  depends  on  exospheric  temperature,  T  ,  and  altitude,  and  required 
coefficients  were  first  calculated  as  a  function  of  molecular  speed  ratio  for 
the  angle  of  attack  range  considered  in  Refs,  3  (2°  and  4°).  These  are 
plotted  in  Fig.  73.  The  molecular  speed  ratio  as  a  function  of  altitude  and 
exospheric  temperature  is  presented  in  Fig.  74.  From  the  data  of 
Figs.  73  and  74,  plots  of  predicted  aerodynamic  coefficients  are  presented 
in  Fig.  75  as  a  function  of  altitude  for  angles  of  attack  of  2°  and  4°  and 
exospheric  temperature  of  1000°  and  1500°K.  Since  the  predicted  variation 
of  ±1#  5  about  a  mean  value  of  52,  5  (±3  percent),  a  constant  value 

was  used  for  this  coefficient.^  The  predicted  characteristics  of  Xc^e  were 
not  explicitly  used  in  the  wind  calculations,  but  were  used  as  a  guide  in  eval¬ 
uating  Eq.  (51)  (Section  B). 


Reproduction  of  curve  presented  in  Ref.  3. 

^When  -*■  0,  the  molecular  speed  ratio  goes  to  co  and  the  coefficients 
Ca/CjsJq.  can  be  calculated  by  an  area  rule  which  results  in  a  value  of 
1  rad  =  57.  2958°. 
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Fig.  74.  Molecular  Speed  Ratio  vs  Altitude  and  Temperature 
(Based  on  Walker- Jacchia  1964  atmosphere  model 
as  modified  by  Bruce.  ) 
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D.  LOGACS  SIDE  FORCE  PROCESSING 


The  equation  for  LOGACS  accelerometer  counts  is  (see  Volume  I, 
Section  II) 


N.  = 
1 


*l\f 


ra4t 


A  cos 

X 


+ 


t.-l 
1 


A  sin 

y 


*t) 


dt  +  E  +  B 


(52) 


where 

N.  =  incremental  accelerometer  counts  recorded  each  second 

r  =  turntable  radius  arm 
a 

=  angular  rate  of  turntable 

A^  =  acceleration  along  the  -x  axis  due  to  aerodynamic  drag 

A  =  acceleration  along  the  -y  axis  due  to  aerodynamic  side  force 
y  and  control  valve  firings 

E  =  error  term  due  to  angular  rate  and  acceleration  induced  by 
control  firings  and  aerodynamic  torques;  the  mean  value  of  E 
is  zero  with  a  variance  proportional  to  sin  and  the  number 
of  control  firings 

K  =  instrument  scale  factor 

B  =  instrument  bias 


The  instrument  scale  factor  and  bias  were  determined  from  the 

rotating  mode  data  at  apogee.  Since  the  expected  (average)  value  of  E  is 

zero,  it  is  dropped  from  further  consideration.  Since  the  recording  interval 

is  1  sec,  the  integration  can  be  dropped  assuming  the  variables  are  constant 

•  2 

over  the  interval.  The  term  r  $>  is  a  known  constant  for  a  given  accelerom- 

cl  t 

eter  mode.  Applying  these  considerations,  Eq.  (52)  can  be  rearranged  as 


A  cos  +  A  sin  <J>, 
x  t  y  t 


•2 
r  $> 
a  t 


(53) 


L 
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Since  <t>  and  the  mode  event  times  are  known,  <j>  is  known  except 

v  t 

for  the  initial  condition  errors  due  to  turntable  mode  transients.  Included 
with  the  LOGACS  data  were  event  times  when  the  turntable  passed  the  0°  and 
180°  positions  (identified  as  Switch  A  and  Switch  B  data).  The  complete  data 
set  of  event  times  was  processed  to  determine  the  mean  and  variance  of  the 
time  difference  between  predicted  and  measured  event  times  (fore  and  aft 
positions).  The  angular  error  was  determined  as 


A$t  =  4>tAt 


The  angular  means  obtained  are  tabulated  as  follows: 

Switch  A  Switch  B 

Mode  1  2.99°  3.22° 

Mode  2  3.19°  3.45° 


The  standard  deviation  about  the  mean  was  less  than  0.  1°  and  was  probably 
due  to  telemetering  errors  rather  than  any  deviation  of  event  times.  It  is 
seen  that  the  data  indicate  an  angular  error  between  Switch  A  and  B  of  «0.  2°. 
The  Mode  1  startup  transient  resulted  in  a  3.  1°  phase  error  during  Mode  1 
and  an  additional  0.2°  during  Mode  2,  which  was  caused  by  the  motor  speed 
change  transient.  For  the  analysis  presented  here,  a  constant  3.2°  phase 
error  was  used. 

An  undetermined  phase  angle  error,  would  result  in  an  error 
in  the  estimated  forces  as  follows: 


rA 


F 

ax 

cos  $ 

e 

sin  $ 

e 

F 

ax 

A 

F 

L  axJ 

-sin  $ 

e 

cos  $ 

e 

F 

.  ay. 

For  small  angles  this  reduces  to 


F  =F  +  F  sin  4 
ax  ax  ay  e 


F  =  -F  sin  4>  +  F 

ay  ax  e  ay 


Since  Fax  is  approximately  an  order  of  magnitude  larger  than  F  ,  little 
error  results  in  the  estimated  value  of  F  .  However,  the  estimation  error 
for  F&y  is  not  negligible.  If  the  above  expression  for  F^  is  used  in  Eq.  (48), 
it  is  seen  that  there  is  a  nearly  one-to-one  relationship  between  4>e  and  the 
resulting  p  calculated. 

The  method  used  to  determine  A  and  A  of  Eq.  (53)  was  to  assume 

x  y  M 

a  model  of  each  function  as  depicted  in  Fig.  76.  A^  was  assumed  to  vary 
linearly  between  estimated  values  for  A^  at  the  fore  and  aft  phase  angles. 

Ay  was  assumed  to  be  constant  over  the  same  phase  angle  interval.  This 
results  in  a  linear  model  with  time-varying  coefficients  where  the  right-hand 

side  of  Eq.  (53)  is  the  corrected  data  set.  The  parameters  to  be  estimated 

A  A 

are  the  table  points  of  A^  and  A^.  This  model  was  mechanized  in  Post 
Flight  Reconstruction  Program  (PFRP),  which  is  part  of  Modularized  Vehicle 
Simulation  Program  (MVS),  and  is  described  in  Ref.  35.  The  PFRP  esti¬ 
mation  model  chosen  was  the  least-squares  version. 

Seven  sets  of  Mode  1  and  2  perigee  data  were  chosen  for  processing. 
Three  were  before  and  four  followed  the  magnetic  storm  activity.  The  specific 
revolutions  were:  Rev.  38,  which  included  the  time  interval  going  toward 
perigee;  Rev.  39,  coming  from  perigee;  Rev.  41,  which  bracketed  perigee; 
and  Revs.  43,  53,  54,  and  56,  which  all  included  data  of  time  intervals  when 
anomalies  in  drag  acceleration  and  yaw  control  response  were  observed. 

In  the  course  of  processing  the  data  sets,  it  was  observed  that  the 
estimated  values  of  side  acceleration  had  a  consistent  oscillation,  producing 
a  result  greater  in  magnitude  for  phase  angles  from  -n  to  zero  than  for  phase 
angles  between  zero  andTr.  The  source  of  this  characteristic  has  not  been 
identified,  but  could  be  the  result  of  a  phase  angle  error,  or  equivalently  an 
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A,  cos  <J>f  +  sin  :  8  "  \ 


-r  0  *  2f  3»  4t  5t  6f . 22* 

PHASE  ANGLE  4>f,  rod 

Fig.  76.  LOGACS  Side  Force  Processing 
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instrument  misalignment  error,  and/or  an  accelerometer  sensitivity  to 
cross -axis  acceleration.  The  order  of  magnitude  of  the  oscillation  indicates 
the  error  magnitude  is  about  0.5°  angle  or  1  percent  accelerometer  cross¬ 
axis  sensitivity.  To  compensate  for  the  oscillation  effects,  adjacent  values 
of  Ay  estimates  were  averaged  and  tagged  with  the  midpoint  phase  angle 
argument  (zero  or  it).  To  construct  time  histories,  a  linear  interpolation  of 
t  e  average  Ay  values  between  arguments  was  used. 

The  results  of  the  processings  of  the  seven  revolutions  are  presented 
in  Figs.  77  through  83.  The  time  axis  is  referenced  to  the  Mode  1  start  time 
which  is  listed  also  as  TRGO.  These  can  be  directly  correlated  with  data  pre'- 
sented  in  Section  II,  Volume  I.  The  units  for  the  x  and  y  axis  accelerations 
are  in  terms  of  accelerometer  counts  and  are  normalized  so  that  one  count 
equals  10  g.  To  convert  to  force,  the  acceleration  is  scaled  by  the  vehicle 
weight  (1912  lb),  or  the  counts  should  be  scaled  by  1.  912  X  10‘4.  Included  in 
these  plots  are  the  pitch  and  yaw  control  force  data,  also  in  units  of  accelerom¬ 
eter  counts.  These  are  directly  from  the  data  of  Volume  I,  Section  II,  scaled 
by  0.5  to  reflect  the  calibration  results  discussed  in  Section  E. 

Appendix  C  presents  the  Rev.  56  data  set  used  to  produce  the  results 
presented  in  Fig.  83.  This  is  typical  of  all  the  data  sets.  Figure  C-l  pre¬ 
sents  the  corrected  raw  data  [right-hand  side  of  Eq.  (53)]  and  Fig  C-2 
presents  the  residuals  of  the  fit  (right-hand  side  minus  the  left-hand  side  after 

i mating  Ax  and  Ay).  The  large  noise  level  is  primarily  caused  by  the 
effects  of  control  valve  firings. 

E-  PARAMETER  CALIBRATION 

The  initial  attempt  to  process  the  force  data  of  Section  D  to  arrive  at 
angle- of- attack  histories  was  to  use  Eq.  (48)  to  calculate  (3,  use  Eq.  (51)  to 
calculate  X^,  and  use  in  Eq.  (50)  to  calculate  a.  This  led  to  extremely 

erratic  results  for  X^  [some  calculated  values  were  ahead  and  seme  were 
behind  the  vehicle,  dependent  upon  polarity  of  the  denominator  of  Eq.  (51)] 
and  corresponding  erroneous  calculations  of  a. 
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7-15 


A 


liSiyimiisMiiH 


o 


lijHiipjlijjiliiiiljiiiiililllljjjjjiljl 

liliiSimiifl 

-**•* ^  *•♦•1  uni  iitm ism  iip# i!iu  •  ixi « 1 1 *x.  ««<■  ■  . •  *  iirn iwi  ■  t « *>■  ■  .  * .  t « 


i 

i 


II  fiiiaifi  sklBHi  ■■  BUilttii  mm  iiiiiiiii:  .4tti  ilia  liiiil 

m$a  ijiiissiasiffi  mm  | 

•*”i 3.1. ^  Ltii iji*  lin  jin uui n*H  mans.. . 


iilii 
SSHHliil 

p^wnm 


Shmh 

iMl» 


mm 


lijiliilliiiifiilKliliiHsi 

iSHiiiiEliiii! 


iim  fiHUim  mi  hhj  iHjiiHx  yHiHUj  ih>-i 

|nff:  IfKlIrf* ,  ftvtinftl  rimiixn  tnvmlji  flilllint :r!!:| 

IKtfi  ftffifffB  Bir  :4fe  tHflHIB  BHHM  HHIHHi  HfflIHH  HHI IfHI H  r  ;;l 

I rHtl  rrffllttfllfflii.  tl  fSnHfr  (ftnrttH llnlllttl lirtl Ifni llllliitllffxljlrn? !lt:!l 


i  jif it  iitti  rni;  ikzxi  mu  nisi  uni  mill 

I  HR!  mm  Iffli  rnti  Ixlf  I  nfll  tfifll 


WMMM\ 


«itK 

mmm& 


iSiilliill  yiiil 

ft  HI  lift!  fltlil 

mm 

■:4 


TIME,  sec 


Fig.  78.  Revolution  39  (TRGO  =  208640) 


7-16 


*! 

c[ 

s| 

• 

! 

*1 

l 

4! 

1 

i 

i 

■ 

• 

> 

4 

:::::: 

ilill 

ii 

sissiti: 

iiii 

!llll: 

linillli; 

22222 

!8 

t.  0 

1 

1 

Bias 

■* 

s: _ ;::: 

SB 

s;  !* 

:::::::: 

\ww 

..11; 22121 
tHHHx 

m 

H 

m 

pj)l; 

saas 

i!:l:l!i 

il;; 

ii:: 

;Sjs 

ftttimt: 
1222222*2; 
W  sir 

2222222X1. 

1 

iii;; 

BBS 

TjTTT;". 

l 

m 

iTZr. 

txttttt: 

iijjiii 

11 

llll: 

111 

3  0 

HK 

1-0 

1 

isSiHii 

as 

llll : 

J* 

.31. , 

1B8 

m 

|j| 

Ililll 

iii  iiii! 

aas 

jillffi! 

HEnm 

is 

mu 

111 

'I. 7. 37. *7*2* *: 

Hill 

•Hi: 

7.  it. 

m 

:I7r:!t27, 

::::: 

:;;;; 

1 

1 

ilnlitlli!; 

!!!”  SSill  11 

**..112. 

aihi:: 

nil;.:: 

ii 

X:;. 

ijlSni: 

HUM 

.12 

3jjj 

ili 

:::: 

ji 

11  11  Ill 

nit. 

Hliijj:, 

tl?. i;tr 

IHililii 

i 

TTtttfTTfl 

— 

ta 

23X11 LZ. 

gfnir 

.on::. 

ii! :::::: iiii ii 

3. .. 

iliii 

m 

.X..2 ....] 

11:11 

;27.: 

iitij 

::::: 

am 

11 

ajni:::: 

iBhIRm 

ma 

dfillfl 

iiflllPl 

Hill 

*Cv*t : i •*: 

:s 

.311777... 

aa 

;;;;; 

* 

1*33: ■**’•  •’ 

:::xxxx: 

11 

HiHg. 

III! 

as 

ill! 

»a 

227X212222 

.173. .721. 

*j«2|p||j 

III 

l/r 

1 

aas 

sa 

aa 

P 

f 

EFJ:| 

..a:  u: 

Sag 

ill 

Iiii 

a 

Hi 

1;  1; 

> 

to 

2272222222  n 

..... ..n.  iz 

aa 

2272  22211 

as 

iHH 

txj.rrii 

ki;: 

llll 

m 

ilil 

jlSjjgi 

II 

ifflsn 

,x. iz 

:T t  # 

ICtaaaC 

iiii 

:S3J 

i;;;j :::: 

;;;;! 

s:ss: 

r:r* 

S:27'*o 

«*•  .  i* 

••  •*  •• 

il 

1 

1 

llll: 

llll 

111 

11 

if 

ia 

2232 

d!i 

a 

I!'J* 

A 

Itiii 

ill 

ili 

iili 

il 

m 

IB 

iiii 

\li\£ 

a 

m 

lift 

a 

U 

1 

jjjj; 

3o: 

m 

m 

ill 

ii! 

Li 

ii 

il'illl! 

as 

11 

a 

ill 

ili 

a 

I; 

MA 

a 

1 

as 

m 

IPiH 

m 

m 

a 

a 

m 

1 

a 

i 

Kill 

ill 

[UU 

Hill 

film 

nitrn 

iii 

m 

11 

itiii 

il 

Iilii 

a 

a 

m 

a 

li 

itHi 

:ttS 

ita 

itiH 

iji 

:itc 

m 

iili 

a 

m 

iilii 

::::: 

iliii 

s:m 

as 

li  !:rr 

i 

i 

i 

Hiii 

sa 

m 

IH 

a 

a 

il 

iliii 

;;:;i :) 

JHj 

tiii 

tlffi 

ran  it  | 

B  itHi 

: it 

tJUUii 

*r*r 

ijia 

1H:: 

rfeti  ?! 

I!  lift 

Ej'Hi 

mil 

ftwltrl 

fill;!; 

HOirt 

-.IS 

iff*! 

tlrt* 

lift! 

m 

H 

m 

ai 

Hi! 

m 

H 

i 

Hi! 

urn 

»ii 

itn 

QU] 
nil  j 

Ufl! 

m 

nm 

as  ii 

HI 

iili 

an 

as 

ia 

m 

m 

li 

ii 

11 

•;Hi 

ITfil 

a 

Hil! 

!M3 

1! 

iilii 

;;;;; 

..... 

WV 

?  600 

I 

|  400 
( 

8  200 

0 

F 

MBS 

fill 

1 

1! 

tin 

uui 

rttli 

as 

ill 

m 

la 

a 

iiii! 

'iiii 

11 

mm 

until 

1 

mi 

wi 

aa 

ifui 

liHilK 

tffiml 

ill 

ai 

ia 

a 

a 

ii 

iili; 

li;:1 

iiii! 

***** 

*••*« 

:::« 

a 

il 

m 

lliiii 

iilHii 

tiHt:* 

as 

liliil 

■ 

nia 

niiii 

aa 

« 

1 

li 

i 

F 

il 

i 
li 
1 
«* 

ii 

[Mil 

$$ 
ft  hi 

i; 

m 

tun 

jffii 

im! 

m 

m 

lil 

Bill 

£2? 

iili 

iiii 

Iii! 

ill! 

M 

as 

as 

sa 

ii 

sa 

Sil! 

11 

iiii 

llll 

iiii 

a 

a 

iii 

lliiii 

11 

h 

a 

II;;; 

iilii 

a 

KS 

ilil! 

a 

ii 

iliii 

Iilii 

|i'|: 

iili; 

iilii 

iliii 

mil 

it!!] 

!tm 

illg 

111!!? 

ill':) 

ill! 

sa 

ii 

ill 

i 

i 

il 

ii 

il 

Hill 

Iiii 

iiii 

Hi 

iii 

m 

ii 

a 

Iiii 

iili! 

JS 

sin 

'iiii 

’ 1 

ill;; 

P! 

ih’li 

’ill!!! 

111! 

1! 

111! 

Ei 

iliii 

m 

as 

im 

iiii 

a 

iiii 

iilii 

as 

sa 

11 

ill! 

1 

iilii 

Iiii 

as 

sa 

Iiii 

ill 

mi: 

Iilii 

il'- 

a 

11  if 

llll 

li 

i 

ii 

ill 

111 

as 

na 

;!ii 

iliii 

iii; 

iili! 

mil 

4 

M 

aa 

iin 

li 

1 

tun 

1 

ill 

Hill 

Hi 

Ilil 

iiiii 

ii 

iiii 

iilii 

iilii 

11 

llll! If 

mini 

ii 

ill! 

i  in 

dt 

nil! 

11 

Hill 

as 

iliii 

Iili 

mil 

iiii 

ii 

iii;; 

iitij 

ii 

llffli 

m\\ 

111! 

llll 

i 

li 

il 

iii 

III! 

Ili! 

tiffi 

m 

ii 

ill 

Iilii 

| 

lil 
El  III 

m 

in 

11 

Iiii 

i 

a 

iili! 

iilii 

a 

■■ 

iili 

iii) 

iii 

m 

ill!! 

il 

jjs;: 

\m 

r 

! 

L 


7-17 


7-18 


TIME,  sec 


Fig.  81.  Revolution  53  (TRGO  =  282240) 


TIME,  sec 


Fig.  82.  Revolution  54  (TRGO  =  288384) 
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Fig.  83.  Revolution  56  (TRGO  =  298624) 
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Tho  next  approach  was  to  use  several  different  scale  factors  on  the 

control  force,  F  ,  and  several  different  values  of  Y  in  Eq.  (51)  to  observe 
cy’  eg  M 

the  resulting  time  histories  of  X  ,  a,  and  3  for  the  data  sets  of  each  of  the 

epe 

seven  revolutions.  Tho  best  results  were  obtained  when  a  scale  factor  of  0.5 

was  used  on  control  force  and  the  predicted  value  of  Y  was  used.  Except 

for  Rev.  53,  the  resulting  time  history  of  X  was  near  constant  and  did  not 

’  epe 

reflect  the  altitude  dependence  predicted  by  the  curves  of  Fig.  75.  There 

does  appear  to  be  a  possible  temperature  (T  )  sensitivity  in  that  the  results 

for  Revs.  38,  39,  41,  and  43  imply  a  constant  value  of  2  ft  for  X  e,  while 

Revs.  54  and  56  appear  to  show  a  somewhat  lower  value  for  X  .  The  above 

epe 

numbers  are  approximate,  based  on  an  eyeball  estimate.  The  principal 


characteristic  from  this  evaluation  is  that  to  first  order  X 


is  a  constant 


epe 

and  does  not  depend  on  altitude  for  the  range  of  altitude  (75  to  110  nmi)  con¬ 
sidered  here.  This  result  implies  that  either  the  moment  coefficients  and 
resulting  aerodynamic  center  of  pressure  as  a  function  of  molecular  speed 
ratio  are  in  error  (Fig.  73),  or  the  molecular  speed  ratio  versus  altitude 
(Fig.  74)  is  in  error  for  the  altitude  regions  considered. 

With  a  constant  value  for  X  ,  the  calculation  of  angles  of  attack 
using  the  moment  relationship  [Eqs.  (49)  and  (50)]  is  in  error  by  a  possible 
scale  factor,  depending  on  the  value  of  ^CpC  chosen.  To  further  investigate 
the  proper  choice  of  numerical  values  to  use,  we  made  a  consistency  check 
by  comparing  the  results  obtained  for  (3  from  Eq.  (48)  to  those  obtained  from 
Eq.  (49).  For  this  test,  a  value  of  2  ft  was  used  for  X  As  before,  the 

value  of  0.5  for  control  force  scaling  and  the  predicted  value  for  Y  pro¬ 
duced  the  most  consistent  results.  The  comparison  of  Rev.  54  and  Rev.  56 
also  showed  that  a  slightly  smaller  value  for  X  ^  would  have  produced  a 
closer  comparison;  i.  e.  ,  for  these  revolutions  the  force  equation  produced  a 
slightly  larger  calculated  (3. 
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Some  characteristics  of  equation  sensitivities  to  data  errors  can  be 
explained  by  the  following  comments.  The  sensitivity  of  the  calculated  angle 
of  attack  to  lateral  eg  is  determined  from  Eq.  (49)  as 


dP 


c)Y 


=  -26° /ft  -  -2.  16° /in. 


eg 


mo 


An  error  in  predicted  eg  thus  acts  as  a  bias  of  calculated  angle  of  attack. 

Comparing  Eqs.  (48)  and  (49)  with  respect  to  control  force  scaling,  one  sees 

that  the  control  force  is  scaled  by  a  factor  of  X  /X  in  Eq.  (49)  N4.  5). 

1  gv  epe  M 

Also,  the  polarity  of  the  control  force  term  is  different  in  each  equation. 

Therefore,  if  a  control  scale  factor  error  exists,  one  equation  will  give  a 

greater  value  and  the  other  a  smaller  value  than  the  true  value  in  the  absence 

of  other  errors.  Additionally,  the  control  force  is  generally  a  factor  of  2  or 

more  smaller  than  total  side  force  so  that  Eq.  (48)  is  relatively  insensitive 

to  a  scale  factor  error  of  the  control  force.  The  coefficient  C  .  /CXT  is 

A  No 

common  to  both  expressions. 

With  respect  to  the  control  force  scale  factor,  it  had  been  determined 
(Volume  I,  Section  II)  that  the  control  firings  were  a  pulse  of  minimum  dura¬ 
tion  (20  msec).  Also,  the  nominal  control  force  level  was  assumed  to  be 
0.5  lb.  Therefore,  the  processing  for  this  earlier  work  was  based  on  an 
impulse  of  0.01  lb-sec  per  gas  valve  firing.  The  pulse  referred  to  here  was 
the  minimum  electrical  current  pulse.  Since  the  valve  does  not  respond 
instantaneously  to  the  current  pulse,  the  full  impulse  is  not  realized.  Test 
results  at  LMSC  had  indicated  that  only  63  percent  of  the  ideal  impulse  is 
achieved  when  valves  are  pulsed  at  the  minimum  duration  for  force  levels  of 
10  lb.  No  tests  have  been  run  for  the  0.5 -lb  force  level.  With  these  consid¬ 
erations,  a  0.  5  scale  factor  is  within  expected  tolerances  of  the  control 
system. 


t 
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It  is  of  interest  to  compare  the  time  histories  of  a  and  p  derived 
from  the  data  with  those  predicted  from  the  orbit  characteristics.  The  pitch 
angle  of  attack  o  is  related  to  the  flight  path  angle  and  attitude  as 

a  -  G  -  V  (54) 


where 

0  =  pitch  attitude  determined  from  horizon  sensor  data 

V  =  flight  path  angle  determined  from  the  Best-Fit  Ephemeris  (BFE) 

Therefore,  when  the  vehicle  is  going  toward  perigee,  a  >  0;  when  coming  from 
perigee,  a  <  0.  For  the  yaw  angle  of  attack  (3,  the  predicted  model  chosen  is 
one  which  assumes  the  atmosphere  rotates  with  the  earth.  For  a  polar  orbit 
and  the  yaw  attitude  reference  of  Agena,  the  following  relationship  for  pre¬ 
dicted  p  results 


P  =  -(v  +  Rwe  cos  \  [sgn  \]) 


(55) 


where 


+  = 

R  = 

u>  = 
e 

\  = 
[sgn  \]  = 


yaw  attitude  with  respect  to  the  orbit  plane  (sum  of 
+  4>0,  see  Fig.  72) 

radius  from  center  of  earth 

rotation  rate  of  the  earth 

latitude 

sign  of  latitude  rate 


Therefore,  when  the  vehicle  is  going  south,  p  >  ;  when  going  north,  p  <  ^ . 

The  data  processing  results  of  the  seven  revolutions  are  shown  in 
Figs.  84  and  85.  Presented  are  the  predicted  and  processed  data  points, 
where  the  predicted  curves  are  smooth  lines  and  the  processed  data  are 
identified  with  a  symbol.  For  p,  both  the  results  from  the  force  equation 
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ALTITUDE,  nmi  beta,  d«, 

FORCE 


L 


Btt 

SSSS 


[Eq.  (48)]  and  the  moment  equation  fEq.  (49)]  are  plotted.  The  high  degree 
of  correlation  between  predicted  and  processed  data  for  Revs.  38,  39,  and  41 
is  apparent.  Also  apparent  is  the  good  correlation  of  results  (except  for  part 
of  Rev.  53)  obtained  by  the  two  separate  p  algorithms  on  all  seven  revolutions. 
Finally,  it  is  observed  that  the  differences  of  the  p  angles  from  those  pre¬ 
dicted,  which  are  of  magnitudes  greater  than  7°  and  of  both  polarities,  are 
contained  in  the  data,  and  that  the  transitions  from  one  polarity  to  the  other 
occur  in  very  short  time  intervals. 

F.  MISSION  PITCH  AND  YAW  ANGLE-OF- ATTACK  TIME 
HISTORIES 

Equations  (49)  and  (50)  require  time  history  data  of  control  force  and 
axial  force  only  to  calculate  angles  of  attack.  Therefore,  these  equations  can 
also  be  used  for  the  modes  when  the  accelerometer  orientation  was  fixed, 
provided  the  control  system  was  on  its  deadband  (Iu>  =  0).  Since  the  time 
history  data  were  available  from  Volume  I,  Section  II,  no  further  data  reduc¬ 
tion  was  necessary.  To  provide  consistent  time  arguments  on  the  control 
force  and  the  axial  force  data,  we  filtered  the  latter  to  produce  data  at  the 
100-sec  data  rate  of  the  control  data  processing.  In  general,  the  attitude  was 
on  the  control  deadband  when  the  drag  acceleration  level  reached  a  100-count 
level,  Therefore,  for  this  analysis  all  acceleration  data  below  the  100-count 
level  were  deleted  from  further  consideration  (in  general  this  limits  the 
analysis  to  altitudes  below  110  nmi).  When  the  data  duty  cycle  is  low,  the 
confidence  level  of  the  processed  control  force  data  is  also  low  (see  Volume  I, 
Section  II).  Therefore,  control  force  values  which  were  processed  for  duty 
cycles  less  than  90  percent  were  deleted  from  consideration  here.  With  the 
above  data  editing,  the  remaining  data  were  processed  for  angle -of-atlack 
histories  using  Eqs,  (48)  and  (49).  In  this  processing,  a  control  force  scale 
factor  of  0.5,  a  X  of  2  ft,  and  the  predicted  eg  values  were  used. 


Duty  cycle  is  defined  as  the  ratio  of  the  number  of  raw  data  points  contained 
in  a  filter  span  to  the  number  of  points  of  a  filter  span. 
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The*  results  of  this  processing  are  plotted  and  presented  in  Fig.  C-3 
of  Appendix  C.  Included  in  this  plot  are  the  predicted  angle -of-attack 
histories,  th  ■  processed  angle -of-attack  histories,  the  100-sec  filtered  axial 
acceleration,  and  the  best-fit  ephemeris  values  of  altitude  and  axial 
acceleration. 

It  can  be  observed  from  these  plots  that  there  exists  good  correlation 

of  predicted  and  processed  angle -of -attack  histories  before  the  magnetic 

storm.  There  are  some  observed  points  in  the  early  portion  of  each  rev. 

(going  toward  perigee'  that  seem  to  be  erroneous.  Two  problems  exist 

for  this  region.  The  first  is  that  the  control  system  may  not  have  settled  on 

the  control  deadband;  therefore,  F  is  zero  until  the  attitude  reaches  the 

cy 

deadband,  and  then  F  is  large  when  it  first  reaches  the  deadband  to  stop  the 
vehicle  angular  rate.  The  second  problem  is  associated  with  force  levels. 
When  the  average  forces  are  small,  the  induced  noise  on  accelerometer  data 
due  to  control  firings  is  greater  in  proportion;  consequently,  less  accuracy 
can  be  expected  for  these  calculated  angles  of  attack.  Nevertheless,  good 
correlation  exists  before  the  storm,  and  large  yaw  angles  of  attack  are 
observed  after  the  storm.  The  pitch  angle  of  attack  (sensitive  to  vertical 
winds)  seems  relatively  unaffected  by  the  magnetic  storm  activity. 

For  these  plots  the  value  of  0[Eq.  (54)]  used  was  1°  and  the  value 
of  vfEq.  (55)j  used  as  0°.  This  implies  that  there  existed  a  yaw  reference 
error  of  1°  or  a  V  error  of  0.46  in.  Either  are  within  expected  tolerances 

eg  r 

and  cannot  be  separated  by  this  analysis.  Also,  control  force  scaling  and  the 

coefficient  C./C  can  influence  the  results  to  some  extent. 

A  me 

G.  YAW  WINDS  IN  GEOGRAPHICAL  COORDINATION 

In  order  to  determine  the  source  of  the  large  winds,  we  decided  to 
first  display  the  wind  components  in  geographical  coordinates.  Figure  86 
presc  nts  a  plan  view  of  the  general  orbit  characteristics  for  this  mission  (the 
orbit  inclination  was  essentially  90°).  Perigee  occurred  in  the  general  latitud 
region  of  45°  N.  The  altitude  below  110  nmi  was  generally  contained  in  the 
sector  from  the  North  Pole  to  the  equator.  It  therefore  was  convenient  to 
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present  the  orbit  trajectory  (geographical  coordinates)  looking  down  from 
the  North  Pole  in  terms  of  colatitude  and  longitude.  For  each  calculated  wind 
component,  a  line  normal  to  the  trajectory  was  constructed  whose  magnitude 
was  proportional  to  the  difference  between  calculated  and  predicted  p,  and  the 
direction  is  based  on  the  polarity  of  the  difference.  An  arrow  was  placed  at 
the  end  of  the  line  to  indicate  the  wind  direction. 

Figure  87  presents  the  results  from  the  parameter  calibration  of 
Revs.  38,  39,  41,  and  43  (Fig.  84  time  histories)  in  the  described  geographi¬ 
cal  coordinates.  Similarly,  Fig.  88  presents  the  results  of  Revs.  53,  54, 
and  56  (Fig.  85  time  histories).  The  rev.  number  associated  with  each 
is  also  indicated  on  the  figures;  an  arrow  displays  the  trajectory  direction. 

The  scale  used  for  p  is  indicated,  where  one  major  division  is  equal  to  10°  p 
(=45  00  fps  wind  magnitude). 

Similarly,  the  time  history  of  Fig.  C-3  is  presented  in  geographical 
coordinates  in  Figs.  89  through  92.  v  To  avoid  trajectory  overlap,  the  data 
are  presented  by  mission  day  (from  23  May  to  26  May).  It  is  readily  observed 
that  very  little  deviation  from  predicted  wind  components  existed  prior  to  the 
magnetic  storm  activity  (Rev,  5  through  Rev.  41).  Following  the  storm, 
large  wind  magnitudes  of  both  polarities  are  evident  in  the  northern  latitudes, 
particularly  for  the  data  of  26  May. 

The  26  May  data  (Fig.  92)  are  of  greatest  interest.  The  overlap 
observed  is  due  to  the  orbit  altitude  decay  with  time,  so  that  the  110-nmi 
altitude  occurred  before  reaching  the  North  Pole.  Some  of  this  overlap  is 
also  due  to  the  northerly  rotation  of  perigee  during  the  mission.  Figures  93 
and  94  amplify  the  data  presented  in  Fig.  92.  The  colatitude  was  reduced  to 
45°  and  the  /?  plot  scale  was  extended  so  that  two  major  divisions  equals  10°. 
Figure  93  plots  the  data  for  the  trajectory  going  toward  the  pole  while  Fig.  94 
plots  the  data  from  the  pole  to  a  colatitude  of  45°.  Also  included  are  lines  of 


The  method  of  plotting  these  data  was  to  plot  the  first  observed  wind  com¬ 
ponent  for  a  particular  revolution  and  then  to  continue  the  trajectory  plot  to 
a  colatitude  of  95°.  If  no  points  were  observed,  then  the  corresponding 
trajectories  were  not  plotted. 
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Fig.  87.  Wind  Components  in  Geographical  Coordinat 


(EAST) 


-3 


ig.  93.  Expanded  View  of  Fig.  92  (Trajectory 
Toward  Pole) 


TRAJECTORY  FROM  THE  POLE 
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constant  altitude  of  110,  100,  95,  and  90  nmi  in  Fig.  93,  and  lines  of  85,  81, 
and  77  nmi  in  Fig.  94. 

H.  WIND  ANALYSIS  CONCLUSIONS 

The  results  clearly  confirm  the  theory  that  the  atmosphere  normally 
rotates  with  the  earth  up  to  and  including  110  nmi  altitude.  However,  the 
vehicle  forces  produced  by  responding  to  this  wind  component  are  highly 
dependent  on  the  location  of  the  vehicle  center  of  gravity  and  the  type  of 
vehicle  control  system.  In  low  altitude  regions  (below  80  nmi),  the  magni¬ 
tude  of  the  lateral  forces  generated  can  exceed  the  magnitude  of  geopotential 
harmonic  terms.  With  the  configuration  of  this  experiment,  lateral  accel¬ 
erations  of  5  to  15  Kg  were  experienced  due  to  aerodynamic  and  control  forces. 

Following  the  magnetic  storm,  the  results  of  the  data  analysis  con¬ 
firm  the  hypothesis  of  large  velocity  wind  components  normal  to  the  orbit 
plane  (apparently  restricted  to  the  polar  regions).  It  therefore  can  be  con¬ 
cluded  that  wind  components  exist  in  the  in-track  coordinate  as  well.  Con¬ 
sequently,  the  analysis  of  drag  acceleration,  and  its  variation,  must  include 
the  effect  of  wind  as  well  as  the  effect  of  density  variation.  For  example,  a 
5000-fps  wind  velocity  in  the  in-track  coordinate  could  result  in  a  40-percent 
variation  in  drag  acceleration. 

The  aerodynamic  coefficients  are  dependent  on  the  predicted 
molecular  speed  ratio,  which  in  turn  is  dependent  on  the  model  of  temperature 
versus  altitude  (Fig.  74).  In  the  processing  of  the  moment  equations,  as 
noted  in  Section  F,  the  variation  with  altitude  of  the  aerodynamic  moment 
coefficient  was  not  demonstrated.  This  suggests  that  the  predicted  molecular 
speed  ratio  versus  altitude  is  in  error.  If  the  predicted  molecular  speed 
ratio  is  in  error,  then  the  predicted  aerodynamic  drag  coefficient  versus 
altitude  is  also  in  error.  This  in  turn  affects  the  accuracy  of  the  derived 
densities  from  drag  acceleration  data.  Additional  analysis  will  be  required 
to  support  this  conclusion. 
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SECTION  VIII 


THERMOSPHERIC  CONVECTIVE  INSTABILITY 

A.  INTRODUCTION 

It  has  been  known  for  some  time  that  local  winds  up  to  100  m/sec 
exist  in  the  thermosphere  even  during  magnetically  quiet  times.  Some  of 
these  disturbances  are  undoubtedly  related  to  the  traveling  gravity-heat 
conduction  type  of  wave  phenomena  (Refs.  36  and  37).  Recently,  vapor  trail 
measurements  at  low  latitudes  during  geomagnetically  active  times  reveal 
strong  wind  components  of  approximately  200  m/sec  at  180-km  altitude 
(Ref.  38).  The  direction  and  speed  of  these  winds  seem  to  agree  quite  well 
with  the  hypothesis  that  they  are  generated  by  the  existence  of  polar  high 
pressure  regions  during  geomagnetic  activity  (Ref.  39).  Despite  the  lack  of 
direct  proof  of  the  existence  of  such  polar  high  pressure  regions  in  the 
thermosphere,  it  has  been  conjectured  that  they  are  generated  by  particle 
and/or  magnetic  energy  absorption  below  160  km  (Refs.  40  and  41),  Monte 
Carlo  calculations  have  shown  that  particle  energy  deposition  occurs  essen¬ 
tially  in  a  limited  layer  of  the  thermosphere  (Ref.  25).  Although  the  low 
latitude  winds  are  good  indirect  evidence  of  thermospheric  heating  in  the  polar 
regions  during  magnetic  storms,  direct  evidence  is  as  yet  not  available. 
Theoretically,  if  such  heating  exists  and  is  strong  enough  a  direct  consequence 
would  be  the  production  of  strong  atmospheric  motion  since,  under  suitable 
conditions,  a  fluid  layer  heated  from  below  is  convectively  unstable. 

B.  OBSERVATION 

Polar  anomalous  winds  were  observed  on  26  May  1967  by  the  1967- 
SOB  satellite.  The  measurements  were  recorded  when  the  satellite  passed, 
in  near  circular  orbit,  within  approximately  30°  of  the  North  Pole.  Since 
the  satellite  has  a  perigee  of  145  km  at  approximately  40°N  latitude,  the 
observations  were  recorded  in  a  thermospheric  layer  of  approximately  60  km, 
i.e,,  from  220  to  160  km.  This  particular  satellite  ephemeris  is  rather 
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important  because  it  establishes  the  polar  character  of  the  anomalous  winds, 
for  if  it  were  global  in  character  they  would  have  been  observed  in  the  low 
latitude  regions  where  the  satellite  passed  through  a  layer  of  similar 
altitude. 

Approximately  8  hr  after  the  peak  of  the  magnetic  storm  at  the  polar 
regions,  wind  patterns  were  seen  to  persist  for  15  hr  with  very  little  change. 
The  horizontal  components  of  wind  velocity  exhibit  well-organized  direction 
changes  with  a  characteristic  length  of  2000  to  4000  km.  Since  the  horizontal 
component  changes  very  smoothly  within  one  horizontal  characteristic  length, 
the  altitude  change  of  the  satellite  within  that  horizontal  length  allows  us  to 
determine  the  minimal  vertical  characteristic  length,  which  turns  out  to  be 
1  5  to  30  km.  The  maximal  horizontal  speed  observed  was  approximately 
1  km/sec,  while  the  maximal  vertical  speed  was  100  ±  100  m/sec.  There 
seems  to  be  an  increase  of  maximal  horizontal  speed  with  altitude  of  less  than 
one  e-fold.  The  full  wind  pattern  unfortunately,  cannot  be  mapped  because  of 
the  lack  of  information  on  the  in-track  wind  component.  Finally,  the  growth 
phase  of  the  wind  pattern  seems  to  last  one  or  two  satellite  orbits. 

From  the  physical  point  of  view,  there  are  several  significant 
features  which  must  be  considered  if  one  attempts  to  understand  the  meaning 
of  the  measurements. 

The  anomalously  high  horizontal  wind  speeds  (~1  km/sec)  at  a 
typical  altitude  of  200  km,  where  the  sound  speed  is  approximately  0.8  km/sec, 
would  exclude  the  interpretation  of  such  wind  phenomenon  as  gravity  waves, 
which  are  essentially  stable  small  amplitude  oscillations  akin  to  sound. 
Although  large  wavelength  gravity  waves  may  conceivably  have  a  phase  velocity 
of  order  of  1  km/sec,  it  must  be  emphasized  that  the  observation  is  on  the 
molecular  flow  speed  and  not  the  phase  velocity  of  the  disturbance.  Hence,  it 
is  eminently  clear  that  we  cannot  explain  the  observation  on  the  basis  of  small 
amplitude  stable  oscillations.  On  the  other  hand,  the  necessity  of  small 
amplitude  disturbances  (~50  m/sec)  to  grow  to  anomalously  large  speeds  finds 
a  rather  natural  explanation  if  we  are  to  interpret  the  observations  in  terms 
of  some  hydrodynamic  instability  of  the  thermosphere  that  allows  such  growth 
to  occur. 
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The  close  association  between  the  wind  observations  and  the 
magnetic  storm  as  well  as  evidences  of  atmospheric  heating  during  magnetic 
disturbances  point  rather  directly  at  the  cause  of  such  a  thermospheric 
instability,  namely,  excessive  heating  of  the  thermosphere  at  a  particular 
level.  If  one  adopts  such  a  hypothesis,  the  polar  character  of  the  anomalous 
winds  immediately  finds  a  rather  natural  explanation  since  the  polar  caps  are 
directly  accessible  to  energetic  particles,  which  can  deposit  their  energy  in 
the  thermosphere. 

C.  CONVECTIVE  INSTABILITY 

The  observational  evidence  of  the  previous  sections  strongly  suggests 
that  the  anomalous  winds  can  be  interpreted  as  convective  motion  due  to  sta¬ 
tionary  heating.  Since  it  is  likely  that  ene  'gy  is  deposited  in  the  thermosphere 
below  160  km  during  a  magnetic  storm  (Refs.  40  and  41),  we  envisage  a  layer 
of  atmosphere  heated  from  below  160  km  as  the  basis  of  our  model.  Con¬ 
vective  instabilities  in  a  fluid  heated  from  below  have  been  studied  extensively 
in  the  astrophysical  and  meteorological  regime.  Detailed  studies  of  convective 
instability  have  been  made  by  Lord  Rayleigh  (Ref.  42),  Pellew  and  Southwell 
(Ref.  43),  and  Chandrasekhar  (Ref.  44)  for  incompres  sible  fluids.  More 
recently,  convective  instability  in  a  polytropic  gas  has  been  studied  by  a 
number  of  authors  [see,  for  example.  Ref.  45)].  Most  of  these  are  considered 
for  conditions  appropriate  to  the  granulation  phenomenon  in  the  solar  photo¬ 
sphere,  where  convective  heating  is  a  dominant  feature.  In  the  earth's 
thermosphere,  where  the  density  is  essentially  exponential,  the  conditions 
are  somewhat  different.  Unlike  the  solar  photosphere,  the  earth's  thermo¬ 
sphere  is  only  occasionally  heated  strongly  from  below,  if  at  all;  hence,  we 
expect  unstable  situations  to  occur  only  rarely  in  the  thermosphere.  Further¬ 
more,  it  is  not  clear  that  a  temperature  inversion  can  be  set  up  in  the 
thermosphere  to  generate  convective  instability. 

To  study  the  generation  of  an  unstable  temperature  inversion,  we 
first  recall  that  steady  heat  conduction  with  only  vertical  motion  is  clearly 
possible  in  the  thermosphere.  Indeed,  one -dimensional  models  with  very 
slow  vertical  motion  have  been  considered  (Ref.  46)  in  order  to  account  for 
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observed  time  delays  between  atmospheric  heating  and  the  density  response 

of  the  thermosphere.  For  moderate  energy  deposition,  these  models  show 

no  tendency  to  form  unstable  inverse  temperature  gradients.  The  reason 

for  this  occurrence  is  rather  clear  because  the  vertical  motion  due  to  the 

expansion  of  the  heated  air  quickly  spreads  the  heat,  so  that  after  a  few  hours 

only  a  general  rise  of  the  temperature  profile  is  evident.  On  the  other  hand, 

if  a  larger  amount  of  energy  is  deposited  over  a  long  period,  say  8  hr,  the 

situation  becomes  quite  different.  The  temperature  increase  near  the  level 

of  steady  heating  now  overtakes  the  tendency  to  spread  the  temperature  by 

mass  motion,  and  a  local  temperature  inversion  in  the  temperature  profile 

may  result  in  the  model.  Figure  95  shows  the  temperature  profile  calculated 

5  2 

from  a  Thomas-Ching  model  in  which  a  total  energy  of  3  X  10  erg/cm  is 
deposited  over  8  hr  as  a  sheet  source  at  150  km.  The  above  amount  of  energy 

deposited  over  an  area  within  30°  of  the  poles  is  equivalent  to  a  ring  current 

1 8 

energy  dissipation  rate  of  7  X  10  erg/sec  over  8  hr.  Figure  96  shows  the 

ring  current  energy  dissipation  rate  of  the  same  26  May  1967  magnetic  storm. 

18 

Note  that  7x10  erg/sec  over  8  hr  is  only  a  small  fraction  of  the  total  ring 
current  loss.  From  Fig.  95  we  see  that  a  temperature  inversion  occurs  in  a 
layer  of  approximately  60  km  for  several  hours. 

In  order  to  gain  a  proper  interpretation  to  the  temperature  inversion 
shown  on  Fig.  95,  as  well  as  subsequent  figures,  several  points  should  be 
considered.  First  of  all,  the  Thomas-Ching  one -dimensional  model  calculates 
essentially  the  response  of  the  atmosphere  to  heating;  i.  e.  ,  it  predicts  incre¬ 
ments  of  temperature  profile  as  a  result  of  energy  input.  The  total  tempera¬ 
ture  profile  is  a  superposition  of  some  arbitrary  static  profile  and  the  incre¬ 
ment  profile.  In  order  to  show  that  the  temperature  increment  profile 
immediately  above  the  heating  layer  is  sufficiently  large  to  produce  a  tempera¬ 
ture  inversion  against  a  steep  static  temperature  profile,  we  have  chosen  to 
superimpose  the  increment  profile  and  the  CIRA  2  (Ref.  14)  profile,  which  has 
the  steepest  temperature  gradient.  Since  there  is  no  clean  separation  between 
the  heat  sources  maintaining  the  two  profiles,  the  temperatures  shown  on  the 
figures  may  be  too  high.  However,  the  point  is  that  our  increment  profile  is 
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ALTITUDE,  km 


SHEET  SOURCE  AT  150  km 
Q  =  3  x  I05  erg  /  cm2  OVER  8  hr 


TIME  2/4/  6  8  10  12 


Fig.  95.  Thomas -Ching  Model  Temperature  Profile  of  CIRA  2  Atmosphere 
Heated  by  Sheet  Source  at  150  km  Depositing  3  X  10®  erg/cm2 
Spread  over  Polar  Caps,  equivalent  to  7  X  1 0 1 8  erg/sec 
Energy  Deposition  Rate  over  8  hr 


Fig.  96.  Ring  Current  Energy  Dissipation  Rate  of 
26  May  1967  Storm 
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sufficiently  strong  to  maintain  an  inversion  against  the  steepest  static  CIRA 
model  temperature  profile.  Second,  the  total  amounts  of  energy  input  used 
in  the  model  calculations  represent  upper  limits  to  the  energy  input  in  the 
production  of  temperature  inversion  since  we  spread  the  total  energy  input 
over  areas  within  30°  of  both  polar  caps.  In  reality,  it  is  only  necessary 
to  produce  a  temperature  inversion  over  an  area  whose  linear  dimension 

is  a  couple  of  the  observed  characteristic  lengths.  Thus,  temperature 

1 8 

inversion  can  be  generated  by  an  energy  input  much  less  than  7  X  10  erg / 
sec  over  8  hr.  A  more  appropriate  quantity  is  energy  deposition  per  unit 
area,  which  is  shown  on  the  figures  as  the  quantity  Q.  Finally,  it  should  be 
emphasized  that  a  temperature  inversion  is  necessarily  unstable  so  that  an 
instability  occurs  to  destroy  the  temperature  inversion.  Thus,  the  occurrence 
of  an  inverse  temperature  profile  in  a  one -dimensional  model  does  not  mean 
that  such  an  inversion  exists  in  nature  but  simply  that  a  three-dimensional 
model  is  required  to  describe  the  motion  of  the  atmosphere  adequately  in  such 
a  case. 

Since  the  generation  of  our  convective  instability  is  predicated  on 
the  production  of  a  temperature  inversion  in  a  given  atmospheric  layer,  we 
could  vary  the  energy  deposition  rate  to  obtain  a  qualitative  criterion  between 
the  magnitude  of  the  storm  ring  current  loss  rate  and  the  likelihood  of  observ¬ 
ing  the  anomalous  wind.  Figures  97  through  99  show  the  temperature  profile 

predicted  by  the  Thomas -Ching  one -dimensional  model  corresponding  to  ring 

18  18  18 

current  loss  rate  of  1, 17  X  10  ,  2.  34  X  10  ,  and  4.  67  X  10  erg/sec  over 

8  hr  and  within  30°  of  the  polar  caps,  respectively.  It  is  seen  that  we  are 

unlikely  to  observe  anomalous  winds  unless  the  above  ring  current  loss  rate 

18 

is  larger  than  2.34  X  10  erg/sec.  Such  a  loss  rate  is  to  be  compared  with 
18 

7  X  10  erg/sec  of  the  historic  26  May  1967  storm,  which  generated  the 
observed  winds.  Finally,  Fig.  100  shows  that  the  heating  source  thickness 
can  be  increased  up  to  10  km  if  we  increase  the  heat  input  slightly. 

To  investigate  the  possible  mechanism  for  the  generation  of  insta¬ 
bility,  we  shall  henceforth  assume  that  an  unstable  temperature  profile  with 
a  layer  of  inverse  temperature  gradient  as  shown  in  Fig.  95  exists.  The 
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Fig.  97. 


Thomas -Ching  Model  Temperature  Profile  of  CIRA  2  Atmosphere 
Heated  by  Sheet  Source  at  150  km  Depositing  3  X  10^  erg/cm^ 
Spread  over  Polar  Caps,  Equivalent  to  1.  17  X  10^8  erg/sec 
Energy  Deposition  Rate  over  8  hr 
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Fig.  98.  Thomas -Ching  Model  Temperature  Profile  of  CIRA  2  Atmosphere 
Heated  by  Sheet  Source  at  150  km  Depositing  1X105  erg/cm^ 
Spread  over  Polar  Caps,  Equivalent  to  2.  34  X  1 0 1 8  erg/sec 
Energy  Deposition  Rate  over  8  hr 


8-9 


TEMPERATURE,  *K 


Fig.  99.  Thomas -Ching  Model  Temperature  Profile  of  CIRA  2  Atmosphere 
Heated  by  Sheet  Source  at  150  km  Depositing  2  v  10^  erg/cm2 
Spread  over  Polar  Caps,  Equivalent  to  4.67  X  1 0 1 8  erg/sec 
Energy  Deposition  Rate  over  8  hr 
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350 


A=  10  km  SOURCE,  Q  = 
B:  10 km  SOURCE,  Q -* 
C=  50km  SOURCE,  0  = 


Fig.  100.  Effect  of 


Heat  Source  Thickness 


characteristics  of  the  instability  can  then  be  described  by  the  equations  of 
continuity,  motion,  and  heat  flow  for  the  realistic  case  of  a  rotating  exponen¬ 
tial  atmosphere.  The  full  formalism  and  solutions  shall  be  exhibited  in  the 
next  section.  Since  a  number  of  the  features  turn  out  to  be  not  significantly 
different  from  that  of  the  simpler  stationary  convection,  we  shall,  for 
illustration  purposes,  first  develop  the  description  of  thermospheric  con¬ 
vection  in  a  nonrotating  exponential  atmosphere. 

The  general  set  of  equations  governing  atmospheric  motion  of  a 
layer  without  heat  sources  are  as  follows: 


Continuity: 

div  (pv)  =  -  dp/dt 

(56) 

Motion:  "v”  • 

V"v  +  ( 1/p)  Vp  -  "g  =  YV^v”  -  d'v'/d  t 

(57) 

Heat  Flow:  v  • 

VT  -  (1/pc  )  div  (KVT)  =  -  3T/3t 

P 

(58) 

3  — 

In  the  above,  pis  the  density  (gm/cm  ),  v  is  the  velocity  (cm/sec),  p  is  the 

2 

pressure  (dyne/cm  ),  and  T  is  the  temperature  (°K).  The  parameters  are 

2 

the  kinematic  viscosity  V  (cm  /sec),  heat  capacity  as  constant  pressure  c^ 

(erg/gm-°K),  and  heat  conductivity  K  (erg/cm/sec-  °K),  Throughout  the 

analysis,  y,  c  ,  and  K  shall  be  assumed  constant,  equal  to  their  values 
P 

evaluated  at  200-km  altitude.  A  more  sophisticated  treatment  would  include 
altitude  dependence  of  these  parameters  as  well  as  a  tensor  form  for  the 
kinematic  viscosity. 

At  the  onset  of  the  instability  we  assume  that  the  velocity  is  small 
so  that  the  above  equations  can  be  linearized  by  dropping  all  second-order 
terms  in  the  perturbations  of  density,  pressure,  temperature,  and  velocity. 
The  perturbed  quantities  shall  be  denoted  by  the  subscript  1,  and  the  unper¬ 
turbed  quantities  shall  be  denoted  by  the  subscript  0.  Assuming  that  the 
zeroeth-oraer  quantities  obey  the  equation  of  state 

pQ  =  (k/m)  PqT q  (59) 
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the  zeroeth-order  solutions  are  simply 


To=To-zAT/h 


(60) 


(61) 


where  Tg,  pg,  and  pg  are  constants  given  by  an  appropriate  CIRA  model.  In 
Eq.  (60),  we  have  specified  our  solution  to  one  appropriate  to  our  inversion 
region  of  thickness  h  and  maximum  inversion  AT.  The  coordinate  z  is 
measured  from  the  lower  surface,  i.e.  ,  0  <  z  <  h.  Since  AT/Tg  is  usually 
small,  e.  g. ,  Fig.  95,  we  have  ignored  the  zAT /h  term  in  Eq.  (61).  Hq  in 
Eq.  (61)  is  the  scale  height  and  is  approximately  53  km  for  our  conditions. 

For  the  perturbed  quantities,  we  note  that  the  density  and  tempera¬ 
ture  are  related  by  the  law  of  thermal  expansion 


With  the  use  of  Eqs.  (59)  through  (62),  the  first-order  equations  become, 
after  disregard  of  9p^/9t  (Refs.  42  and  44) 


9v  9v  , 

(63) 

(64) 

2_  AT  3T1 

T1  "  -p0Cp  h  vz  +  p0Cp  9t 

(65) 
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2 

where  V  is  the  Laplacian  operator.  Equations  (62)  through  (65)  now  give  a 
complete  perturbative  approach  to  instability.  The  procedure  is  to  look  for 
solutions  that  represent  a  spatial  disturbance,  i.  e.  ,  exp  (ik  •  x),  but  have  an 
exponentially  increasing  time  dependence  ewt(Rew  >  0).  Such  solutions  are 
unstable  because  the  amplitude  of  the  disturbance  will  grow  with  time  until 
the  state  of  the  system  becomes  quite  different.  These  unstable  solutions, 
however,  can  be  separated  into  two  classes  depending  on  whether  Imw  =  0  or 
Imw  i  0.  Those  with  Imw  =  0  are  convective  instabilities  because  the  fluid 
motion,  though  growing  in  time,  is  confined  to  a  volume  in  space.  Those 
with  Imw  /0  are  traveling  wave -type  instabilities  (overstability).  For  the 
incompressible  case,  Pellew  and  Southwell  (Ref.  43)  have  shown  that  only  the 
convective  mode  exists.  Chandrasekhar  (Ref.  44)  h?„s  shown,  however,  that 
the  overstable  mode  is  also  allowed  in  the  constant  density  case  if  we  include 
the  effects  of  a  rotating  atmosphere.  In  the  following,  we  shall  show  that 
these  results  are  essentially  unchanged  in  the  case  of  the  more  realistic 
exponential  atmosphere. 

Since  the  measured  wind  pattern  seems  to  be  stationary,  in  this 
section  we  shall  study,  for  illustrative  purposes,  the  simpler  convective 
instability  in  a  nonrotating  atmosphere.  Furthermore,  the  stable  and  unstable 
states  are  separated  by  a  marginal  state  of  stationary  convection,  i.  e. , 

Rew  =  0.  If  the  growth  rate  is  small,  the  characteristics  of  the  marginal 
state  are  essentially  the  same  as  the  unstable  state;  hence,  for  the  sake  of 
simplicity  we  shall  first  study  the  solutions  of  Eqs.  (63)  through  (65)  with 
d/dt  =  0.  The  complete  case  is  considered  in  the  next  section. 

With  d/dt  =  0,  we  decouple  Eqs.  (62)  through  (65)  to  obtain 

vx  =  ‘  <^/dy  (66) 


and 


v 

y 


-  d$/dy  +  b^/dx 


(67) 
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where 


=  o 

(68) 

$  =  dv  /dz  -  v  /Hn 
z  z  0 

(69) 

-PqV^U/Po)^  PoVz) 

(70) 

and 


n  dz  r0 


p0  az 


VKh 


_2 


=  vWv 

z  1  z 


(71) 


2  2  2  2  2 

In  these  equations,  V.  =  <3  /dx  +  d  /dy  .  Thus,  if  we  solve  Eq.  (71)  for  v  , 

1  z 

all  the  other  unknowns  could  be  obtained  by  solving  Eqs.  (65)  through  (70). 
The  boundary  conditions  are 


v  (z  =  0)  =  v  (z  =  h)  =  0 
z  z 


(72) 


If  our  inversion  layer  were  bounded  by  free  surfaces  then  we  would  require, 

2  2 

in  addition,  d  v  /dz  to  vanish  there  because  dv  /dz  and  dv  /dz  must  vanish 
’  z  x  y 

at  the  surfaces.  However,  there  are  no  surfaces  in  our  problem  and  because 

of  viscous  drag,  our  solution  for  v^  and  v^  must  be  joined  to  a  boundary  layer 

type  of  solution  outside  of  0  £  z  £  h.  The  boundary  condition  in  the  x  and  y 

coordinates  are,  of  course,  not  determined  by  our  problem  unless  a  particular 

convection  patte m  is  postulated.  Whatever  the  convection  pattern,  the 
2 

operator  Vj  in  Eq.  (71)  indicates  quite  clearly  that  the  solution  is  a  linear 
combination  of  exp  (iK  •  "r),  where  r  is  a  vector  in  the  x-y  plane.  Thus  the 
boundary  condition  (72)  determines  the  dominant  horizontal  wave  number  k, 
which  in  turn  is  related  to  the  size  of  the  convection  cells.  It  should  be 
pointed  out,  however,  that  a  complete  solution  of  the  convection  pattern 
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cannot  be  obtained  unless  we  know  more  about  the  horizontal  boundary 

conditions  and  the  boundary  layers  in  the  vertical  direction. 

Transforming  to  the  dimensionless  variable  4  =  z/h,  we  seek  v  in 

z 

the  form 


v 

z 


(73) 


where  f(4)  is  determined  by 


f(4)  =  -k2h2Rf(U 


(74) 


R,  the  Rayleigh  number,  is  defined  by 


(75) 


It  is  customary  in  the  treatment  of  convection  in  compressible  fluids 
(Refs.  45  and  47)  to  treat  R  as  a  constant  although,  strictly  speaking,  this 
is  incorrect  since  p^  changes  with  4.  We  shall  show  that  large  changes  in 
R  induce  only  small  changes  in  kh. 

Before  we  consider  the  solution  of  Eq.  (74),  it  may  be  instructive 
to  repeat  the  solution  of  the  incompressible  case,  which  corresponds  to  the 
infinite  Hq  limit,  i.e., 

-  k2h2  j3  f(|)  =  -k2h2Rf(|)  (76) 


The  general  solution  to  Eq.  (76)  is  a  linear  combination  of  cosh  and 
sinh  p£,  where  pis  determined  by  the  algebraic  equation 


Since  the  various  cosh  and  sinh  solutions  are  linearly  independent,  the 
solution  satisfying  the  boundary  condition  is 

£( 4 )  =  C  sin  mTT^  ,  m  i  0  (78) 

i.  e.  ,  (jl  =  iimir.  Thus,  in  terms  of  the  wave  length  in  units  of  h,  one  finds 

H  =  2ir/kh 

Equation  (77)  becomes  a  relation  between  the  convection  cell  size  qh  and  the 
Rayleigh  number  R,  which  depends  on  the  local  properties  of  the  atmosphere 


R  - 


(79) 


The  case  for  m  =  1  is  plotted  as  Curve  C  on  Fig.  101.  Note  that  no  instability 
can  occur  unless  R  is  greater  than  R  critical  where 

F  critical  =  27  m^^/4  (80) 

The  critical  Rayleigh  numbers  for  m  =  1  and  2  are  also  shown  on  Fig.  101. 

The  case  of  an  exponential  atmosphere  is  considerably  more  com¬ 
plex.  We  can  seek  a  solution  f(4)  to  Eq.  (74)  in  the  form 

f(|)  =  C  e^  sin  miT^  (81) 
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which  satisfies  the  boundary  conditions  (72).  Substitution  of  Eq.  (81)  into 

i 

Eq.  (74)  yields  two  equations  from  which  X  and  kh  can  be  determined. 
These  are 


/,  h  \r/x2  2  2,  2,  2,2  .  2  2.2, 

(X  -  j[(X  -  m  tt  -  k  h  )  -  4m  it  X  ] 


,  -,w,2  2  2  ,2,  2, 

+  2X(  X  -  m  tt  -kh) 


("'Ho) 


22  ,  2,  2 

-mu  -  k  h 


=  G(  X)  =  0  (82) 


r,,2  2  2  2,2  .  2  2,2, 

[(X  -  m  tt  -  a  )  -  4m  tt  X  J 


Kf 


2  2  ,2,2 

-  m  n  -  k  h 


a  2  2, 

8m  tt  X 


(x-£ 


2  2  2  2  2  2  2 
(X*  -  m  n  -  kTi  )  =  -RkV  (83) 


Equation  (82)  determines  X  as  a  function  of  k,  and  Eq.  (83)  determines  the 
relation  between  cell  size  and  the  Rayleight  number.  It  is  a  simple  matter 
to  show  that  Eq.  (82)  has  five  real  roots  for  X.  The  bounds  for  these  roots 
are  given  by 


.  ,  2  2,  2,  2,1/2 

+  (mTr  +  k  h  )  <  xi  <  <* 


,  ,  ,  _  ,  .  2  2,,  2,  2,1/2 

h/Hp  £  X^  ^  +  (m  tt  +kh) 


0  <  X3  <  h/HQ 


(m2TT2  +  k2h2)1/2  <  X4  <  -h/HQ 


,  ,  ,  -  2  2  ,  ,2,  2,1/2 

-  co  <  X,.  <  -  (m  it  +  k  h  ) 
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9 


22  221/2 

These  inequalities  are  derived  under  the  assumption  1i/Hq  <  (m  tt  +k  h  ) 

which  is  quite  well  satisfied  because  1  £  1i/Hq  £  2  in  our  case.  Since  the 

observations  do  not  show  a  strong  dependence  of  maximum  wind  speed  with 

altitude,  solutions  X,,  X-,  X,,  and  \r  can  be  excluded.  Solution  X_,  on  the 

1  2  4  5  3' 

other  hand,  seems  to  be  a  good  candidate  for  the  physical  situation,  since  it 
gives  an  increase  of  maximum  wind  speed  of  less  than  one  e-fold.  Under  the 
approximation 


uixj  „  ,  2  2  ,  .2,  2.1/2 

Ii/Hq  «  (m  tt  t  k  h  ) 


we  obtain 


(h/H0)(m2TT2  +  k2h2) 

3(m2iT2  +  k2h2)  -  2 f h / H 0 ) 2 


(84) 


(85) 


Figure  102  shows  the  rigorous  solution  for  G(X)  =  0  for  the  case  m  =  1  and 
2  2  2 

k  h  «tt  .  It  is  seen  that  Eq.  (85)  is  a  good  representation  of  X  even  if 

2  2  2  ^ 

h/Hp  2.  We  shall  see  presently  that  k  h  «  tt  is  an  exceedingly  good 
approximation  for  realistic  values  of  R  in  our  case. 

Using  the  definition  of  n  and  the  representation  of  X^  given  by 
Eq.  (85),  we  obtain  the  relation  between  q  and  R  from  Eq.  (83).  Curves  A 
and  B  in  Fig.  101  show  this  relationship  for  h/H^  equal  to  1  and  2  in  the 
m  1  mode,  respectively.  Comparison  with  the  constant  density  case 
(Curve  C)  shows  that  the  exponential  atmosphere  case  for  solution  X^  differs 
little  from  it.  For  other  solutions  and  Xc),  the  differences 

between  compressible  and  incompressible  atmospheres  are,  however,  quite 
drastic.  The  existence  of  the  solution  is  not  surprising  for  it  is  well  known 
that  compressible  hydrodynamics  often  has  a  mode  of  motion  approximating 
incompressible  hydrodynamics. 
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Fig.  102.  Solution  of  G(\)  =  0  for  \  Near  Zero  and  m  =  1 
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D. 


THERMOHYDRODYNAMIC  INSTABILITY  IN 
ROTATING  THERMOSPHERE 


In  the  previous  section,  we  have  examined  the  features  of  steady 
convection  in  a  stationary  exponential  atmosphere.  In  this  section,  we  shall 
examine  the  more  realistic  case  of  a  rotating  exponential  atmosphere. 

Further,  the  question  of  growth  rate  and  the  existence  of  growing  solutions 
will  be  studied  in  more  detail. 

The  first-order  perturbation  equations  corresponding  to  a  rotating 
exponential  atmosphere  within  a  layer  of  inverse  temperature  profile  are 
given  by  Eqs.  (63)  to  (65)  with  the  slight  modification  that  a  Coriolis  force 
term  Zv  X  IT  is  added  to  the  right-hand  side  of  Eq.  (64)  corresponding  to  the 
effects  of  an  atmosphere  rotating  with  angular  velocity  vector  il.  Since  we 
are  interested  in  local  winds  in  the  polar  regions,  S2  shall  be  assumed  parallel 
to  the  z  axis,  and  effects  due  to  curvature  of  the  earth  shall  be  ignored. 

Linear  dependence  on  3/3t  of  the  variables  v,  T^,  and  p^  introduces  the 
common  time  factor  c  (w  complex)  in  our  equations.  The  resultant  set  is 


3v  3v  , 

x  .  y  ,  1  3  ,  . 

9x  9y  pQ  3z  ^0Vz 


0 


(86) 


_l!Zi 

P0  3* 


w)v  +  2S2v 
x  y 


(87) 


1  dPl  2 

~a —  =  (Y^7  '  w)v  "  2S2v  (88) 

p0  dv  Y  y  x 


=  <YV2  -  w)vz  +  gpT  j  (89) 
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(90) 


._2.„  _  at 

'  1  wPocPTi  =  -pocP  h  vz 

”/>'o 

With  t  e  and  the  usual  assumption  of  constant  Rayleigh  number,  a 

tedious  reduction  gives  an  equation  for  v 


(91) 


Proceeding  as  before,  we  consider  solutions  of  the  form 


vvt  ik*  r  ... .  t  ,,  .  . 

\  z  -  e  e  f(|)  ,  |  =  z/h  ,  r  =  (x,  y) 

and  obtain  the  dimensionless  equation  for  f(|) 


,  2  2 
+  4S2q 


f(|) 


qw 


a2f(^) 


(92) 


2  2 

where  the  parameters  are  defined  by  p  =  h  p^c^/K,  q  =  h  /y,  s  =  h/H^,  and 
a  =  hk.  It  is  a  simple  but  tedious  matter  to  show  that  Eq.  (92)  has  a  solution 


of  the  form 


f(|)  =  e ^  sin  muf,  (93) 

One  then  obtains  two  conditions  relating  the  unknowns  a,  \,  and  w.  These 
relations  are  not  particularly  enlightening  because  of  their  complexity;  thus 
we  shall  consider  the  various  limit  situations  in  order  to  study  the  effect  of 
various  te rms. 
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First,  we  consider  the  question  of  the  existence  of  the  overstable 
mode.  In  the  case  of  a  nonrotating  atmosphere,  i.  e.  ,  n  vanishes,  then 
Eqs.  (92)  and  (93)  yield  two  conditions 

2  2  2  2  .<.,2  2  2  2  .r,,  .2  2  2  2,  „  2  2W,  ,» 
(X  -  m  tt  -  a  -  wp)j  (\  -  m  tt  -a  -  wq)[(\-  s)  -  m  tt  -a  ]  -  4m  n  X(  X  -  s ) } 

-  4m2T72X  j  X[  (  X  -  s )2  -  m2  tt  2  -  a2 ]  +  (  X  -  s )( X2  -  m2n2  -a2  -  wq)  (  =  -Ra2  (94) 

and 

.  r..  ,2  2  2  2ir~/.2  2  2  2.  ,  .  . , 

X[(  X  -  s)  -  m  tt  -  a  J[2(X  -  m  tt  -a)+w(-p-q)J 

,  /»  .r/,2  2  2  2  . .  .2  2  2  2  2  2.2,  _  /nc, 

+  (X  -  s)[(X  -  m  tt  -  a  -  wq)(  X  -  m  tt  -  a  -  wp)  -  4m  tt  X  |  =  0  (95) 

From  Eqs.  (94)  and  (95),  we  obtain 

X[w(-p  -  q)  +  2(X2  -  m2r2  -  a2)]{f(X-  s)2-  m2r2  -  a2]2  +  4m2n2(X-s)2  =  Ra2(\-s) 

(96) 

Since  all  parameters  are  real  in  Eq.  (96),  w  cannot  be  complex  when  Eq.  (96) 
is  considered  as  an  equation  for  w.  This  implies  that  no  overstable  modes  are 
allowed  in  an  exponential  atmosphere  without  rotation  because  an  overstability 
necessarily  requires  that  Rew  >  0  and  Imw  t  0.  On  the  other  hand,  if  Q  t  0, 
the  conditions  on  w  would  be  quartic  in  w  co  that  complex  solutions  with 
Rew  >  0  and  Imw  i  0  are  possible.  While  the  actual  case  involves  coupled 
quartic  equations  for  w,  the  X— 0  limit  leads  to  simpler  equations  which 
indeed  have  a  solution  Rew  >  0  and  Imw  i  0.  Similar  results  were  obtained 
by  Chandrasekhar  and  Pellew  and  Southwell  by  a  somewhat  different  analysis. 

Next,  we  note  that  the  existence  of  solutions  with  Rew  >  0  has  been 
tacitly  assumed  in  the  previous  section.  In  order  to  prove  the  existence  of 
such  growing  solutions  in  the  nonrotating  case,  we  consider  Eq.  (96)  together 
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with  the  physical  solution  \  =  l/3s.  Since  the  expression  in  braces  is 
2  2  2  2 

positive  and  m  tt  +  a  >  X.  ,  it  is  easy  to  see  that  growing  solutions  can  be 
obtained  in  either  of  two  cases 


Case  1 :  -p  -  q  <  0 

and 

Ra2(  X.  -  s) 

> 

2\(X2  -  nA2  -  a2)|  *  •  •  [  | 

Case  2:  -p  -  q  >  0 

and 

Ra2(  X.  -  s) 

< 

2X(X2  -  m2ir2  -  a2){-  •  •  }| 

In  the  case  of  the  thermosphere  at  160  to  220  km,  -p  -  q  <  0  holds, 
hence  condition  1  must  be  satisfied.  Using  the  observed  value  of  q  =  18, 
i.e.  ,  a  =  1/3  and  values  of  parameters  given  in  the  next  section,  we  obtain 
the  growth  time 


t 


1.1  (|-R  -  2 ,  1 1 5 j  1  hr 


where  R  is  measured  in  units  of  1000.  Hence,  if  R  a  9500  is  contained  in 
our  range  of  estimated  R  values,  the  observed  1  to  2  hr  growth  time  can  be 
accounted  for. 

Next,  we  shall  consider  the  effects  of  atmospheric  rotation  on  the 

relation  between  the  horizontal  characteristic  length,  a,  and  the  Rayleigh 

number.  For  this  purpose,  we  study  two  limits  of  Eq.  (92),  namely  the  case 

of  small  growth  rate,  i.  e.  ,  w  —  0  and  the  case  of  large  growth  rate,  i.  e. , 

2 

w  »  a  . 

The  w  —  0  case  yields  the  conditions 


o\/\2  2  2  2,r,. 

2M  X.  -  m  tt  -  a  )[ (  X 


.2  2  2  2 ,  .  ,r,v  2  2  2  2.2  ..2  2  2  ,  _2,  _ 

s)  -mir  -  a  J  +  (  X.  -  s  )[(\  -  m  tt  -a  )  -  4\  m  tt  +Q  J  =  0 


and 

r/  .2  2  2  2,2  .,2  2  2ir , .  ,2  2  2  2, 

[(X.  -  m  h  -  a)  -4\mTr][(\-s)  -  m  tt  -a] 

-8m2TT2k(\-  s)(\2  -  mZTT2  -  a2)  +  Q2[(\-  s)2  -  nA2]  =  -Ra2  (98) 
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2  9  2 

where  the  dimensionless  parameter  Q  is  2fth  /  X.  =  5  for  X  =  10  cm  /sec  at 

200-km  altitude.  Equations  (97)  and  (98)  are  to  be  compared  to  Eqs.  (82)  and 

2  4 

(83)  for  then— 0  limit.  Since  Q  «(mir)  ,  the  solution  X^  for  Eq.  (97)  is 
essentially  unchanged,  i.e.  ,  X  »  l/3s.  The  relation  between  R  and  n, 

Eq.  (98),  is  shown  on  Fig.  103.  Thus,  for  very  small  growth  rate,  the 
characteristic  length  of  the  convective  instability  is  not  essentially  changed 
even  if  we  include  the  effects  of  atmospheric  rotation. 

For  the  case  of  large  growth  rate,  the  situation  becomes  much  more 
complicated.  The  essential  features,  however,  can  be  studied  in  the  constant 
density  limit.  Letting  s  vanish  in  Eq.  (92),  we  obtain  the  solution 

f(4)  =  sin  mTi| 


and  the  condition 


2  2  2  >r  2  2  2  .2,  2  2  2.  ^2  2  2,  „  2,  2  2  2 
(m  tt  +  a  +  pw)[(m  tt  +  a  +qw)  (m  u  +a  )  +  Q  m  i  ]  =  Ra  (m  n  +a  +  qw) 


(99) 


For  values  of  ft  and  q  appropriate  to  the  thermosphere  near  200-km  altitude, 
2  4  4 

Q  «  m  it  is  very  well  satisfied;  hence,  we  lose  no  generality  in  simplifying 

2 

Eq.  (99)  by  dropping  the  Q  term.  In  terms  of  the  dimensionless  charac¬ 
teristic  length  n=  2Tr/a,  one  finds 


R 


16tt4 

4 


2  2 

m  it  +  pw 
4„2 


1,2  +  sw  n2)(i 


(100) 


This  is  to  be  compared  to  Eq.  (79),  which  is  a  special  case  of  Eq.  (100)  with 
w  =  0. 

Since  pw  and  qw  are  both  positive,  Eq.  (100)  indicates  that  for  a 
given  R  the  cell  size  q  obtained  is  smaller  than  that  given  by  Eq.  (79).  For 
the  observed  growth  rate  of  approximately  1.  5  to  3  hr,  the  reduction  of  q  due 
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to  this  effect  is  quite  small  since  pw  and  qw  are  both  smaller  than  tt  .  On 
the  other  hand,  if  X.  i  0,  this  tendency  can  be  compensated  for  in  the  opposite 
direction,  depending  on  the  value  of  h/Hg.  Large  uncertainties  in  the  param¬ 
eters  making  up  R  do  not  make  a  determination  of  growth  rate  and  cell  size 
profitable,  although  the  previous  consideration  does  show  that  1.  5  to  3  hr 
growth  time  is  quite  consistent  with  q  ~  20  as  is  observed. 

E.  CHARACTERISTIC  FEATURES  AND  OBSERVATION 

At  the  end  of  the  previous  section,  we  show  that  the  instability 
growth  time  can  be  accounted  for.  In  this  section,  we  shall  compare  other 
characteristic  features  of  our  theory  with  the  observations.  The  first  and 
crucial  question  that  has  to  be  answered  is,  of  course,  whether  a  convectively 
unstable  solution  can  be  maintained  by  a  reasonable  amount  of  energy  deposi¬ 
tion.  Part  of  the  answer  lies  in  the  production  of  a  temperature  inversion, 
which  is  model-dependent.  Figure  96  shows  that  an  inversion  can  be  main¬ 
tained  in  at  least  one  atmospheric  model  by  an  energy  deposition  roughly 
equal  to  one  third  the  energy  loss  from  the  ring  current.  Given  the  existence 
of  a  temperature  inversion,  the  criterion  of  instability  is  given  in  the 
incompressible  case  by 


R  = 


gyoh4  /Ap 

KhV  (tJ 


> 


4  4 

27m  tt 

4 


(101) 


The  small  difference  between  the  incompressible  case  and  the  solution  of 
the  exponential  atmosphere  can  be  ignored.  From  Fig.  96  it  is  seen  that  the 
first  mode  is  excited  at  R  >  65  7  and  the  second  mode  is  excited  at  R  >  10600, 
On  the  other  hand,  our  range  of  R  values  must  encompass  R  =  9500  in  order 
to  account  for  the  observed  growth  time. 

Tc  see  if  instability  can  be  excited,  we  shall  use  the  atmospheric 
parameters  appropriate  to  200 -km  altitude  during  extremely  high  geomagnetic 
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activity  (CIRA  10).  c^,  K  are  given  in  terms  of  pg  and  Hg  (Ref.  48).  The 

kinematic  viscosity  y  is  given  by  Kohl  and  King  (Ref.  49).  We  have 

0^  =  107  erg/gm-°K 

K  =  1.  06  X  10^  erg/cm-sec-°K 
9  2 

y  =  10  cm  /sec 

From  Fig.  95,  AT/Tg  =  1/10.  If  we  let  Pg  vary  over  the  range  of  altitudes 
from  160  to  220  km,  the  Rayleigh  number  is  bounded  by 

4250  <  R  <  22,000 

with  a  typical  value  at  190  km  of  R  =  9500.  Thus,  the  parameters  favor  the 
generation  of  the  first  mode  of  the  instability.  The  existence  of  the  second 
mode  of  instability  is  possible,  but  the  characteristic  lengths  are  much 
smaller. 

If  we  are  to  interpret  the  observations  in  terms  of  the  excitation  of 
convection  in  the  first  mode,  we  find  that  the  large  horizontal  characteristic 
length  can  be  explained.  From  the  curves  on  Fig.  101,  we  determine  the 
horizontal  characteristic  length  to  be  between  960  and  2000  km  for  the  above 
variation  of  R.  If  the  second  mode  is  excited  the  characteristic  length  is  of 
the  order  150  km.  On  the  other  hand,  if  we  include  the  effect  of  growth  rate 
a.  decrease  of  10  percent  in  the  above  figures  is  possible  for  w  slightly  above 
2  hr.  As  is  evidenced  in  Fig.  103  and  Eq.  (99),  the  decrease  in  horizontal 
characteristic  length  due  to  the  effects  of  atmospheric  rotation  is  very  small. 

Next,  without  constructing  the  full  set  of  solutions,  we  can  also  find 
the  magnitude  ratio  between  horizontal  and  vertical  velocity  components.  This 
is  easily  obtained  by  substitution  of  the  solution  for  v 

v  =  v  n  e^£  sin  e^*  r  (102) 

z  zO 
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into  the  equation  of  continuity  (63).  Then  it  is  evident  that 


v  _  v  - 
xO  _  yO  tt  _  g 

\t>  ‘  vz0  ’  kh  '  2 


(103) 


From  Fig.  96  again,  q  is  bounded  by  16  £  *1  £  30,  so  we  have 

%0/vz0  '  10  <104> 

Equation  (104)  agrees  quite  well  with  observation.  Thus,  we  see  that  the 
most  significant  features  of  the  observations  are  quite  consistent  with  the 
interpretation  of  lowest  mode  steady  convection  generated  by  local  heating. 

Finally,  one  must  examine  the  validity  of  comparing  convection 
characteristics  of  a  low  velocity  perturbation  theory  with  observed  charac¬ 
teristics  of  an  extremely  high  speed  wind  pattern.  To  this  end,  we  note  that 
our  predicted  features  originate  from  solution  for  from  Eqs.  (82)  and  (83), 
which  as  we  recall,  also  hold  if  the  growth  rate  of  the  instability  is  small. 

This  immediately  means  that  the  characteristic  length  L  is  essentially  time- 
independent.  Further,  the  time  dependence,  ewt  (Imw  =  0),  in  vx^  and  v^g 
cancel  in  Eq.  (103),  so  the  speed  ratio  Eq.  (104)  remains  unchanged.  On  the 
other  hand,  if  the  growth  rate  is  large,  the  effect  is  to  decrease  the 
characteristic  length. 

F.  FULLY  DEVELOPED  CONVECTIVE  INSTABILITY 

In  previous  sections  we  have  attempted  to  interpret  features  of  the 
observed  wind  in  terms  of  theoretically  predicted  features  of  a  convective 
instability  analysis  valid  during  the  onset  phase  when  the  velocities  are  small. 
Although  many  of  the  features  show  good  agreement  between  theory  and  experi¬ 
ment,  we  still  must  be  cautious  about  the  validity  of  the  comparison  since  the 
observed  wind  speed  components  are  so  large  that  the  nonlinear  •  V\T 
dominate.  Thus,  in  order  for  our  interpretation  to  be  valid,  it  is  necessary 
to  show  that  solutions  of  the  type  given  by  Eq,  (102)  exist  in  the  case 
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2  — 

yV  v«v  •  y  v,  so  that  the  characteristic  features  of  such  solutions  can  be 
maintained  when  the  amplitude  of  a  disturbance  grows  to  large  values. 

In  compressible  hydrodynamic  flow,  if  the  velocity  is  close  to  the 
velocity  of  sound,  compressibility  becomes  important,  and  the  number  of 
degrees  of  freedom  rapidly  increases  until  turbulence  sets  in.  Qualitatively, 
the  onset  of  a  turbulent  state  as  distinguished  from  a  well-organized  state  of 
motion  is  determined  by  a  critical  Reynolds  number  N  =  lv/y  where  1  is  the 
dimension  of  the  problem,  v  is  the  observed  velocity,  and  y  is  the  viscosity. 
For  our  maximum  speeds  N  =  1000.  If  incompressible  flow  between  two 
planes  a  distance  h  apart  is  any  guide,  the  critical  Reynolds  number  for  onset 
of  turbulence  is  N  =  7700  (Ref.  50).  Thus,  a  steady,  organized  motion  is 
quite  possible  even  for  wind  speeds  of  the  order  of  1  km/sec,  although  the 
sonic  speed 


c  =  (,gH0)1/2 


(105) 


with  3/2  <  v  £  5/3  near  200  km,  is  approximately  800  m/sec. 

Assuming  the  existence  of  organized  motion,  we  can  obtain  an 
approximate  set  of  equations  valid  for  the  high  speed  case.  We  shall  assume 
that  pressure,  density,  and  temperature  changes  are  still  small.  Ignoring 


dv  dv  ,  , 

_ x  _ y  Id,  . 

"dx~  +  ~dy  +  dz”  (p0Vz) 


0 


(106) 


Vv  +—  Vp 

po 


+  ?PT .  =0 


(107) 


v  •  VT  =  0  (108) 

In  Eq.  (107),  we  shall,  in  addition,  ignore  v  (dv  /dz)  in  accordance  with  the 

z  z 

observations. 
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Although  there  is  little  resemblance  between  Eqs.  (107)  and  (63)  or 
between  Eqs.  (108)  and  (64),  the  above  set  of  equations  have  the  following 
partial  solutions  which  are  akin  to  Eq.  (73) 


v  =  u  (z)  sin  k  •  r  ,  r  =  (x,y) 
z  z 


(109) 


V 

u  (z) 

X 

X 

V 

u  (z) 

y 

y 

cos  k  •  r 


(110) 


»1 

’p(z)" 

T>. 

t(z) 

cos 


(111) 


Substitution  of  Eqs.  (109)  through  (111)  into  Eqs.  (106)  through  (108)  gives 


ku  +kuu  -  u  -  t  —  +  — —  p  =  0 
xx  y  x  y  zdz  Pf 


kuu  +  ku  -u  - ,  p  =  0 

x  x  y  yy  zdz  p  * 


du  2k 
0 

du  2k 

% 


(112) 


(113) 


kuu  +  kuu  +  — —  4^  -  gBt  =  0 
x  x  z  y  y  z  pn  dz 


(114) 


k  u  +  k  U  =  -—-4“  (Pnu  ) 
XX  y  y  p0  dz  ^0  z 


(115) 


2t(k  u  +  k  u  )  =  u  r 
xx  y  y  z  dz 


(116) 
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(117) 


Defining  w  =  p.u  ,  we  find  that  Eqs.  (114)  and  (115)  yield 

u  Z 


t 


=  ±<tQ/w2)w 


2 


where  t^  =  max  (|t|  )  and  w^  =  max(w),  Now  Eqs.  (112)  through  (116)  can  be 
decoupled  to  obtain  a  nonlinear  differential  equation  for  w 


'  w„  ' 


w  = 


dw  d  w 

dz  ,  2 
dz 


w 


H0  1  2 
0  dz 


+  2k 
H 


dw 

dz 


(118) 


where  k  =  k^  +  k^  as  before.  Similar  to  the  determining  equation  (74), 

Eq.  (118)  is  to  be  solved  with  boundary  conditions  w(0)  =  w(h)  =  0.  The 
2 

coefficient  of  w  on  the  left-hand  side  shall  be  considered  constant. 

The  general  solution  of  Eq.  (118)  is  very  difficult  because  of  the 

2 

high  degree  of  nonlinearity.  In  our  case  k«l/h,  hence  the  term  2k  (dw/dz) 

can  be  ignored.  On  the  other  hand,  if  H_  becomes  large  compared  to  h/2,  the 

2  U  2  2 

term  (l/HgUdw/dz)  also  becomes  small  compared  to  (2/HQ>(d  w/dz  ).  Hence, 

if  such  asymptotic  approximations  hold,  we  have  a  simpler  equation 


± 


dw  d2w  d3w  2w  d2w 
dzdz2  dz3  »D  dz2 


(119) 


A  particular  solution  of  Eq.  (119)  satisfying  w(0)  =  w(h)  =  0  is. 


w(z)  =  Wq  sin  mwz/h 


m  =  1,  2, 


where  m  satisfies  the  relation 


2  2  k  h  g^H0p0t0 
m  w  =  - 

W0 


,  m  =  1,  2,  *  *  * 


surprisingly, 

(120) 

(121) 
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Using  Wq  =  max  (pgUz)  and  assuming  Pq  to  be  approximately  constant,  we 
have  from  Eq.  (121)  for  the  lowest  mode,  m  =  1, 


max(u  ) 
z 


■( 


k2h2 


1/2 


tngPH 


0 


(122) 


From  our  previous  sections,  we  had  shown  that  kh/rr  «1;  hence  max(uz)«c 
even  in  the  case  of  fully  developed  instability  when  tg  —  Tg  if  we  compare 
Eq.  (122)  with  Eq.  (105). 

Next,  from  the  equation  of  continuity  (106),  and  Eq.  (120)  we  obtain 
an  estimate  of  the  horizontal  velocity  component  for  m  =  1 


max(u*,y)  =  (l^gHo) 


1/2 


(123) 


Thus  we  see  that  our  horizontal  wind  speed  components  approach  the  speed  of 
sound  when  tg/Tg  approaches  v.  In  such  a  case,  is  no  longer  small  com¬ 
pared  to  T„.  Nevertheless,  maxfu  )  can  be  a  substantial  fraction  of  the 
v  0  x,y 

sound  speed. 

Finally,  it  is  clear  from  Eqs.  (109)  and  (120)  that  a  solution  of  the 
convective  type  does  exist  for  the  highly  nonlinear,  high  wind  speed  case. 
Although  we  have  only  shown  that  the  solution  exists  in  the  asymptotic  limit 
OfHg  »h/2,  the  full  equation  (118),  together  with  the  consistency  condition 
Wq  =  max(w)  can  be  solved  numerically  to  obtain  w(z)  and  Wq.  Further, 

Eq.  (123)  shows  that  our  solution  is  valid  for  cases  when  the  wind  speed  is  a 
substantial  fraction  of  the  speed  of  sound. 
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APPENDIX  A 


ACCELEROMETER  DATA 


The  data  for  Revs.  5  through  67  are  presented  in  Figs.  A-l  through 
A-63,  respectively.  The  actual  data  points  are  indicated  by  a  cross  and  the 
dubbed  points  by  triangles.  The  data  in  the  plots  have  been  converted  to  fps^ 
by  multiplying  the  instrument  counts  by  -g  X  10“^  (g  =  32.174  fps^);  the  abcissa 
is  TRG  time  in  seconds. 

The  time  tag  on  the  accelerometer  data  was  referenced  to  a  clock  on 
board  the  spacecraft,  0.0  sec  being  some  time  after  launch.  In  order  to 
adjust  the  data  time  on  the  abcissa  of  the  plots  (TRG  time)  to  universal  time, 
a  time  bias  of  67,  507. 0  sec  must  be  applied  to  the  data.  Referenced  to 
0.0  hr  of  the  launch  date  22  May  1967,  the  universal  time  is 


t 


universal 


=  TRG  +  67,507.0  sec 


The  data  shown  in  this  report  are  available  on  cards  and  on  magnetic 
tape  in  various  forms.  One  of  the  forms  available  on  magnetic  tape  is  time- 
sequenced  card  images  in  a  format  which  is  readily  convertible  to  computer 
systems  other  than  the  Aerospace  system.  The  tape  was  written  at  556  bits 
per  inch  in  the  Binary  Coded  Decimal  (BCD)  mode  in  one  data  file,  and  one 
physical  record  consists  of  80  Binary  Coded  Integer  (BCI)  characters 
(240  bytes).  The  format  for  each  record  is  as  follows: 


Columns 


Contents 


1  -  16  Blank 

17  -  23  TRG  time  reference  in  seconds  (integer  right  justified) 


A-  1 


Columns 


Contents 


24  -  25  Blank 

26  -  39  Accelerometer  counts  (variable  of  the  form  E  14.7) 

40  -  79  Nonessential 

80  1,  2,  3,  or  4  (1,  2,  3  indicate  actual  data;  4  indicates 

dubbed  data) 

The  accelerometer  data  on  this  tape  have  been  processed  according 
to  the  method  outlined  in  Volume  I  of  this  report.  Since  the  data  appear  as 
instrument  counts,  they  must  be  multiplied  by  -g  X  10"?  to  convert  to  fps^0 
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Fig.  A-24.  LOGACS  Accelerometer  Data  for  Rev.  28 
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Fig.  A- 30.  LOGACS  Accelerometer  Data  for  Rev.  34 
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Fig.  A-34,  LOGACS  Accelerometer  Data  for  Rev,  38 
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Fig.  A-48.  LOGACS  Accelerometer  Data  for  Rev.  52 
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Fig.  A-62,  LOGACS  Accelerometer  Data  for  Rev.  66 
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Fig.  A-63.  LOGACS  Accelerometer  Data  for  Rev.  67 
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APPENDIX  B 


DENSITY  DATA 


Densities  reduced  from  the  LOGACS  data  were  processed  through 
The  Aerospace  Corporation  Computer  Printout  of  Densities  (CPD)  program, 
to  interpolate  for  density  values  at  preselected  values  of  time,  altitude, 
density,  and  other  parameters.  CPD  accepts  as  input  the  accelerometer 
data  (fps^)  and  an  ephemeris,  outputs  the  density  derived  from  these  two 
quantities,  computes  the  values  of  solar  and  geomagnetic  indices  from  tabular 
input,  and  compares  LOGACS  densities  with  densities  from  selected  model 
atmospho  res . 

To  study  the  behavior  of  the  lower  altitude  portions  of  the  LOGACS 
density  data,  we  prepared  semi-logarithmic  plots  of  density  versus  altitude. 
These  plots  cover  an  altitude  range  from  approximately  140  to  250  km.  The 
data  were  arranged  in  files  so  that  two  plots  per  rev.  were  obtained:  one 
toward  perigee  (250  to  140  km)  and  one  away  from  perigee  (140  to  250  km). 
The  data  for  Revs.  5  through  67  (both  toward  perigee  and  away  from  perigee) 
are  presented  in  Figs.  B- 1  through  B-98,  respectively. 


B-l 


\w. - rar - tw, — m-. - 22 v.  2<an mr 

altitude  c KH) 


.  B-l.  LOGACS  Density  Data  for  Rev.  5,  Toward  Perigee 
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Fig.  B-2.  LOGACS  Density  Data  for  Rev.  5,  Away  from  Perigee 
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Fig.  B-3.  LOGACS  Density  Data  for  Rev.  6,  Toward  Perigee 
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Fig.  B-4.  LOGACS  Density  Data  for  Rev.  6,  Away  from  Perigee 
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Fig.  B-5.  LOGACS  Density  Data  for  Rev.  7,  Toward  Perigee 
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6.  LOGACS  Density  Data  for  Rev.  7,  Away  from  Perigee 
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Fig,  B-8.  LOGACS  Density  Data  for  Rev,  9,  Away  from  Perigee 
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Fig.  B-12.  LOGACS  Density  Data  for  Rev.  11,  Away  from  Perigee 
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g.  B-16.  LOGACS  Density  Data  for  Rev.  1  3 ,  Away  from  Perig 
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Fig.  B-17.  LOGACS  Density  Data  for  Rev.  14,  Away  from  Perige 
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Fig.  B-18.  LOGACS  Density  Data  for  Rev.  18.  Toward  Perigee 
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Fig.  B-20.  LOGACS  Density  Data  for  Rev.  20,  Away  from  Perigee 
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Fig.  B-2t  .  LOG  AC’S  Density  Data  for  Rev.  21,  Toward  Perigee 
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Fig.  B-22.  LOGACS  Density  Data  for  Rev.  21,  Away  from  Perigee 
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Fig.  B--2^.  I.OGACS  Density  Data  for  Rev.  22,  Toward  Perigee 
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Fig.  B-24.  LOGACS  Density  Data  for  Rev.  24,  Toward  Perigee 
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Fig.  B-25.  LOGACS  Density  Data  for  Rev.  24,  Away  from  Perigee 
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Fig.  B-27.  LOGACS  Density  Data  for  Rev.  25,  Away  from  Perigee 
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Fig.  B-28.  LOGACS  Density  Data  for  Rev.  26,  Toward  Perigee 
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Fig.  B-29.  LOGACS  Density  Data  for  Rev.  26,  Away  from  Perig 
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-30.  LOGACS  Density  Data  for  Rev.  27,  Toward  Perigee 


I  ip.  B-31.  LOGACS  Density  Data  for  Rev.  27,  Away  from  Perigee 
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Fig.  B-32.  LOGACS  Density  Data  for  Rev.  28,  Toward  Perigee 
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Fig.  B-33.  LOGACS  Density  Data  for  Rev.  29,  Toward  Perigee 
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Fig.  B-35.  LOGACS  Density  Data  for  Rev.  30,  Towa 
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g.  B-40.  LOGACS  Density  Data  for  Rev.  32,  Away  from  Perigee 
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Fig.  B-43.  LOGACS  Density  Data  for  Rev.  34,  Toward  Perigee 
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Fig.  B-44.  LOGACS  Density  Data  for  Rev.  34,  Away  from  Perigee 
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Fig.  B-45.  LOGACS  Density  D 
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B-47.  LOGACS  Density  Data  for  Rev.  36,  Away  from  Perigee 
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Fig.  B-48.  LOGACS  Density  Data  for  Rev.  37,  Toward  Perigee 
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Fig.  B-50.  LOGACS  Density  Data  for  Rev.  ^8,  Toward  Perigee 


B-51 


ALT  I  TUDE  C  K.H) 


Fig.  B-52.  LOGACS  Density  Data  for  Rev.  39,  Toward  Perigee 
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Fig.  B-53.  LOGACS  Density  Data  for  Rev.  39,  Away  from  Perigee 
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Fig.  B-55.  LOGACS  Density  Data  for  Rev.  40,  Away  from  Perigee 
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B-58.  LOGACS  Density  Data  for  Rev.  43,  Toward  Perigee 
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Fig.  B-59.  LOGACS  Density  Data  for  Rev.  43,  Away  from  Perigee 
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Fig,  B-6?.  LOC. ACS  Density  Data  for  Rev.  45,  Away  from  Perigee 
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.  B-64.  LOG  ACS  Density  Data  for  Rev.  46,  Away  from  Perig 


Fig.  B-65.  LOGACS  Density  Data  for  Rev.  47,  Toward  Perigee 
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Fig.  B-66.  LOGACS  Density  Data 
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Fig,  R-76.  LOGACS  Density  Data  for  Rev.  S3,  Away  from  Perigee 


B-77 


777 


140. 160. 180. 

•»  ALTITUDE  (KM) 


2407 


Fig.  B-77,  LOG  ACS  Density  Data  for  Rev.  54,  Inward  Perigee 
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Fig.  B-78.  LOGACS  Density  Data  for  Rev.  54,  Away  from  Perigee 
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Fig.  B-79.  LOGACS  Density  Data  for  Rev.  55,  Toward  Perige 
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Fig.  B-HO.  LOGAC'S  Density  Data  for  Rev.  55,  Away  from  Perigee 
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Fig.  B-82.  LOGACS  Density  Data  for  Rev.  56,  Away  from  Perigee 
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Fig.  R-84.  LOGACS  Density  Data  for  Rev.  57,  Away  from  Perigee 
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Fig.  B-85.  LOGACS  Density  Data  for  Rev.  58,  Away  from  Perig 
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Fig.  B-86,  LOGACS  Density  Data  for  Rev.  59,  Toward  Perigee 
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Fig.  B-87.  LOGACS  Density  Data  for  Rev.  59,  Away  from  Perigee 
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Fig.  B-88.  LOGACS  Density  Data  for  Rev.  61,  Away  from  Perig 


Fig.  B-89.  LOGACS  Density  Data  for  Rev.  62,  Toward  Perigee 
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Fig.  B-90.  LOGACS  Density  Data  for  Rev.  62,  Away  from  Perig 


Fig.  B-91.  LOGACS  Density  Data  for  Rev.  63,  Toward  Perigee 
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Fig.  B-92.  LOGACS  Density  Data  for  Rev.  63,  Away  from  Perigee 
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Fig.  B-93.  LOGACS  Density  Data  for  Rev.  64,  Toward  Perigee 
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LOGACS  Density  Data  for  Rev.  64,  Away  from  Perig 
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P’ig.  B-95.  LOGACS  Density  Data  for  Rev.  65,  Toward  Perigee 
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Fig.  B-96.  LOGACS  Density  Data  for  Rev.  65,  Away  from  Perigee 
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Fig.  B-97.  LOGACS  Density  Data  for  Rev.  66,  Toward  Perigee 
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Fig.  B-98.  LOGACS  Density  Data  for  Rev.  66,  Away  from  Perigee 
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APPENDIX  C 

LOGACS  WIND  DATA 
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Fig.  C-3.  Pitch  and  Yaw  Angle-of- Attack  Time  Histories 
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Fig.  C-3 .  Pitch  and  Yaw  Angle-of-Attack  Time  Histories  (Continued) 
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